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Block  20.  Abstract  (concluded) 

I")  The  determination  of  ingredient  percentages  with  solids  particle- 

sizes  held  constant. 
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a -2}  The  in-situ  determination  of  ammonium  perchlorate  and  aluminum 

particle  sizes  with  ingredient  percentages  held  constant.  ■ ' 
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jy  The  simultaneous  determination  of  ingredient  percentages  and 
solids  particle  sizes. 

The  methodology  was  developed  and  demonstrated  for  low-burning-rate  polybuta- 
diene acrylic  acid  type  propellants  and  high-burning-rate  hydroxyl-terminated 
polybutadiene  propellants  containing  ultrafine  ammonium  perchlorate.  Quanti- 
tative procedures  were  developed  for  determining  ferric  oxide,  polybutadiene 
acrylic  acid  polymer,  ammonium  perchlorate,  and  aluminum  in  polybutadiene 
acrylic  acid  propellants,  and  ammonium  perchlorate,  aluminum,  and  a ballistic 
modifier  in  high  rate  hydroxyl-terminated  polybutadiene  propellants.  The 
procedures  are  generally  applicable  to  all  types  of  composite  propellants. 
Emphasis  was  placed  on  the  eatabliahment  of  procedures  directly  applicable 
to  a propellant  manufacturing  process. 

Propellant  aamples  were  analyzed  nondeatruct ively  in  most  caaes  with  an 
estimated  relative  atandard  deviation  and  a relative  error  for  ingredient 
determinations  of  1Z  to  2Z.  Ultrafine  ammonium  perchlorate  agglomeration  in 
high-rate  propellants  reduced  the  precision  of  aluminum  determinations.  The 
estimated  relative  standard  deviation  in  this  case  was  4Z  to  5Z.  The  total 
analysis  time  for  four  replicates  of  a propellant  batch  was  15  to  30  min  with 
the  manual  instrumentation  uaed.  Multiple  linear  calibration  methods  were 
uaed  to  correct  for  matrix  effects.  Stable  reference  standards  were  used  to 
compensate  for  instrumental  fluctuations;  corrections  were  made  for  variable 
emission  line  absorption  by  the  Mylar  films  used  on  the  sample  holders. 
Statistical  procedures  were  developed  for  placing  Joint  confidence  intervals 
on  the  actual  ingredient  percentages  of  a production  propellant  batch. 
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INTRODUCTION 


Solid  composite  propellants  are  used  in  the  propulsion  systems 
of  many  types  of  Array  rcckets  and  missiles.  These  composite  propellants 
depending  on  the  particular  application,  are  composed  of  various  combi- 
nations of  a rubber-base  binder;  an  oxidizer  such  as  ammonium  perchlor- 
ate; a fuel  such  as  aluminum  powder;  a ballistic  modifier  such  as  ferric 
oxide  or  ferrocenes;  and  an  aliphatic  or  aromatic  ester  type  plasticizer 
The  propellant  mechanical  properties  are  controlled  primarily  by  the 
type  of  binder  system  used  and  by  the  binder-solids  interaction  char- 
acteristics. The  propellant  ballistic  and  rheological  properties  are 
strongly  affected  by  the  particle  sizea  of  the  solids  and  the  types  and 
percentages  of  the  ballistic  modifier  and  plasticizer.  The  propellant 
burning  rate  at  fixed  pressure  is  a particularly  significant  ballistic 
parameter  and  the  particle  size  of  the  ammonium  perchlorate  plays  a 
significant  role  in  rate  adjustment  and  control. 

Clearly,  both  the  propellant  ingredient  percentages  and  the  solid 
particle  sizes  must  be  carefully  controlled  during  propellant  manu- 
facture to  insure  that  the  finished  propellant  will  have  acceptable 
performance  and  reproducible  ballistic,  mechanical,  and  rheological 
properties.  Although  uncured  composite  propellants  can  be  analyzed 
by  a comb in at i on  of  existing  wet-chemical  and  instrumental  methods, 
these  methods  lack  the  speed  and  selectivity  required  for  routine 
quality  control  applications  in  propellant  manufacturing.  Moreover, 
existing  instrumental  methods  are  not  suitable  for  controlling  the 
particle  sizes  of  propellant  solids  after  they  are  incorporated  in  the 
propellant.  Cured  propellants  are  very  difficult  to  analyze  by 
wet-chemical  methods  because  of  the  intractable  nature  of  the  cured 
binder. 

X-ray  fluorescence  spectrometry  has  been  used  by  the  Army  Propul- 
sion Directorate  of  the  US  Array  Missile  Research  and  Development 
Command  {1)  for  many  years  in  propellant  research  applications.  E^rly 
applications  of  X-ray  fluorescence  spectrometry  to  composite  propellant 
analysis  were  also  reported  by  the  Thiokol  Corporation  (2].  The 
in-house  research  conducted  prior  to  initiation  of  this  project  demon- 
strated that  the  X-ray  fluorescence  method  can  be  advantageously  used 
as  a tool  to  monitor,  control,  and  improve  the  quality  of  production 
propeilanta.  X-ray  fluorescence  spectrometry  [3]  is  especially  attrac- 
tive for  propellant  analysis  because  of  ita  speed,  high  degree  of  pre- 
cision, and  the  fact  that  samples  can  l*e  analyzed  nondestruct ively 
without  prior  chemical  treatment.  Furthermore,  X-ray  fluorescence 
spectrometry  is  the  only  known  technique  that  is  capable  of  in-situ 
propellant-solids  particle  size  measurements  (1),  Because  of  these 
unique  features,  X-ray  fluorescence  analysis  of  uncured  production 
ptopeiianta  prior  to  motor  caating  enables  a decision  to  accept  or 
reject  the  batch  to  be  made  thereby  preventing  aubsequent  costiy  motor 
rejections,  if  unexplained  propellant  problems  arise  later,  the  cured 
propellant  can  be  analyzed  by  a similar  nondestructive  procedure. 
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This  project  was  conducted  as  part  of  the  US  Army  Materials  Testing 
Technology  Program,  The  objective  was  to  develop  a rapid,  precise,  and 
accurate  X-ray  fluorescence  method  of  analysis  for  general  application 
to  all  types  of  composite  propellants  used  in  Army  missile  systems. 
Emphasis  was  placed  on  the  development  of  techniques  directly  applicable 
to  propellant  manufacture.  The  method  was  specifically  applied  to 
polybutadiene  acrylic  acid  (PBAA) , low  burning  rate  propellants,  and 
hyciroxyl-terminated  polybutadiene  (HTPB)  high  burning  rate  propellants 
because  these  types  afforded  the  best  combinations  of  variables  needed 
to  develop  the  required  analytical  procedures.  The  developed  proce- 
dures can  be  readily  applied  with  little  or  no  modification  to  other 
types  of  composite  propellants.  Three  experimental  cases  were  con- 
sidered as  follows: 

a)  The  determination  of  propellant  ingredient  percentages  with 
solids  particle  sizes  held  constant. 

b)  Tlie  determination  of  ammonium  perchlorate  and  aluminum 
particle  sizes  with  ingredient  percentages  held  constant. 

c)  The  simultaneous  determination  of  ingredient  percentages 
and  particle  sizes. 

Appropriate  calibration  procedures  were  developed  to  handle  each  of 
these  experimental  cases.  Both  cured  and  uncured  propellants  were 
analyzed. 


If.  EXPERIMENTAL  METHODS 

A.  Instrumentation 

A wavelength-dispersive  universal  vacuum  X-ray  spec- 
trometer marketed  by  Philips  Electronic  instruments  was  used.  The 
fiat-crystal  X-ray  optical  system  of  the  spectrometer  is  shown  in 
Figure  1.  The  spectrometer  has  four  sample  compartments  each  of  which 
can  be  Individually  rotated  above  the  primary  X-ray  beam.  With  the 
Inverted  optical  system  the  bottom  surface  of  the  sample  it  irradiated. 
Either  a Philips  FAQ  60/1  (1600  W)  chromium  target  X-ray  tube,  or  a 
Philips  FAQ  60/i  (1900  W)  tungsten  target  X-ray  tube  was  used  depending 
on  the  analysis  requirements.  Both  X-ray  tubes  were  powered  by  a 
3-kVA  water-cooled  generator.  Tie  voltage  to  the  generator  was 
stabilized  with  a 5-kVA  line  voltage  stabilizer. 

Other  spectrometer  components  consisted  of  a 10.2  cm  * 0.S1  mm 
parallel  plate  entrance  collimator,  sodium  chloride  (200),  pentaeryth- 
rltol  (002),  and  echylenediamine  D-tartrate  (020)  analyzing  crystals, 
a coarse  exit  collimator,  and  a gas-flow  proportional  detector.  Hie 
associated  electronic  circuit  panel  (Type  12206/0)  has  a decade  scaler 
and  a single  channel  pulse  height  analyzer.  The  X-ray  optical  path 
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was  flushed  with  helium  for  uncured  propellant  analysis  and  evacuated 
(<  0.1  torr)  for  cured  propellant  analysis.  More  specific  instrumental 
operating  conditions  are  given  aiong  with  each  experiment  described 
later. 

The  choice  of  X-ray  tube  depends  on  the  elements  to  be  analyzed. 

A comparison  of  emission  line  intensities  for  propellant  elements  of 
interest  in  tills  investigation  using  tungsten  and  chromium  target  X-ray 
tuhes  is  given  in  Tahie  i.  The  chromium  target  tube  is  a better  choice 
for  the  analysis  of  sulfur,  chlorine,  and  aluminum;  therefore,  it  is  a 
more  generally  useful  tube  for  light-element  propellant  analyses. 


TABLE  i.  COMPARISON  OF  MISSION  LINE  INTENSITIES  FOR  PROPELLANT 
ELEMENTS  USING  TUNGSTEN  AND  CHROMIUM  TARGET  X-RAY  TUBES 


Emission 

Line 

Intensity 

(counts/sec) 

Rat  io 
(Cr/W) 

Chromium  Tube 

Tungsten  Tube 

FeK 

Cl 

3418 

25,090 

0.14 

C1K 

u 

55,990 

25,390 

2.21 

SK 

t 

1823 

766 

2.  38 

AIK 

7435 

2926 

2.54 

Sample:  Uncured  propellant  reference  standard. 

Chromium  tube:  50  kV,  28  mA  constant  potential. 

Tungsten  tuhe:  50  kV,  50  mA  constant  potential. 

The  sodium  chloride  crystal  was  selected  primarily  for  the  deter- 
mination of  the  low  percentage  of  sulfur  in  the  PBAA  polymer.  Either 
the  ethy iened inmine  D-tartrate  (EDDT)  or  the  pentaerythr itol  (PET) 
crystal  must  be  used  for  the  aluminum  determination.  The  relative 

reflectivities  of  these  crystals  for  several  light  elements  K emission 

<» 

lines  are  shown  in  Tahlc  2.  The  PET  crystal  is  a better  choice  when 
maximum  emission  line  intensity  is  required.  The  interpianar  d-spacing 
of  the  PET  crystal  varies  with  temperature,  however,  so  that  the 
emission  Uae  angle  most  he  adjusted  to  detect  the  emission  line  peak 
as  the  crystal  t raperature  varies. 

Pulse  height  analysis  was  used  to  reduce  background  and  increase 

the  p*. ak-to-background  ratios  for  aluminum  and  sulfur  determinations. 

Pulse  height  discrimination  was  less  effective  for  sulfur  K measure- 

u 

menis,  however,  because  the  chiorine  K fluorescence  from  the  sodium 
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TABLE  2.  COMPARISON  OF  THE  REFLECTIVITIES  OF  EDDT  AND  PET  CRYSTALS 

FOR  THE  K EMISSION  LINES  OF  SEVERAL  LIGHT  ELEMENTS 

a 


Emission 

Line 

Element  or 
Compound 

Intensity 

(counts/sec) 

Ratio 
PET/ EDDT 

Tube* 

Setting 

EDDT 

PET 

AiK 

a 

Ai 

9035 

16,975 

1.88 

40  kV,  30  mA 

SiK 

a 

Si 

3986 

7197 

i .81 

40  kV,  30  mA 

FK 

a 

NaH.PO, ‘H.O 
l 4 1 

1024 

1673 

i . 63 

40  kV,  30  mA 

SK 

a 

s 

13,537 

20,109 

1.49 

35  kV,  30  mA 

*Cr  target  tube  operated  at  constant  potential. 


chloride  crystal  was  also  passed  by  the  pulse  height  analyzer,  in  both 
cases  the  pulse  height  analyzer  primarily  eliminated  scattered  short 
wavelength  radiation  from  the  X-ray  tube  continuum.  There  was  no 
problem  with  spectral  line  interference.  A typical  pulse  height  dis- 
tribution curve  for  aluminum  K radiation  measured  with  the  gas-flow 

a 

proportional  detector  is  shown  in  Figure  2.  The  effect  of  Increasing 
the  gas-flow  detector  gain  (voltage)  on  Intensities  of  the  analytical 
emission  lines  is  shown  in  Figure  3.  The  detector  voltage  was 
operated  in  the  plateau  region  for  each  element. 

Only  peak  X-ray  intensity  measurements  were  made;  that  is,  no 
correction  was  made  for  background  radiation.  The  use  of  X-ray 
Intensity  measurements  only  at  the  peaks  of  the  analytical  emission 
lines  is  justified  because  of  the  calibration  method  used  and  the  fact 
that  experimental  conditions  were  chosen  to  give  large  peak-to-background 
ratios.  Typical  peak-to-background  ratios  for  the  analysis  of  light 
elements  in  PBAA  propellants  are  given  in  Table  i.  The  peak  inten- 
sities were  generally  kept  beiow  23,000  counts/sec  because  at  higher 
Intensities  the  response  was  nonlinear  due  to  the  dead  time  of  the 
linear  amplifier  of  the  electronic  circuit  panel. 

B.  Propellant  Preparation 

Aii  propellants  were  made  In  500-g  batches  (0,47-ilter 
size)  in  a vertical  double-sigma  biade  Baker-Perkins  type  mixer.  The 
mixer  bowl  was  heated  to  a temperature  of  55*  to  60*C»  Then  the  liquid 
components,  except  for  the  binder  curing  agent,  were  added  and  blended 
together.  The  propellant  solids  were  then  added  incrementally.  After 
the  last  solids  addition,  the  propellant  was  mixed  for  2 hr.  Then  the 
curing  agent  was  added  and  mixing  was  continued  for  20  min  under  vacuum. 
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Figure  2.  Pul**  height  distribution  curv*  (differential) 

for  alualnun  K radiation  and  a gaa-flov  detector 
(1 

eat  at  1600  Vdc. 
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)r  plateaus  for  propellant  elements 


TABLE  3.  PEAK-TO-BACKGROUND  RATIOS  FOR  THE  ANALYSIS 
OF  LIGHT  ELEMENTS  IN  PBAA  PROPELLANTS 


Emission 

Line 

Uncured  Propellant 

Cured  Propellant 

Without  PHA 

With  PHA 

Without  PHA 

With  PHA 

AIK 

a 

6.0 

180.1 

i5 . 3 

275.5 

SK 

a 

5.9 

10.3 

5.8 

9.6 

ClK 

a 

176.5 

265.7 

290.0 

457.4 

FeK 

a 

56.0 

85.5 

47.2 

72.0 

The  resulting  uncured  propellant  slurry  having  a viscosity  of  1 to  15 
kllopolse  was  taken  from  the  mixer  and  analyzed  directly.  Cured 
propellant  samples  were  prepared  by  vacuum  casting  the  propellant 
slurry  into  an  appropriate  2.5-  to  C>.l-cm  diameter  container  ^uch  as 
a mailing  carton  or  Teflon  tube  and  then  curing  the  propellant  at 
55°  to  60°C  for  severe  1 days.  This  is  a generalized  procedure.  Some 
deviations  were  made  to  accommodate  specific  formulations. 

Production  propellants  are  made  in  generally  the  same  manner 
except  that  the  batch  sizes  are  much  larger.  Typical  production  mixer 
sizes  are  1135-  to  2271-liter  capacities. 

C.  Mylar  Film  Corrections 

Uncurec*  propellant  samples  are  analyzed  in  a circular 
aluminum  holder  fitted  with  a thin  Mylar  film.  The  propellant  is 
pressed  against  the  film  and  held  in  place  by  gravity.  As  shown  in 
Figure  1,  the  primary  X-rays  from  the  X-ray  tube  pass  through  the  Mylar 
film  and  excite  the  elements  In  the  propellant  surface  to  fluoresce. 

The  character i-.tic  X-ray  fluorescence  emission  lines  than  pass  through 
the  Mylar  film  where  they  are  dispersed  by  the  analyzing  crystal  and 
detected  by  the  gas-flow  proportional  counter.  The  Mylar  film  must  be 
strong  enough  to  support  the  propellant  sample,  but  thin  enough  to 
transmit  a high  percentage  of  the  Incident  fluorescent  radiation.  The 
transmittance  for  a given  film  thickness  varies  significantly  with  the 
wavelength  of  the  fluorescent  radiation,  becoming  smaller  as  the  wave- 
length Increases.  An  experimental  determination  of  the  transmittance 
of  several  emission  lines  of  interest  through  two  different  Mylar  film 
thicknesses  Is  shown  In  Table  4.  Aluminum  radiation  Is  strong) 

absorbed  even  by  the  3.8-pm  film;  whereas  iron  K radiation  Is  absorb*.  1 

very  little.  The  3.8-pra  Mylar  film  was  used  for  most  of  this  work 
because  of  Its  higher  transmittance. 
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TABLE  4.  TRANSMISSION  OF  SOFT  X-RAYS  BY  MYLAR  FILMS 


Emission 

Line 

Transmittance 

Wavelength,  A 

6.4-pm  Mylar 

3.8-pm  Mylar 

FeK 

a 

1.94 

0.98 

0.99 

C1K 

a 

4.73 

0.76 

0.85 

SK 

a 

5.37 

0.68 

0.79 

AIK 

a 

8.34 

0.24 

0.43 

The  absorption  of  the  emission  lines  of  chlorine,  sulfur,  and 

aluminum  by  the  film  would  have  little  effect  on  uncured  propellant 

analysis  if  the  film  thickness  remained  constant  from  sample  to  sample. 

Unfortunately,  the  thickness  varies  sufficiently  among  Mylar  samples  to 

introduce  significant  intensity  measurement  errors  if  a correction  for 

the  film  thickness  variation  is  not  made  for  chlorine  K , sulfur  K , 

a a 

and  aluminum  K intensity  measurements.  As  indicated  by  the  data  in 

a 

Table  4,  no  correction  need  be  made  for  the  iron  K intensity 
measurements. 

An  empirical  procedure  to  correct  for  the  effects  of  variable 
Mylar  film  thickness  on  the  measured  intensities  of  chlorine  K^, 

sulfur  K^,  and  aluminum  radiations  was  previously  reported  by  Ailey 

and  Higgins  (4).  The  same  procedure  was  used  here  except  that  the 
accuracy  of  film  thickness  corrections  was  improved  by  using  a better 
model  derived  from  X-ray  absorption  theory,  and  by  reanalyzing  standard 
samples  to  get  a better  fit  of  the  experimental  data  to  the  model.  In 
this  procedure  aluminum  correction  factors  are  determined  directly 

by  analyzing  identical  aluminum  standards,  and  the  chlorine  and 

sulfur  K correction  factors  are  calculated  using  the  aluminum  K data, 
a u 

The  equations  for  calculating  the  chlorine  K and  sulfur  K correction 

factors  are  as  follows: 


log  C - 0.2078  log  CA1  (1) 

log  Cs  - 0.3189  log  CA1  (2) 
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where'  C.^,  and  are  the  correction  factors  for  the  elements  indi- 
cated as  subscripts.  For  routine  appl!cntions  where  large  numbers  of 
samples  are  analyzed,  it  is  convenient  and  facilitates  the  analysis  to 
have  the  correction  factors  tabulated  as  shown  in  Table  5. 

An  example  of  the  procedure  for  determining  Mylar  film  thickness 

correction  factors  for  aluminum  K , sulfur  K , and  chlorine  K radia- 
nt a a 

tions  is  given  in  Table  6.  Application  of  the  factors  to  propellant 
data  will  be  illustrated  later.  In  practice,  identical  high  purity 
aluminum  standards  are  sequentially  analyzed  in  the  holder  that  will  be 
used  for  the  stable  reference  standard  and  in  each  sample  holder.  In 
this  example  a total  of  four  sampler,  of  the  same  propellant  batch  will 
be  analyzed  as  a replication  of  duplicates.  The  aluminum  intensity 

from  each  aluminum  standard  after  transmission  through  the  Mylar  film 
on  the  holder  was  measured  in  seconds  to  collect  500,000  total  counts. 
The  aluminum  correction  factor  for  each  sample  holder  is  ohtained 

by  dividing  the  seconds  for  the  reference  standard  holder  by  the 

seconds  for  the  particular  sample  holder.  The  correction  factors  for 

sulfur  K and  chlorine  K emission  lines  are  rhen  calculated  using 
a a 

Equations  (1)  and  (2),  or  they  are  obtained  from  the  tabulated  values 
in  Table  5. 


D.  Reference  Standards 

There  are  several  inherent  short-  and  long-term  sources 
of  variation  in  X-ray  spectrometry  that  affect  the  measured  intensities 
of  X-ray  fluorescent  emission  lines.  Some  of  these  are  as  follows: 

1)  Fluctuations  in  the  X-ray  tube  output. 

4 * 2)  Mechanical  errors  of  positioning  samples  and  the 

goniometer. 

3)  Electronic  drift. 

To  insure  the  highest  possihle  analytical  precision  and  accuracy,  some 
type  of  standard  must  be  employed  to  compensate  for  the  short-  and 
long-term  sources  of  variation.  The  preferred  approach  for  propellant 
analysis  (the  one  used  in  this  work)  is  to  analyze  a reference  standard 
In  conjunction  with  the  propellants.  The  reference  standard  must  he 
stable  chemically  and  toward  repeated  exposure  to  the  primary  X-ray 
beam.  It  must  be  affected  hy  the  sources  of  variation  in  essentially 
the  same  manner  as  the  propellant.  Another  Important  advantage  of 
using  a stahle  reference  standard  is  that  calibration  rr  otdures,  which 
are  somewhat  involved,  are  valid  over  long  periods  of  time. 

A substantial  amount  of  effort  was  devoted  to  the  development, 
preparation,  and  evaluation  of  suitable  reference  standards  for  use 
with  uncured  PBAA  and  HTPB  propellants,  and  cured  PBAA  propellants. 
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Sulfur  K 


.900  0.967  0.978 

.905  0.969  0.979 

.910  0.970  0.981 

.915  0.972  0.982 

.920  0.974  0.983 

.925  0.975  0.984 

0.930  0.977  0.985 

0.935  0.979  0.986 

0.940  0.980  0.987 

0.945  0.982  0.988 

0.950  0.984  0.989 

0.955  0.985  0.990 

0.960  0.987  0.992 

0.965  0.989  0.993 

0.970  0.990  0.994 

0.975  0.992  0.995 

0.980  0.994  0.996 

0.985  0.995  0.997 

0.990  0.997  0.998 

0.995  0.998  0.999 

1.000  1.000  1 . 000 

1.005  1.002  1.001 

1.010  1.003  1.002 

1.015  1.005  1.003 

1.020  1.006  1.004 

1.025  1.008  1.00S 

'1.030  1.009  1.006 

1.035  1.011  1.007 

1.040  1.013  1.008 

1.045  1.014  1.009 

1.050  1.016  1.010 

1.055  1.017  1.011 

1.060  1.019  1.012 

1.065  1.020  1.013 

1.070  1.022  1.014 

1.075  1 ?023  1.015 

1.080  1.025  1.016 

1.085  1.026  1.017 

1.090  1.028  1.018 

1.095  1.029  1.019 


Aluminum  K 

a 


Chlorine  K 


1.100  1.031  1.020 
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TABLE  6.  EXAMPLE  OF  MYLAR  FILM  THICKNESS  CORRECTION  PROCEDURE 


Samp le 
Holder 

Mylar  Aluminum 
Standard 

(sec/500,000  counts) 

Aluminum 

Correction 

Factor 

Sulfur 

Correction 

Factor 

Chlorine 

Correction 

Factor 

Reference 

Standard 

21.66 

Sample  1 

21.95 

0.987 

0.996 

0.997 

Sample  2 

22.21 

0.975 

0.992 

0.995 

R aferenoe 
Standard 

21.72 

Sample  3 

21.47 

1 .012 

o 

o 

*— < 

1.002 

Sample  4 

20.89 

1.040 

1.013 

1.008 

The  compositions  of  the  standards  that  were  developed  and  used  through- 
out the  program  are  shown  in  Table  7.  It  was  necessary  to  develop  a 
standard  for  each  propellant  type  and  physical  state  so  that  analytical 
emission  line  intensities  from  the  standards  reasonably  approximated 
the  line  intensities  of  the  same  elements  from  the  corresponding 
propellants.  In  the  case  of  hTPB  propellant  only  the  uncured  propellant 
was  analyzed  during  this  program. 
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TABLE  7.  REFERENCE  STANDARD  COMPOS iTIONS  (WEIGHT  7.) 


P3AA  Propellant 

Uncur'd 
High  R..te 
HTPB  Propellant 

Ingredient 

Uncured 

Cured 

Sodium  Chloride 

35.0 

54.3 

55.0 

Aluminum  Powder 
(10  to  35  pm) 

22.0 

14.0 

20.0 

a-Cel lulose 

41.0 

30.6 

20.0 

Iron  (II)  Oxide 

0.8 

0.7 

Zinc  Sulfide 

1.2 

0.7 

5.0 

Sodium  chloride  was  used  in  place  of  ammonium  perchlorate  in  the 
reference  standard  because  ammonium  perchlorate  decomposes  when  exposed 

i 

I 

A 
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to  the  primary  X-ray  beam  for  periods  in  excess  of  approximately 
30  min.  Zinc  sulfide  was  used  in  place  of  PBAA  polymer  to  provide 
the  sulfur  emission,  because  it  was  desirable  to  prepare  a solid 

standard,  a-cellulose  is  an  excellent  binding  material  composed  of 
light  elements,  and  it  is  stable  toward  extended  X-ray  exposure.  Each 
standard  was  made  by  blending  the  mixture  for  15  min  on  a SPEX  No.  8000 
mixer/mill,  and  then  pressing  the  hlended  mixture  into  a pellet  in  a 
3.18-cm  die  at  a pressure  of  207  MPa.  1'ie  pellet  was  then  bonded  to 
a plexiglas  disc  of  the  same  diameter  to  facilitate  handling.  When 
not  being  used,  the  standard  pellets  are  stored  in  a desiccator.  The 
standards  have  been  found  to  have  excellent  long  terra  storage  stability. 

The  ability  of  the  reference  standards  to  compensate  for  variations 
in  the  intensity  of  the  X-ray  tube  primary  radiation  is  demonstrated  for 
the  cured  PBAA  standard  in  Table  8.  This  is  a very  severe  test,  because 
the  X-ray  tube  kV  end  mA  settings  were  purposely  varied  over  a very  wide 
range.  In  actual  practice  the  variations  will  be  many  orders  of  magni- 
tude smaller.  The  constancy  of  the  X-ray  intensity  ratio  for  each 
element  in  Table  8 as  the  X-ray  tube  settings  were  changed  shows  the 
excellent  compensating  ability  of  the  standard.  The  other  standards 
are  equally  effective. 

The  stability  of  the  uncured  PBAA  propellant  standard  upon 
extended  X-ray  exposure  is  shown  in  Table  9.  Again,  the  essentially 
constant  X-ray  intensity  ratio  for  each  element  demonstrates  that  the 
constituents  of  the  standard  are  not  degraded  by  the  X-ray  exposure. 

The  use  of  compounds  in  the  reference  standard  for  stahility 
reasons  that  are  not  present  in  the  propellant  causes  an  undesirable 
experimental  effect  as  shown  In  Tahle  10.  It  Is  well  known  that  the 
exact  wavelength  and  hence  reflecting  angle  for  a given  emission  line 
depends  somewhat  on  the  electronic  environment  of  the  atom.  Thus  the 
peak  angles  for  <. hlorine  radiation  from  the  standard  and  propellant 

differ  by  0.10°  20.  Therefore,  for  the  most  precise  work  the  goniometer 
should  be  set  to  the  peak  angle  of  both  the  standard  and  propellant 
during  quantitative  ammonium  perchlorate  determinations.  The  difference 
hetween  the  peak  angles  of  sulfur  radiation  from  the  standard  and 

propellant  was  found  to  he  insignificant. 

E.  Spectrometer  and  Propellant  Variables  Evaluation 

Prior  to  the  development  and  application  of  calibration 
procedures  for  quantitative  determinations  and  the  analyses  of  large 
numbers  of  samples,  an  extensive  evaluation  of  propellant  mixing  and 
spectrometer  operating  variables  was  made  to  establish  optimum  condi- 
tions, The  main  variables  that  were  evaluated  are  shown  In  Table  11. 
These  are  the  ones  that  were  considered  most  likely  to  affect  the 
precision  and  accuracy  of  quantitative  determinations. 
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TABLE  9.  STABILITY  OF  UNCURED  PBAA  PROPELLANT 
REFERENCE  STANDARD  TOWARD  X-RAYS* 


Cumulat ive 
Time  (hr) 

Intensity  Ratio 

Iron  Ka 

Chlorine  K 

c X 

Sulfur  K 

a 

Aluminum  Ka 

0 

1.147 

1.230 

0.858 

1.577 

1.52 

1.149 

1.2  34 

0.852 

1.578 

2.02 

1.149 

1.234 

0.852 

1.567 

2.5!i 

1.146 

1.236 

0.856 

1.569 

3.38 

1.149 

1.238 

0.855 

1.579 

3.97 

1.146 

1.239 

0.858 

1.577 

4.57 

1.145 

1.243 

0.861 

1.587 

5.17 

1.146 

1.243 

0.862 

1.581 

5.67 

1.145 

1.243 

0.867 

1.573 

6.17 



1.148 

1.242 

0.861 

1.581 

♦Tungsten  target:  50  kV,  45  mA. 


TABLE  10.  EFFECT  OF  ELEMENT  ELECTRONIC  ENVIRONMENT  ON  MEASURED 
ANGLES  FOR  SELECTED  EMISSION  LINES 


Sample 

Analyzing 

Crystal 

Emission 

Line 

Goniometer 
Angle  (deg  20) 

A 20 

Sodium  Chloride 

NaCl 

Chlorine  K 

114.05 

0.10 

Ammonium  Perchlorate 

a 

113.95 

Potassium  Sulfate 

NaCl 

Sulfur  K 

144.57 

0. 18 

PBAA  Polymer 

ct 

144.75 

Aluminum  Oxide 

EDDT 

Aluminum  K 

142.68 

0.09 

Aluminum 

Cl 

142.77 
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The  significance  of  the  selected  variables  was  determined  by 
statistical  analysis  of  data  obtained  from  factorial  experimental  designs. 
Although  statistical  significance  between  the  levels  was  found  for 
several  of  the  variables,  in  most  cases  the  difference  was  of  no  practical 
significance.  There  was  no  practical  difference,  for  example,  between 
the  levels  of  the  mixing  variables.  Analysis  under  vacuum  has  an  adverse 
effect  on  uncured  propellant  results;  whereas  cured  propellant  can  be 
satisfactorily  analyzed  under  vacuum.  There  is  a small  effect  of  sample 
temperature  on  analysis  results;  it  is  preferable  to  allow  the  sample 
to  cool  to  nearly  ambient  temperature  before  analysis.  The  time  at 
which  the  sample  is  analyzed  after  it  has  been  loaded  into  the  sample 
holder  was  not  found  to  affect  results  significantly;  however,  it  is  good 
analytical  procedure  to  analyze  the  samples  as  soon  as  possible  after 
they  are  prepared.  Rotation  of  the  sample  in  its  own  plane  had  little 
effect  on  the  precision  of  uncured  propellant  analyses,  but  did  improve 
the  precision  of  some  cured  propellant  analyses,  particularly  if  the 
propellant  tended  to  be  inhomogeneous.  Whether  the  propellant  had  been 
deaerated  before  analysis  had  no  effect  on  results,  but  the  use  of 
deaerated  samples  is  nevertheless  a better  choice, 

F.  Analysis  Procedure 

i.  Uncured  Propellant 

The  sampling  components  used  for  the  majority  of 
the  uncured  propellant  analyses  are  shown  in  Figure  4(a).  A 3.8-pm 
Mylar  film  was  placed  on  the  bottom  of  each  circular  aluminum  holder 
to  support  the  reference  standard  and  propellant  slurry  samples.  The 
Mylar  film  was  then  supported  with  a metal  disc.  The  propellant  sample 
was  pressed  against  the  Mylar  surface  using  a plexiglas  backing  disc. 
Before  the  samples  were  loaded,  however,  the  Mylar  film  thickness 
correction  factors  were  established  as  shown  in  Table  6.  One  of  the 
aluminum  sta:iards,  machined  from  a piece  of  bar  stock,  is  shown  in 
Figure  4(a).  Properly  mounted  aluminum  foil  is  also  suitable  for  use 
ns  a standard  for  film  thickness  corrections. 

Although  the  aluminum  holders  work  very  well  for  uncured  propellant 
analysis,  they  must  be  cleaned  after  each  analysis.  This  is  time 
consuming  and  therefore  not  very  desirable  for  rapid,  routine  quality 
control  analysis  of  production  propellants.  Consequently,  toward  the 
end  of  the  program  disposable  Chemplex  No,  1530  (Chemplex  industries, 

Inc.)  plastic  sample  cups  were  purchased  and  evaluated.  These  sampling 
components  are  shown  in  Figure  4(b).  Mylar  film  corrections  and  sample 
loadings  were  accomplished  in  the  same  manner  as  for  the  a'uminum  holders. 
The  propellant  samples  in  this  case  were  pressed  against  the  3.8-pm 
Mylar  surface  using  the  large  end  of  a No.  10  cork  stopper.  The  loaded 
sample  cups  are  supported  by  the  aluminum  holders  during  analysis,  and 
are  subsequently  discarded.  Results  using  the  disposable  sample  cups 
are  comparable  to  those  using  the  aluminum  holders;  therefore,  they  are 
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(1)  ALUMINUM  STANDARD  FOR  MYLAR  CORRECTION 

(2)  PLEXIGLAS  DISC  FOR  PRESSING  PROPELLANT  INTO 
INTO  SAMPLE  HOLDER 

(3)  REFERENCE  STANDARD  MOUNTED  ON  PLEXIGLAS  DISC 
<41  CIRCULAR  ALUMINUM  SAMPLE  HOLDER 

(S)  UNCURED  PROPELLANT  AND  CONTAINER 
(S)  PLATF  AND  METAL  DISC  TO  SUPPORT  MYLAR  FILM 
ON  HOLDER  DURING  SAMPLE  LOADING 

Figure  4 (a).  Components  for  uncured  propellant  analysis 
oy  X-ray  spectrometry. 
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(1)  ALUMINUM  STANDARD.  (2)  NO.  10  CORK  STOPPER, 

<31  C HEMP  LEX  NO.  1530  DISPOSABLE  HOLDER  A NO 
SNAP-ON  RING.  <41  ASSEMBLED  NO.  1530  HOLDER. 

<5)  NO.  1530  HOLDER  LOADED  WITH  PROPELLANT 
SAMPLE.  <S)  PROPELLANT  SURFACE  IN  HOLDER. 

AND  (71  CIRCULAR  ALUMINUM  SAMPLE  HOLDER. 

Figure  4 (b).  Components  for  uncured  propellant 
analysis  by  X-ray  spectrometry. 
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recommended  for  the  analysis  of  uncured  propellant.  A certain  amount  of 
art,  however.  Is  required  to  install  the  Mylar  film  on  the  cup  properly. 

Samples  were  analyzed  in  a helium  environment  either  in  duplicate 
or  triplicate  in  conjunction  with  the  appropriate  reference  standards 
shown  in  Table  7.  Because  only  peak  X-ray  intensity  measurements  were 
made  for  each  emission  line,  either  a fixed  count  or  fixed  time  measure- 
ment technique  was  appropriate.  The  fixed  count  technique  was  used  for 
the  majority  of  measurements  in  this  program  so  that  the  counting  error 
variance  component  could  be  readily  separated  from  the  other  sources  of 
error.  The  numbe>  of  seconds  required  to  collect  a preselected  fixed 

count  for  the  reference  standard,  t , and  slurry  samples,  t , were 

s c 

measured  in  rapid  succession  at  the  peak  analytical  goniometer  angle 
for  each  ingredient.  Average  X-ray  intensity  ratios  for  each  propellant 
batch,  in  this  case  for  a single  replication  of  duplicate  samples,  were 
calculated  by 


where  is  the  X-ray  intensity  ratio  averaged  over  samples  for  the  ith 

Ingredient  in  the  jth  propellant  botch.  An  example  of  a typical  uncured 
propellant  analysis  illustrating  the  calculation  procedure  and  applica- 
tion of  the  Mylar  correction  factors  is  shown  in  Table  12. 

2.  Cured  Propellant 

Cured  propellunt  samples  were  analyzed  in  a manner 
similar  to  that  described  for  the  uncured  samples  except  that  Mylar  film 
was  not  used,  and  as  mentioned  earlier,  the  samples  could  be  analysed  in 
vl ther  helium  or  vacuum.  Neither  the  cured  nor  uncured  samples  can  be 
analyzed  In  air  because  it  strongly  absorbs  the  long  wavelength  fluores- 
cent emission  lines  from  the  propellant  light  elements,  particularly 
the  aluminum  K line. 

•i 

Some  of  the  sample  preparation  and  handling  components  for  cured 
propellant  analysis  are  shown  In  Figure  5.  The  propellant  was  Initially 
cured  in  either  wax-coated  mailing  cartons  or  Teflon  tubing.  Afi.er 
curing,  the  container  material  was  removed  and  the  propellant  was  slieed 
into  O.feWra  thick  discs  using  a guillotine  built  especially  for  this 
purpose.  A microtome  blade  was  used  for  cutting;  it  resulted  in  very 
smooth  propellant  surfaces  without  pulling  large  ammonium  perchlorate 
particles  from  the  binder.  After  slicing  the  propellant,  discs  of 
3.18-cm  diameter  were  punched  out  to  fit  the  circular  holder  for 
analysis.  The  surfaces  of  the  propellant  discs  were  analyzed  as 
described  for  the  uncured  propellant  samples, 
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TABLE  12.  EXAMPLE  OF  UNCURED  PROPELLANT  (PBAA)  ANALYSIS  DATA  TREATMENT 


(1)  METAL  PUNCH  FOR  CUTTING  3.2-cm  PROPELLANT 
DISC.  (2)  GUILLOTINE.  (31  PROPELLANT  SAMPLE. 

(4)  PLEXIGLAS  BACKING  DISC.  (SI  CIRCULAR  ALUMINUM 
SAMPLE  HOLDER.  AND  (6)  MAILING  CARTON  IN  WHICH 
PROPELLANT  IS  CURED. 


Figure  5.  Components  for  cured  propellant  analysis 
by  X-ray  spectrometry. 


III.  CALIBRATION  METHODOLOGY 


A.  General  Considerations 

X-ray  spectrometry,  like  many  other  physical  methods,  is 
not  absolute.  That  is,  the  concentrations  and  particle  sizes  of  ingre- 
dients in  unknown  propellant  samples  must  be  determined  with  the  aid  of 
calibration  curves  or  equations  which  relate  elemental  emission  line 
intensities  with  sizes  and  concentrations.  Although  very  precise  measure- 
ments of  X-ray  emission  line  intensities  can  be  made  by  the  X-ray  fluo- 
rescence analysis  procedure,  accurate  elemental  determinations  are  often 
difficult  to  achieve  because  of  sample  matrix  or  interelement  effects. 
Thus,  simple  linear  calibration  procedures  that  suffice  for  a number  of 
other  types  of  spectrometry  are  generally  unsuitable  for  the  X-ray 
spectrometric  analysis  of  complex  mixtures  such  as  propellants.  Mitchell 
[5]  has  described  the  qualitative  and  quantitative  aspects  of  the  effects 
and  the  associated  problems. 

Stated  simply,  a significant  sample  matrix  effect  means  that  the 
intensity  of  characteristic  radiation  from  an  element  in  the  propellant 
depends  not  only  on  its  concentration  and  particle  size  but  also  on  the 
concentrations  and  particle  sizes  of  all  other  elements  in  the  sample. 

The  magnitude  of  these  matrix  effects  can  vary  markedly  depending  on  the 
type  end  relative  quantities  of  elements  in  the  sample.  A suitable  cali- 
bration procedure  must  accurately  compensate  for  these  matrix  effects. 
Three  main  experimental  techniques  have  been  used  in  X-ray  spectrometry 
to  minimize  or  correct  for  matrix  effects.  These  are  as  follows: 

1)  Addition  of  one  or  more  internal  standards  to  the  sample. 

2)  Dilution  of  the  sample  with  an  inert  component. 

3)  Restriction  of  ingredient  calibration  concentrations  and 
particle  sizes  to  narrow  ranges  coupled  with  simple  linear  calibration. 

The  first  two  methods  are  impractical  for  the  rapid,  accurate  analysis 
of  production  propellants;  the  last  method  has  limited  possibilities. 

The  most  practical  calibration  approach  for  propellant  analysis 
is  to  use  mathematical  procedures,  based  on  calibration  mixtures,  to 
compensate  for  matrix  effects  and  to  estimate  unknown  ingredient  per- 
centages. Many  different  theoretical  and  empirical  calibration  proce- 
dures have  been  developed  in  the  past  for  application  to  the  X-ray 
fluorescence  analysis  problem.  Considerable  experience  in  this  labora- 
tory with  mathematical  methods  indicates  that  multiple  regression 
analysis  is  the  best  approach  to  use  for  propellant  applications.  The 
application  of  multiple  regression  methods  to  propellant  analysis  has 
been  described  by  Alley  and  Myers  |6J.  Williams  (7)  gives  a discussion 
of  the  general  problem.  Simple,  and  multiple  linear  regression  methods 
were  considered  in  this  project.  Calibration  by  multiple  linear 
regression  analysis  is  the  better  approach  for  production  propellant 
analysis  applications. 
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B.  Regression  Methods 


1.  Simple  Linear  Regression 

The  simplest  calibration  method  for  X-ray  fluores- 
cence analysis  consists  of  establishing  a relationship  for  the 
regression  of  X-ray  intensity,  R,  from  an  element  on  its  percentage, 

X,  in  a sample  mixture;  namely, 


R = a + bX 


(4) 


As  already  stated,  this  model  is  of  very  limited  use  in  direct  propellant 
analysis  because  of  its  failure  to  account  for  interelement  or  matrix 
effects.  When  it  is  applicable,  the  ingredient  percentage  in  an  unknown 
sample  is  estimated  by  inverting  Equation  (4): 


y - (R  ~ a> 

X b 


(5) 


Propellant  analysis  results  using  Equation  (5)  are  valid  for  only  a very 
narrow  region  about  the  calibration  compositions  used  to  determine  the 
intercept,  a,  and  slope,  b.  If  the  slope  is  zero  or  of  negligible 
magnitude  and  the  slopes  for  both  the  calibration  and  unknown  regression 
lines  are  assumed  to  be  equal,  then  the  following  expression  can  be 
derived : 


X 

u 


X 

s 


R 

u 


R 

s 


(6) 


where  the  subscripts  u and  s refer  to  the  unknown  and  calibration 
propellants,  respectively.  This  expression  has  potential  use  in 
production  propellant  analysis  as  a "go-no-go"  test.  Because  the 
production  propellant  nominal  composition  is  known,  its  actual  composi- 
tion can  be  verified  by  comparison  against  a single  nominal  calibration 
composition  previously  analyzed  to  establish  Rg.  if  R^  agrees  with  Rg 

within  an  acceptable  confidence  region  for  each  Ingredient,  then  the 
production  batch  is  accepted.  On  the  other  hand,  if  the  actual 
propellant  composition  differs  significantly  from  nominal,  the  deter- 
mination of  X will  be  inaccurate  because  of  uncompensated  matrix 
,,  u 
effects . 


2.  Multiple  Linear  Regression 

Assuming  that  the  functional  relationship  between 
the  dependent  and  independent  variables  is  linear  over  the  percentage 
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ranges  of  Che  calibration  ingredients,  the  statistical  model  for  the 
analyses  of  four  ingredients  in  a mixture  when  particle  size  is  held 
constant  is: 


Rij  " Bi0  + Bil  Xlj  + Bi2  X2j  + Bi3  X3j  + Bi4  X4j  + Eij 


(7) 


(i  = 1,2, 3, 4) 

(j  * 1,2, ... ,n) 


where  R^.  is  the  intensity  ratio  for  the  itu  ingredient;  is  the 

percentage  by  weight  of  the  kth  ingredient  in  the  jth  mixture;  the  B^ 

are  the  regression  coefficients;  and  is  the  random  error  associated 

with  R^.  The  percentage  of  each  ingredient  appears  in  the  model  no 

matter  which  one  is  being  determined.  This  accounts  for  mutual  absorp- 
tion and  enhancement  effects  among  the  ingredients.  The  coefficients 
indicate  the  relative  amounts  of  radiation  absorption  and  enhancement 
occurring.  A negative  coefficient  indicates  absorption;  a positive 
coefficient  indicates  enhancement. 


Equation  (7)  can  be  used  to  develop  a set  of  working  expressions 
for  estimating  the  ingredient  concentrations.  Equation  (7)  in  matrix 
notation  is: 


R - b + B X 


(8) 


I 


where  R represents  the  vector  of  intensity  ratios,  b the  vector  of 
intercept  terms,  B the  matrix  of  regression  coefficients,  and  X the 
vector  of  unknown  ingredient  concentrations  to  be  determined  by 
analysis.  The  B^  element  of  B is  the  coefficient  of  X^  in  the  ith 

regression  equation.  inverting  Equation  (8)  to  solve  for  X gives: 


X - B_i  (R  - h) 


(9) 


> 

s 


f' 


assuming  that  B 1 exists.  Equation  (9)  was  then  used  to  estimate  the 
X^’s  from  X-ray  intensity  ratios,  as  calculated  from  Equation  (3),  with 

particle  sizes  held  constant.  Estimates  of  ammonium  perchlorate  and 
aluminum  particle  sizes  at  constant  ingredient  concentrations  were  also 
ohtained  with  Equation  (9)  by  replacing  the  £ with  the  particle  size 
vaiu.?,  W.  Equation  (9)  gives  more  accurate  analyses  than  Equation  (5), 
because  it  contains  terms  that  correct  for  matrix  effects.  Because  of 
the  assumed  linear  relationship,  Equation  (9)  is  valid  only  when 
propellant  ingredient  percentages  vary  o«r  reasonably  narrow  ranges. 
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This  is  perfectly  acceptable  for  the  analysis  of  production  propellants. 
The  actual  relationship  between  radiation  intensities  and  ingredient 
percentages  for  propellants  is  curvilinear,  but  this  can  be  approximated 
by  a linear  relationship  over  limited  concentration  regions. 

The  b.,  coefficients  in  Equation  (8)  were  derived  by  a least 
lk 

squares  analysis  which  involves  minimizing  the  sums  of  squares  of  the 
residual  errors.  For  this  analysis,  the  ingredient  percentages,  X^, 

were  chosen  as  the  independent  or  concomitant  variables  because  they 
can  be  controlled  and  their  measurement  error  is  small  in  comparison 
to  that  in  the  R’s.  From  an  experimental  standpoint  it  would  be 
simpler  to  reverse  the  role  of  the  X's  and  R's  to  estimate  the  X's 
directly.  However,  because  tbe  intensities  are  not  controlled  and  are 
measured  with  non-negligible  error,  this  latter  approach  cannot  be  used 
without  some  danger.  The  general  problem  has  been  discussed  by 
Berkson  [81.  The  inverse  of  the  model  of  Equation  (8)  was  used  in  a few 
cases  during  this  program,  but  the  pitfalls  of  using  such  a model  are 
fully  recognized. 

The  X-ray  fluorescence  method  can  also  be  used  to  determine  ingre- 
dient concentrations  and  particle  sizes  when  both  are  varied 
simultaneously  (9).  The  analysis  is  restricted,  however,  to  the  deter- 
mination ot  a number  of  parameters  equal  to  the  number  of  X-ray 
intensity  measurements  made.  The  particie  sizes  of  ammonium  perchlorate, 
W.,,  and  aluminum,  W^,  were  considered  here.  They  were  determined  as 

the  weight  percentages  of  a fin**  fraction  in  a bimodal  blend  of  fine 
and  coarse  fractions.  As  reported  earlier,  the  actual  average  (weight 
mean)  particle  size  in  the  propellant  can  be  determined  by  referring 
to  an  appropriate  calibration  curve.  The  main  objective  in  this 
project  was  to  determine  particle  size  changes  quantitatively  that 
could  affect  the  properties  of  production  propellants. 

Consider  a set  of  multiple  regression  equations  of  the  following 


Rlj  " 2 (Bik'V  * "l5U2J  + V" 


M + Cij 


(i  - 1,2, 3, 4) 


where  ■ 1.  Equation  (10)  contains  both  particle  size  and  concentra- 
tion terras,  and  can  be  written  in  matrix  notation  as 


R - + B2  W 
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The  intensity  ratio  vector  can  now  be  corrected  for  either  particle  size 
or  concentration  to  solve  for  one  effect  in  the  presence  of  the  other, 
cr  to  give  a combined  parameter  determination  equal  to  the  number  of 
intensity  measurements.  A complete  X-ray  analysis  will  require  partial 
information  about  the  propellant  composition  from  an  external  source. 
Solving  Equation  (11)  for  concentration,  X,  gives 

X = B"1  (R  - B2  W)  , (12) 

and  solving  for  particle  size,  W,  gives 

W = b'1  (R  - B X)  . (13) 


As  mentioned,  the  particle  size,  , is  expressed  as  the  weight  fraction 
of  fine  component  i in  a bimodal  size  mixture. 

C.  Experimental  Designs 

The  selection  or  preparation  of  calibration  mixtures  for 
the  least  squares  estimation  of  the  b^  coefficients  in  Equations  (8) 

and  (10)  is  critical  to  insure  that  the  coefficients  are  accurately 
determined  |6).  Otherwise,  the  analyses  of  unknown  propellant  composi- 
tions might  be  very  inaccurate.  In  addition,  it  is  desirable  to  use 
as  small  a number  of  calibration  mixtures  as  practicable  because  of  the 
amount  of  labor  involved  in  preparing  and  analyzing  them.  The  standards, 
in  any  event,  must  be  representative  of  the  type  of  propellant  that  will 
be  analyzed  with  the  resulting  calibration  equations. 

The  regression  coefficients  can  be  accurately  and  efficiently 
estimated  if  the  calibration  mixtures  are  prepared  in  accordance  with 
a suitable  statistical  experimental  design.  This  generaiiy  eliminates 
the  problem  of  confounding  of  effects  and  high  degrees  of  correlation 
among  ingredient  percentages  which  obviously  must  be  avoided.  Over  the 
years  this  laboratory  has  evaluated  many  different  types  of  experimental 
designs  for  application  to  the  propellant  analysis  calibration  problem. 
Some  of  these  are  factorial  and  fractional  factorial  designs  f 10 ) , 
central  composite  designs  [11),  simplex  lattice  [12)  and  simplex  lattice 
designs  with  reference  mixtures  [13,14],  extreme  vertices  designs  [15), 
and  rotatable  designs  [16)  recently  developed  that  are  especially 
applicable  to  mixture  problems  similar  to  the  propellant  analysis  case. 

The  factorial  and  fractional  factorial  designs  and  the  central 
composite  design  were  used  during  this  program,  but  in  some  cases 
another  design  might  work  as  well,  or  better.  The  appropriate  design 
must  be  selected  by  the  experimenter  depending  on  the  calibration  model 
and  the  objectives  of  the  experiment. 
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IV.  DETERMINATION  OF  INGREDIENT  CONCENTRATIONS 
WITH  PARTICLE  SIZE  CONSTANT 

A.  Low  Rate  PBAA  Propellant 

1.  Experimental 

4 

A half  fraction  of  a 2 factorial  was  used  to  prepare 
calibration  standards  for  the  model  represented  by  Equation  (7).  The 
factors  and  factor  levels  are  shown  In  Table  13,  and  the  design  data  and 
defining  contrasts  are  shown  in  Table  14.  With  this  design,  It  was 
necessary  to  use  the  weights  of  ingredients  as  factor  levels  instead 
of  their  percentages.  The  factor  levels  were  arbitrarily  chosen  to 
result  in  practical  ingredient  percentage  ranges  that  were  still  narrow 
enough  to  permit  a good  fit  of  the  linear  model  to  the  data. 


TABLE  13.  FACTORS  AND  FACTOR  LEVELS  FOR  CALIBRATION  PROPELLANT  BATCHES 
(CONSTANT  PARTICLE  SIZE) 


Factor 

Symbol 

Low  Level  (g) 

High  Level  (g) 

Ferric  Oxide 

A 

4.5 

5.5 

Ammonium  Perchlorate 

B 

660.0 

700.0 

PBAA  Polymer 

C 

125.0 

145.0 

Aluminum 

D 

150.0 

170.0 

The  treatment  combinations  in  Table  14  for  the  first  eight  calibra- 
tion mixtures  form  the  principal  block  of  the  full  factorial.  Following 
convention,  the  high  level  of  each  factor  In  the  treatment  combination 
is  denoted  by  the  presence  of  the  lower  case  factor  symbol,  and  the  low 
level  by  absence  of  the  symbol.  When  a specific  propellant  composition 
is  to  be  analyzed,  for  example  a production  propellant,  its  nominal 
composition  should  be  placed  at  the  midpoint  of  the  design.  This  Improves 
the  precision  of  estimating  the  propellant  ingredient  percentage.  The 
midpoint  composition  (Batch  9)  was  the  one  of  primary  Interest  in  this 
example.  In  addition  to  the  principal  block  and  midpoint,  three  more 
design  points  represented  by  Batches  10,  11,  and  12  were  added  to 
Increase  the  degrees-of-freedom  for  estimating  error. 

The  actual  percentage  compositions  of  the  12  calibration  batches 
that  were  made  in  accordance  with  the  design  and  tested  are  given  In 
Table  IS.  Each  batch  also  contained  a binder  component  that  could  not 
be  analyzed  and  which  is  therefore  not  shown  in  the  table.  As  a con- 
sequence, the  Ingredient  percentages  for  each  batch  are  Independent  of 
each  other,  which  is  a requirement  for  the  model  used.  As  shown  in 
Table  15,  the  blmodal  ammonium  perchlorate  blend,  and  hence  its  particle 
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TABLE  14.  DESIGN  DATA  AND  DEFINING  CONTRASTS  FOR  DETERMINING 
INGREDIENT  CONCENTRATIONS  IN  PBAA  PROPELLANT  (PARTICLE 
SIZES  CONSTANT) 


blend  of  40X  15-pm  ammonium  perchlorate  and  60%  200-ur,i 
i perchlorate. 


size,  was  held  constant  for  all  batches.  Tlu*  particle  size  of  the 
aluminum  powder  was  likewise  held  constant.  Polymers  cannot  normally 
be  analyzed  by  X-ray  spectrometry.  It  was  possible  to  analyze  the 
P BAA  polymer  here  because  ic  contains  a small  percentage  of  sulfur 
(approximately  1%)  in  the  polymer  chain.  Hc.?ver,  a recalibration  for 
P BAA  polymer  must  be  made  if  the  polymer  lot  Is  changed. 

The  analytical  emission  lines  and  instrumental  parameters  used  are 
given  in  Table  lb.  Peak  X-ray  emission  line  intensity  measurements 
were  made  by  a fixed  count  technique,  and  pulse  height  discrimination 
was  used  to  increase  the  peak-to-background  ratios  for  PBAA  polymer  and 
aluminum  analyses. 

2.  Uncured  Propellant  Results 

Four  samples  of  each  calibration  propellant  batch 
were  analyzed.  Duplicate  samples  were  analyzed  In  conjunction  with  the 
uncured  propellant  reference  standard,  and  this  analysis  was  repeated. 

The  individual  intensity  measurements  for  the  four  ingredients  from 
each  of  the  12  calibration  batches  are  recorded  in  Table  17  as  the 
seconds  required  to  collect  the  fixed  counts  listed  in  Table  16.  These 
raw  data  were  used  to  calculate  the  individual  Intensity  ratios  shown 
In  Table  18.  The  mean  X-ray  intensity  ratios  responses,  shown  in 
Table  19,  were  used  along  with  the  ingredient  percentages  in  Table  15 
to  derive  the  partial  regression  coefficients  for  the  simple  and  multiple 
linear  regression  models  used.  The  regression  analyses  were  made  by 
computer  (CD06600)  using  a stepwise  muitiple  regression  program  (17). 

The  regression  data  obtained  by  using  the  simple  regression  model 
of  Equation  (4)  are  given  In  Table  20.  In  addition  to  the  least  squares 
estimates  of  the  coefficients,  the  tahie  lists  the  estimated  standard 
error  of  the*£Bfef f icient , S^;  Student's  "t"  for  evaluating  the  signifi- 
cance of  the  coefficient;  the  standard  error,  S , for  estimating  R. ; 

2 ° 1 

and  the  correlation  index,  R^.  The  correlation  index,  which  should  not 

be  confused  with  the  X-ray  intensity  ratio,  is  the  fraction  of  the  total 
corrected  sum  of  squares  among  calibration  batches  that  is  explained  by 
regression.  it  measures  how  well  the  model  fits  the  data,  and  has  a 
value  of  one  for  a perfect  fit. 

The  fit  of  the  simple  linear  model  is  best  for  the  ferric  oxiue 
determination.  Because  iron  is  relatively  heavy  compared  to  the  otb* r 
elements  analyzed,  it  is  not  affected  as  much  by  matrix  effects.  The 
simple  regression  model  is  unsuitable  for  ammonium  perdue  rate  and  PBAA 
polymer  analyses  because  of  significant  matrix  effects  :hat  nic  no’, 
compensated  for  by  the  model.  The  ferric  oxide  results  for  the  simple 
regression  model  are  plotted  in  Figure  6 to  show  the  scatter  cf  data 
points  about  the  regression  line.  Figure  7 compares  regression  lines 
using  the  models  of  Equations  (4)  and  (6).  Analyses  using  F.quation  (6) 
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TABLE  17.  X-RAY  INTENSITY  MEASUREMENTS  IN  SECONDS 
FOR  FIXED  COUNTS*  FOR  UNCURED  PBAA  PROPELLANT 
ANALYSES  (CONSTANT  PARTICLE  SIZE) 


Ferric 

Ammonium 

PBAA 

Batch 

Oxide 

Perchlorate 

Polymer 

Aluminum 

mm 

24.45 

Wmm 

24.84 

1 

mm 

24.79 

17.21 

12.31 

HH 

24.50 

15.43 

13.45 

20.62 

24.07 

15.38 

13.58 

20.59 

25.16 

17.19 

12.22 

16.99 

22.78 

18.56 

13.73 

18.41 

22.93 

18.14 

13.81 

18.03 

22.81 

17.17 

12.21 

17.03 

23.05 

18.71 

13.84 

18.34 

23.27 

18.61 

13.70 

18.31 

23.17 

BK 1 

H|i| 

17.20 

22.71 

mam 

22.02 

20.17 

■SB® 

■ 

22.27 

20.13 

17.37 

12.23 

17.21 

22.77 

15.55 

15.23 

22.63 

20.26 

15.63 

15.23 

22.84 

20.50 

m mm 

17.19 

24.18 

ISfil 

18.70 

24.63 

A 

14.18 

18.57 

24.80 

17.20 

12.32 

17.23 

24.26 

14.78 

14.10 

18.52 

24.37 

14.71 

14.15 

18.46 

24.75 

wmm 

17.48 

23.38 

19.49 

20.83 

C 

19.15 

21.02 

J 

17.06 

17.19 

23.53 

17.81 

15.61 

19.02 

21.27 

18. 16 

15.15 

19.16 

21.20 

bh 

m 

17.27 

23.13 

I mm 

23.00 

21.28 

ESS 

14.90 

22.86 

21.36 

( 6 

KEl. 

12.26 

17.42 

23.49 

18.77 

14.77 

22.81 

21.07 

18.96 

14.56 

22.67 

21.10 
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Ferric 

Ammonium 

PBAA 

. 

Batch 

Oxide 

Perchlorate 

Polymer 

Aluminum 

17.13 

12.25 

17.22 

22.81 

15.95 

14.57 

19.82 

21.10 

7 

16.08 

14.54 

19.88 

21.19 

/ 

17.21 

12.28 

17.41 

22.82 

15.98 

14.73 

20.11 

20.96 

16.11 

14.55 

20.15 

20.97 

16.95 

12.13 

17.52 

22.58 

18.09 

13.75 

21.71 

22.12 

8 

17.67 

13.96 

21.09 

21.97 

17.04 

12.13 

17.67 

22.61 

17.80 

13.80 

21.46 

21.92 

17.55 

14.06 

21.51 

21.83 

17.18 

12.29 

17.32 

23.20 

16.73 

14.62 

20.75 

21.89 

16.89 

14.58 

20.97 

22.06 

7 

17.15 

12.29 

17.28 

23.34 

16.92 

14.59 

20.67 

21.75 

16.87 

14.52 

20.53 

22.18 

17.50 

12.48 

17.59 

23.05 

16.33 

15.45 

23.82 

20.94 

10 

16.19 

15.25 

23.72 

20.97 

17.72 

12.46 

17.93 

23.18 

16.49 

15.40 

24.04 

21.14 

16.55 

15.  30 

23.98 

21.25 

17.30 

12.27 

17.56 

23.08 

19.19 

14.53 

20.  31 

21.22 

11 

19.06 

14.80 

20.42 

21.45 

17.32 

12.28 

17.79 

22.97 

19.31 

14.75 

20.75 

21.20 

19.33 

14.93 

20.67 

21.43 

17.15 

12.30 

17.43 

23.93 

18.72 

13.63 

21.81 

24.47 

12 

16-58 

13.79 

21.44 

24.20 

17.26 

12.34 

17.66 

24.17 

18.  18 

13.95 

21.52 

24.13 

18.  39 

13.93 

21.60 

24.57 

*Table  16  contains  th**  fixed  counts  used. 
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Batch 


Ferric 

Ammonium 

PBAA 

Oxide 

Perchlorate 

Polymer 

1.1327 

0.8883 

0.8269 

1.1290 

0.8821 

0.8196 

1.1154 

0.9151 

0.8199 

1.1190 

0.9065 

0.8213 

0.9262 

0.8903 

0.9227 

0.9476 

0.8846 

0.9422 

0.9177 

0.8824 

0.9283 

0.9226 

0.8916 

0.9301 

1 .1385 

0.8019 

0.7810 

1.1187 

0.8040 

0.7723 

1.1170 

0.8030 

0.7605 

1.1113 

0.8030 

0.7534 

1.1590 

0.8685 

0.9191 

1.1621 

0.8689 

0.9259 

1.1637 

0.8/40 

0.9302 

1.1693 

0.8708 

0.9335 

0.9290 

0.8145 

0.8968 

0.9397 

0.8089 

0.9129 

0.9579 

0.7892 

0.9038 

0.9394 

0.8130 

0.8969 

0.9032 

0.8280 

0.7508 

0.9123 

0.8255 

0.7555 

0.9046 

0.8299 

0.7637 

0.8956 

0.8420 

0.7684 

1.0740 

0.8409 

0.8690 

1.0653 

0.8427 

0.8663 

1.0770 

0.8337 

0.8657 

1.0683 

0.8441 

0.8640 

0.9370 

0.8820 

0.8068 

0.9592 

0. 8686 

0.8306 

0.9573 

0.8792 

0.8232 

0.9709 

0.8626 

0.8216 

1.0269 

0.8407 

0.8348 

1.0172 

0.8431 

0.8260 

1.0136 

0.8423 

0.8358 

1.0166 

0.8464 

0.8416 

Aluminum 


0.9842 

0.9864 

1.0179 

0.9738 


0.9934 

0.9988 

0.9905 

0.9947 


1.1257 

1.1279 

1.1239 

1.1108 


0.9819 

0.9751 

0.9956 

0.9801 


TABLE  18.  (Concluded) 


Ferric 

Ammonium 

PBAA 

Batch 

Oxide 

Perchlorate 

Polymer 

Aluminum 

1.0716 

0.8080 

0.7383 

1.1006 

10 

1.0809 

0.8181 

0.7415 

1.0989 

1.0746 

0.8089 

0.7458 

1.0965 

1.0707 

0.8144 

0.7477 

1.0909 

0.9015 

0.8443 

0.8645 

1.0878 

11 

0.9077 

0.8289 

0.8593 

1.0758 

0.8969 

0.8323 

0.85:2 

1.0837 

0.8960 

0.8223 

O.P0O8 

1.0717 

0.9161 

0.9024 

0.7992 

0.9779 

12 

0.9230 

0.8921 

0.8128 

0.9889 

0.9494 

0.8847 

0.8206 

1.0016 

0.9386 

0.8861 

0.8177 

C. 9836  1 

TABLE  19.  MEAN  X- 
PBAA  PROPELLANT 


8atch 


1 


RAY  INTENSITY  RATIOS  FOR  Ui  CORED 
ANALYSES  (CONSTANT  PAKiICL'  SIZE) 


Ferric 

Oxide 

R1 

Ammonium 

Perchlorate 

R2 

P8AA 

Polymer 

R3 

Aluminum 

R4 

1.1240 

0.8980 

0.8219 

0.9906 

0.9285 

0.8872 

0.9308 

0.9944 

1.1214 

0.8030 

0.7668 

1.1221 

1.1635 

0.8706 

0.9272 

0.9832 

0.9415 

0.8064 

0.9026 

1.1127 

0.9039 

0.8314 

0.7596 

1.0994 

1.0712 

0.8404 

0.8662 

1.0836 

0.9561 

0.8731 

0.8206 

1.0290 

1.0186 

0.8431 

0.8346 

1.0591 

1.0744 

0.8124 

0.7432 

1.0967 

0.9005 

0.8320 

0.8606 

1.0798 

0 9318 

0.8913 

0.8126 

0.9880 
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Figure  6.  Simple  callbretlon  curve  for  ferric  oxide 
concentration  In  uncured  PBAA  propellants. 
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are  often  referred  to  as  direct  comparison  analyses,  because  the  unknown 
is  compared  directiy  with  the  known  composition  (standard)  that  is  being 
controlied  in  practice.  The  error  of  the  determination  becomes  greater 
as  the  unknown  ferric  oxide  percentage  varies  from  the  0.5%  weight 
nominal  value.  Ihe  same  reasoning  applies  to  the  determination  of  the 
other  ingredients  except  that  the  added  effects  of  matrix  interactions 
must  also  be  considered.  The  ammonium  perchlorate  data  from  Table  20 
are  plotted  in  Figure  8.  The  large  scatter  of  data  about  the  regression 
line  was  indicated  by  the  large  value  of  S and  the  small  value  of 

? e 
R“  in  Table  20.  The  inverted  equation  coefficients  for  estimating 

ingredient  percentages,  as  well  as  the  relative  and  root-mean-square 

errors  (RMSE)  for  estimating  the  ingredient  percentages  are  shown  in 

Table  2i. 

TABLE  2i . SIMPLE  EQUATIONS  FOR  ESTIMATING  INGREDIENT 
CONCENTRATIONS  IN  UNCURED  PBAA  PROPELLANTS 
(CONSTANT  PARTICLE  SIZE) 

*i  ■ dio  + duRi 


Ingredient 

Estimated 

Measured 

Intensity 

Ratio 

Coefficient 

Level 

Coefficient 

Relative 

Error 

(%) 

RMSE 

Ferric  Oxide, 

dio 

-0.05862 

X1 

R1 

dll 

0.54674 

i . i 57 

0.00691 

Ammonium  Per- 

C 

CN 

T3 

27.31403 

chlorate , 

R2 

d22 

48.23926 

1.234 

i. 07700 

PBAA  Poiymer, 

d30 

0.21275 

X., 

3 

R3 

d33 

15.71338 

2.122 

0.34785 

Aluminum,  X, 
4 

d40 

-i. 52741 

R4 

d44 

16.66388 

0.922 

0.19519 

Far  more  accurate  estimates  of  Ingredient  percentages  can  be  made 
using  the  linear  multiple  regression  model  of  Equation  (7)  which  has 
terms  to  compensate  for  matrix  effects.  The  data  for  the  multiple 
regression  response  function  resulting  from  the  use  of  Equation  (7)  are 
given  in  Table  22.  The  standard  errors,  S , and  the  correlation  indices, 

2 e 

R^,  show  an  improvement  over  the  simple  response  function  data  in 

Table  20  for  all  ingredients.  The  Improvement  was  greatest,  however, 
for  the  PBAA  polymer  and  ammonium  perchlorate  determinations. 
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CHLORINE  Ka  INTENSITY  RATIO.  R2 


4 


•t  «7  M M 70 


AMMONIUM  PERCHLORATE  WEIGHT  V Xj 

Figure  8.  Simple  cellbretlon  curve  for  aanonlua  perchlorate 
concentration  In  uncured  PBAA  propellants. 
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TABLE  22.  MULTIPLE  REGRESSION  DATA  FOR  UNCURED  PBAA 
PROPELLANT  ANALYSES  (CONSTANT  PARTICLE  SIZE) 


Coo!  1 i ■'  lent 

— 

) 

! in-,ri’<l  1 rnt 

l.fVfl 

Cocl ! i i i out 

S(> 

l 

C 

K‘ 

i-Vr  r i * • 

hin 

0. 15411 

0.00788 

0.99489 

nxidr 

■’n 

1 . Hr»7  i() 

0.04 4 20 

41  .0183  1 

'*  1 2 

-0.00074 

0.01493 

-0.04848 

hU 

0.00919 

0.01478 

0.48  112 

V. 

-0.008 J2 

0.0.1848 

-0. 40448 

Ammon i um 

Si 

-]  .4  3740 

0.00778 

0.98875 

Ctn  h1nr.itf 

<1.018  12 

0.04471 

0.40081 

^ 1 2 

0.01020 

0.01809 

! . 87894 

h-5 

0.07481 

0.01  >92 

1 . ftOKfin 

-0.00700 

0.0188  1 

-v.  * 

f'HAA  I’nlynuT 

-] .31070 

0.01 1 JO 

0.97944 

h n 

-0.07420 

0.08848 

-1.114  1'. 

hw 

0.02008 

0. 03  1 43 

0.847  18 

hn 

0.08024 

0.07  118 

1.48)80 

h 14 

-0.  00  128 

0.0242  1 

-».  1 1448 

A ! mtti  i num 

hA() 

0.  80788 

0.01  38  4 

0.98  198 

"A! 

-0.  1 13  47 

0.07 444 

- ! .77840 

I),  , 

*4  — 

-0.004  *3 

0.028.’  1 

-0.  1 >.840 

-0.00841 

D.DJW'i 

-0.74701 

hv. 

O.O'ihOS 

0.0371  1 

2.08  7 >0 
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The  regression  equations,  in  algebraic  form,  for  estimating  the 
R_.'s  witli  all  the  X^'s  present  are  as  follows: 

R.  = 0.154  + 1.857  X.  - 0.00074  X.  + 0.00919  X.,  - 0.00832  X, 

1 1 2 3 4 

R.  = -1.437  + 0.0183  X,  + 0.0302  X.  + 0.0256  X-  - 0.0079  X. 

<-  12  3 4 

(14) 

R.  = 1.517  - 0.0743  X.  + 0.0201  X„  + 0.0802  X,  - 0.00328  X, 

3 12  3 4 

R.  = 0.608  - 0.1326  X,  - 0.00442  X„  - 0.00641  X0  + 0.05605  X, 

4 12  3 4 


The  corrected  sums  of  squares  and  cross  products  obtained  during  the 
regression  analysis  are  listed  in  Tables  23(A)  and  23(B). 


TABLE  23(A).  CORRECTED  SUMS  OF  SQUARES  AND  CROSS  PRODUCTS 
(akl)  FOR  UNCURED  PBAA  PROPELLANT  ANALYSES  (CONSTANT 

PARTICLE  SIZE) 


a » 0.02947 

= -0.01256 

aJ3  = -0.04654 

a,.  - 0.05013 
14 

a i - -0.01256 

a2?  = 15.55322 

a23  ■ -8.00854 

a0/  - -7.35425 
24 

a31  * -°*0'654 

a32  = -8.00854 

a33  - 9.27685 

a_,  - -1.01920 
34 

a,,  - 0.05013 
41 

a..  - -7.35425 
42 

a, - - -1.01920 

a. . - 8.17828 
44 

TABLE  23(B).  INVERSE  OF  a^'s  (c  ) FOR  UNCURED  PBAA 
PROPELLANT  ANALYSES  (CONSTANT  PARTICLE  SIZE) 


cn  - 34.72988 

c12  - 1.00642 

c13  - 1.13464 

cn/  - 0.83355 
! 14 

c2I  - 1.00642 

c,2  - 4.29894 

c23  - 4.19776 

c„,  - 4.38275 
24 

c31  - 1.13464 

c 32  - 4.19776 

c 3 3 - 4.20892 

c _ . - 4.29237 
34 

c. .,  - 0.83355 
4i 

c/0  - 4.38275 
42 

c,  - - 4.29237 
43 

c44  " 4,59324 

The  inverse  working  equations  for  estimating  ingredient  percentages 
from  measured  X-ray  intensity  ratios  are  shown  in  Table  24.  In  algebraic 
form  they  are  as  follows: 


X 

X 

X 


1 

2 


X 


3 

4 


-0.1438  + 0.5406  R,  + 0.07935  R„  - 0.0803  Rn  + 0.0867  R, 
1 2 3 4 


38.26  - 0.5767  R,  + 42.57  R„  - 13.11  Rn  + 5.148  R. 

1 2 3 4 

(15) 

8.902  + 0.6984  R,  - 10.48  R„  + 15.69  R~  - 0.4547  R, 
12  3 4 


-7.152  + 1.313  R.  - 2.345  R„  + 0.5705  R_  + 18.40  R, 
12  3 4 


Using  Equations  (15),  the  percentages  of  all  four  ingredients  can  be 
determined  with  a relative  error  of  less  than  1%,  This  type  of  accuracy 
would  provide  for  excellent  control  of  propellant  compositions  during 
a manufacturing  operation.  The  estimated  ingredient  percentages  for 
each  propellant  batch,  using  Equations  (15),  and  the  absolute  (residual) 
errors  of  the  estimates  are  shown  in  Table  25. 


The  estimated  coefficients  in  Equations  (14)  can  be  tested  by 
standard  statistical  methods  to  ascertain  whether  they  are  significantly 
different  from  zero.  Those  that  are  not  can  be  omitted  from  the 
equations  during  the  least  squares  analysis  presumably  to  improve  the 
precision  of  estimating  the  R^j's.  The  stepwise  multiple  regression 

computer  program  readily  performs  this  operation.  The  resulting 
equations,  giving  the  smallest  calculated  errors  of  estimation  and 
containing  only  statistically  significant  coefficients,  are  called  the 
best  set  for  estimating  the  R^'s.  This  best  set  of  regression 

equations  cannot  always  be  inverted,  however,  and  even  when  they  can, 
there  is  no  assurance  that  they  will  produce  a better  set  of  working 
equations  for  estimating  the  X^'s  than  the  set  of  Equations  (15).  The 


best  sets  of  multiple  regression  data  for  this  experiment  are  shown  In 
Table  26.  Comparison  of  these  data  with  those  in  Table  22  for  the  full 
model  shows  that  the  best  sets  provide  only  a slight  improvement. 


All  of  the  coefficients  except  that  for  the  analyte  in  each  of  the 
Equations  (14)  can  be  Interpreted  as  correction  factors  for  the  effects 
of  sample  matrix  ingredients  on  the  X-ray  intensity  ratio.  This  can  be 
illustrated  by  considering  the  equation  for  (ammonium  perchlorate) 

in  Table  26.  If  the  I^'s  are  adjusted  for  the  departure  of  the  individ- 
ual X^'s  from  the  mean  X^  over  all  calibration  standarda,  the  equation 
becomes 
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TABLE  24.  EQUATIONS  FOR  ESTIMATING  INGREDIENT  CONCENTRATIONS 
IN  UNCURED  PBAA  PROPELLANTS  (CONSTANT  PARTICLE  SIZE) 


X.  ,1 . ♦ il . , N,  4 il . ,K  , l <1 . ,K  , ♦ il . , K, 

i id  ill  i , J ill  iii 
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TABLE  26.  BEST  SETS  OF  MULTIPLE  REGRESSION  DATA  FOR  UNCURED  PBAA 
PROPELLANT  ANALYSES  (CONSTANT  PARTICLE  SIZE) 
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ADJ  R2  = R2  - 0.03280  (*3  - X3)  - -2.157  + 0.03762  X2  + 0.03280  X3 

(16) 

The  regression  line  represented  by  Equation  (16)  and  the  ammonium  per- 
chlorate data  points  for  the  calibration  batches  are  given  In  Figure  9. 
The  marked  reduction  in  data  scatter  about  the  regression  iine  for 
ammonium  perchlorate  by  applying  a matrix  correction  for  the  EBAA  polymer 
concentration,  X3*  Is  evident  by  comparing  Figures  8 and  9. 

Another  linear  model  that  might  be  considered  is  one  that  includes 
two-factor  (first  order)  Interactions  such  as  that  shown  in  Table  27. 

This  model  was  fitted  to  the  data  to  determine  whether  the  resulting 
response  function  would  give  a more  precise  estimate  of  the  than 

Equations  (14).  These  regression  equations  cannot  be  inverted  in  a 
straightforward  manner  to  allow  estimation  of  the  ingredient  percentages. 
A significant  two-factor  interaction  means  that  nonparallelism  exists 
between  the  regression  lines  or  planes  involved.  Over  the  relatively 
narrow  ingredient  concentration  ranges  that  would  be  used  for  the 
analysis  of  production  propellants,  the  two-factor  interactions  can 
genera] ly  be  disregarded. 

3.  Analysis  of  Variance 


An  analysis  of  variance  (18)  for  the  uncured  PBAA 
propellant  results  (Table  28)  was  made  to  evaluate  the  various  sources 
of  error,  and  to  permit  later  resolution  of  error  variance  components. 
The  total  sum  of  squares  for  variation  among  the  propellant  batches  has 
been  partitioned  into  the  residual  sum  of  squares,  (SSp)^,  and  the  sum 

of  squares  for  regression,  (SSKeg)^.  The  sum  of  squares  for  regression 
were  determined  as  follows: 


(SSReg) 


.'I 


bik  aikR 


where  a^^  Is  the  corrected  sums  of  cross  products  of  the  R^  with  X^j 

from  the  raw  data.  The  a^p  were  obtained  from  the  corrected  sum  of 

cross  products  (or  the  sum  of  squares  when  1 ■ k)  listed  in  Table  23(A). 
The  residual  sums  of  squares  were  calculated  from 


V,  ■ I 

J"1 


<RU  - v 


so 


TABLE  27.  MULTIPLE  REGRESSION  DATA  WITH  CROSS  PRODUCTS  FOR  UNCURED 
PBAA  PROPELLANT  ANALYSES  (CONSTANT  PARTICLE  SIZE) 


K,  = b.,.+  H..X.  + ...  + h.,X.  + O.^.X,  + ...  + l> . , nx.,x'. 
1 10  ill  1 4 4 I > 1 2 1 1 0 } 'I 


I i i nt 


Cm1!  t ii1  ifnt 
l.iVfl 


C.lH'f  I i ii'iu 


-0.26034 


0.  Of.  I O'. 


■I  . 18490 


t).  I 70.’’. 


•).0|80| 


-0.00901 


1.07180 


1.13040  0.02807  4.48413 


0.01930  1.17041 


-0.07942  0.0)908  -2.0)224 


-0,00038  0.0001  i -4.401 3 1 


0.01)30  2.94983 


-0.03110  0.0  1742  -1.37232 


-0  . 0 0 9 9 4 0.0)1  50  - 2 - 2 14  30 


0.003)7  0.00248  2.103  )2 


-12.00  100 


0.0)103  3.13048 


0.088)7  0.01124  7.87489 


0.12401  4.77212 


-0.00410  0.00244  -1.080)2 


O.OOI87  I -4.81818 


1.14190  0.0  1011  2.12902 


0.0721)  0.02422  ! 2. 4781 1 


-0.  I 1031  0.04712  I -2.  14328 


0.00007  I 10.11 100 


0.00000  0.99744 


0.00004  0.4793 


0.00498  J 0.47730 


TABLE  28.  ANALYSIS  OF  VARIANCE  FOR  UNCURED  PBAA 
PROPELLANT  ANALYSES  (CONSTANT  PARTICLE  SIZE) 


Tlit?  least  squares  analysis  minimized  this  expression.  The  standard 
error,  (S^.).,  was  obtained  from  the  residual  sum  of  squares: 


(19) 


where  n is  the  number  of  calibration  standards  analyzed  (12  in  this 
case),  p is  the  number  of  degrees  of  freedom  for  regression,  and 
n - p - l is  the  number  of  degrees  of  freedom  for  residual  error.  One 
degree  of  freedom  is  required  for  each  term  in  the  regression  model. 
Consequently,  it  is  desirable  to  analyze  a large  number  of  calibration 
batches  so  that  there  will  be  adequate  degrees  of  freedom  remaining 
for  estimating  the  residual  mean  squares  and  standard  errors. 

The  best  estimate  nf  experimental  error  here  is  considered  to  be 
the  interaction  between  calibration  batches,  B,  and  sample  pairs,  T. 

The  sampling  error  was  calculated  from  the  24  duplicate  sample  analyses 
and  represents  sample  repeatability  with  the  instrumental  operating 
conditions  essentially  fixed. 

The  validity  of  the  linearity  assumption  and  the  selected  model 
can  be  checked  by  testing  the  significance  of  the  residual  mean  squares. 
The  residual  errors  for  both  PBAA  polymer  and  aluminum  analyses  are 
significant  at  the  5%  (u  = 0.05)  significance  level.  This  significance 
might  be  partly  due  to  the  fact  that  the  BXT  interaction  is  an  under- 
estimate of  the  true  random  experimental  error.  The  variation  among 
batches,  B,  is  significant  for  all  ingredients.  This  is  an  expected 
result  because  the  compositions  were  purposely  varied  by  the 
experimental  design. 

4.  Cured  Propellant  Results 


X-ray  fluorescence  spectrometry  Is  of  primary 
interest  for  the  rapid,  accurate  analyses  of  uncurcd  propellants  during 
their  manufacture.  At  that  stage,  substandard  compositions  can  either 
be  corrected  or  discarded  without  serious  consequences.  Nevertheless, 
there  are  often  times  when  problems  arise  during  production  or  motor 
storage  and  surveillance  that  an  analysis  of  the  cured  propellant  is 
either  necessary  or  very  desirable.  For  this  reason,  techniques  were 
also  developed  during  this  program  for  the  analysis  of  cured  PBAA 
propellants.  Similar  techniques  can  also  be  applied  to  other  types 
of  cured  propellants. 

Tiie  raw  intensity  measurements  and  the  individual  and  average 
X-ray  Intensity  ratios  for  analyses  of  the  12  calibration  batches 
(Table  15)  are  recorded  in  Tables  29,  30,  and  31.  It  Is  Interesting 
to  compare  the  intensities  of  analytical  emission  lines  from  uncured 
and  cured  propellants  having  the  same  compositions  as  shown  in  Tables 
17  and  29.  The  concentrations  of  ferric  oxide  and  PBAA  polymer  are 
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TABLE  29.  X-RAY  INTENSITY  MEASUREMENTS  IN  SECONDS  FOR  FIXED 
COUNTS*  FOR  CURED  PI1AA  PROPELI.ANT  ANALYSES  (CONSTANT 
PARTICLE  S1XK) 


Ferric 

Ammonium 

PBAA 

Batch 

Oxide 

Perchlorate 

Polymer 

Aluminum 

22.25 

14.13 

19.74 

16.66 

19.71 

14.39 

26.92 

! 5.  16 

I 

20.01 

14.31 

26.80 

15.36 

22.  n 

14.05 

19.98 

16.57 

19.98 

14.28 

27.00 

15.40 

20.04 

14.26 

27.15 

15.61 

22.67 

14.43 

21.01 

18.76 

24.78 

14.75 

26.13 

18.02 

25.06 

14.80 

25.92 

17.82 

22.65 

14.38 

20.86 

18.76 

24.97 

14.73 

26.19 

17.90 

24.68 

14.72 

26.32 

17.84 

22.78 

14.44 

21.08 

18.92 

20.71 

14.91 

32.32 

15.83 

20.61 

15.12 

31.58 

15.57 

5 

22.48 

14.39 

21.21 

18.62 

20.76 

15.02 

31.88 

15.46 

20.15 

15.13 

31.71 

15.21 

22.78 

14.45 

21.89 

18.81 

19.86 

15.21 

27.15 

16.67 

19.72 

15.57 

25.52 

16.59 

22.86 

14.42 

21.56 

18.77 

19.48 

15.37 

25.42 

16.66 

19.76 

15.  19 

25.30 

16.77 

22.71 

14.53 

20.75 

MEM 

23.62 

15.71 

25.23 

IBS 

c 

2 3.75 

15.60 

25.50 

) 

22.47 

14.49 

20.64 

17.19 

24.03 

15.34 

25.76 

14.43 

2 3.98 

15.47 

25.72 

14.24 

22.56 

14.48 

20.90 

17.04 

21.13 

15.  14 

26.78 

14.67 

20.81 

15.20 

26.85 

14.66 

/ 

22.6  3 

14.43 

20.61 

17.08 

20.94 

15.27 

26.56 

14.60 

20.88 

15.17 

26.92 

14.88 
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TABLE  30.  X-RAY  INTENSITY  RATIOS  FOR  CURED  PBAA 
PROPELLANT  ANALYSES  (CONSTANT  PARTICLE  SIZE) 


Batch 


Ferric 

Oxide 


277 

1.1119 

1.1076 

1.1.043 


0.9148 

0.9046 

0.9071 

0.9176 


1.1000 

1.1053 

1.0828 

1.1158 


1.1468 
1. 1552 
1.1735 
1.1569 


0.9615 

0.9562 

0.9351 

0.9370 


1.0681 

1.0846 

1.0807 

1.0838 


1.0156 

1.0215 

0.9996 

0.9932 


1.0733 

1.0579 

1.0602 

1.0713 


Ammonium 

PBAA 

Perchlorate 

Polyme  r 

Aluminum 

0.9819 

0.7333 

1.0989 

0.9874 

0.7366 

1.0846 

0.9839 

0.7400 

1.0760 

0.9853 

0.7359 

1.0615 

0.9783 

0.8040 

1.0411 

0.9750 

0.8106 

1.0527 

0.9762 

0.7965 

1.0480 

0.9766 

0.7926 

1.0514 

0.9685 

0.6522 

1.1952 

0.9550 

0.6675 

1.2152 

0.9580 

0.6653 

1.2044 

0.9508 

0.6688 

1.2241 

0.9503 

0.8062 

1.1284 

0.9281 

0.8578 

1.1338 

0.9382 

0.8482 

1.1266 

0.9493 

0.8522 

1 .1193 

0.9249 

0.8224 

1.2352 

0.9314 

0.8137 

1.2210 

0.9446 

0.8012 

1.1913 

0.9366 

0.8025 

1.2072 

0.9564 

0.7804 

1.1616 

0.9526 

0.7784 

1.1623 

0.9450 

0.7760 

1.1699 

0.9512 

0.7656 

1.1478 

0.9633 

0.7381 

1.1213 

0.9614 

0.7376 

1.1199 

0.9645 

0.7466 

1 . 1037 

0.9739 

0.7413 

1.1113 

0.9760 

0.6499 

1.1519 

0.9756 

0.6580 

1.1483 

0.9710 

0.6789 

1.1544 

0.9671 

0.6751 

1.1572 
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TABLE  30.  (Concluded) 


Ferric 

Ammonium 

PBAA 

Batch 

Oxide 

Perchlorate 

Polymer 

Aluminum 

0.9127 

0.9476 

0.7607 

1.1746 

1] 

0.9063 

0.9494 

0.7603 

1.1718 

0.8871 

0.9613 

0.7303 

1.1528 

0.8986 

0.9524 

0.7465 

1.1685 

0.9291 

0.9897 

0.6894 

1.0526 

12 

0.9195 

1 . 0000 

0.6927 

1.0597 

0.9119 

0.9848 

0.7006 

1.0558 

0.9133 

0.9883 

0.6822 

1.0517 

TABLE  31.  MEAN  X-RAY  INTENSITY 
PROPELLANT  ANALYSES  (CONSTANT 


RATIOS  FOR  CURED  PBAA 
PARTICLE  SIZE) 


Ferric 

Oxide 

R1 

Ammonium 

Perchlorate 

R2 

PBAA 

Polymer 

R3 

Aluminum 

R4 

1.1129 

0.9846 

0.7364 

1.0802 

0.9110 

0.9765 

0.8009 

1.0483 

1.1010 

0.9581 

0.6634 

1.2097 

1.1581 

0.9415 

0.8411 

1.1270 

0.9474 

0.9344 

0.8100 

1.2137 

1.0793 

0. 951 3 

0. 7751 

1.1604 

1.0075 

0.9658 

0.  7409 

1.1140 

1.0657 

0.9724 

0.6655 

1.1530 

0.9012 

0.9527 

0.7494 

1.1669 

0.9184 

0.9907 

0.6912 

1.0550 
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higher  and  the  concentrations  of  ammonium  perchlorate  and  aluminum  are 
lower  in  the  uncured  propellant  surface.  Thus,  the  true  ingredient 
concentrations  in  the  propellant  surface  are  known  only  for  the  cured 
samples.  This  phenomenon  with  the  uncured  propellant  is  due  in  large 
part  to  the  inter facial  tension  between  the  liquid  polymer  and  the 
Mylar  film  on  the  sample  holder.  A practical  consequence  is  that 
theoretical  calibration  procedures  which  depend  on  knowing  the  true 
ingredient  percentage  in  the  analyzed  sample  surface  are  not  applicable 
to  uncured  propellant  analysis. 

Simple  regression  equations  for  cured  propellant  analysis  are 
given  in  Table  32  along  with  several  statistical  parameters  already 
defined.  As  before,  the  best  fit  of  the  model  to  the  data  was  found 
for  ferric  oxide.  If  the  highest  degree  of  accuracy  is  not  required, 
a simple  calibration  of  this  type  is  adequate  for  ferric  oxide  deter- 
minations. The  calibration  curve  for  ferric  oxide  in  Figure  10, 
obtained  from  another  experiment,  shows  that  the  simple  calibration 
procedure  for  ferric  oxide  can  be  used  over  relatively  wide  concentra- 
tion ranges.  The  simpie  inverse  equations  for  estimating  ingredient 
concentrations  are  given  in  Table  33.  Ammonium  perchlorate  can  be 
more  accurately  analyzed  in  cured  propellant  by  this  simple  calibra- 
tion procedure  than  in  uncured  propellant. 

Multiple  regression  equations  for  cured  propellant  analysis  are 
shown  in  Table  34.  All  ingredients  can  be  analyzed  with  a high  degree 
of -precision  using  the  model  of  Equation  (7).  The  corrected  sums  of 
squares  and  cross  products  and  their  inverse,  derived  during  the 
regression  analysis,  are  given  in  Tables  35(A)  and  35(B).  The  analysis 
of  variance,  similar  to  that  for  the  uncured  propellant  analysis,  is 
shown  in  Table  36.  Only  10  of  the  12  calibration  batches  were 
analyzed  in  the  cured  state  hecause  of  proolems  encountered  with  the 
curing  of  two  batches.  This  accounts  for  the  nine  degrees-of-f reedom 
among  batches. 


The  Inverse  working  equations  for  estimating  ingredient  percentages 
are  shown  in  Table  37  along  with  the  estimated  relative  and  root-mean- 
square  errors.  The  relative  errors  for  all  ingredients  are  less  than 
1%.  The  tabulated  Ingredient  percentage  errors  are  in  Table  38.  The 
best  sets  of  multiple  regression  equations  for  cured  propellant  analyses 
are  given  in  Table  39.  These  sets  of  equations  result  in  only  a slight 
improvement  in  the  estimation  precision  for  the  R 's  when  compared 
with  the  regression  equations  in  Table  34,  J 


B,  High-Hate  HTPB  Propellant 


1.  Experimental 

Most  composite  propellants  contain  similar  types, 
but  not  necessarily  the  same,  components.  Therefore,  in  general,  X-ray 
spectometr ic  techniques  developed  for  one  propellant  type  can  be 
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ABLE  33.  SIMPLE  EQUATIONS  FOR  ESTIMATING  INGREDIENT  CONCENTRATIONS 
IN  CURED  P BAA  PROPELLANTS  (CONSTANT  PARTICLE  SIZE) 
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TABLE  SU.  MULTIPLE  REGRESSION  DATA  FOR  CURED  PBAA 
PROPELLANT  ANALYSES  (CONSTANT  PARTICLE  SIZE) 
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TABLE  35(A).  CORRECTED  SUMS  OK  SQUARES  AND  CROSS  PRODUCTS 
(a  ) FOR  CURED  PBAA  PROPELLANT  ANALYSES  (CONSTANT 

PARTICLE  SIZE) 


a = 0.02563 

a12  = 0.04647 

a13  = -0.11970 

a, . - 0.08083 

Ik 

a = 0.04647 

n22  = 14.46875 

a23  = “^*96069 

a..  = -7.56897 
24 

a31  = -0.11970 

:>32  = -6.96069 

a33  = 7.85014 

a,,  = -0.32336 
34 

a/f  = 0.08083 

k I 

a/0  - -7.56897 
42 

a^3  * -0.32336 

= 7.55875 

TABLE  35(B).  INVERSE  OF  a.  . 's  (c.  ,)  FOR  CURED 

kl  kl 

PBAA  PROPELUNT  ANALYSES  (CONSTANT  PARTICLE 
SIZE) 


c « 50.31830 

c^  - -8.25086 

c13  = -6.92344 

c,  . =*  -9.09629 
14 

c21  - -8.25806 

c 22  » 9.82710 

c03  «•  3.01271 

c_,  10.31419 

24 

c3]  = -6.92344 

c 32  = 9.01271 

c33  =*  8.40297 

c_.  * 9.45842 
34 

c..  - -9.09629 
41 

c42  - 10.31419 

c,,  * 9.45842 
43 

e..  - 10.96234 
44 

applied  to  another  typo  with  relatively  minor  modifications.  Of  course, 
a separate  calibration  must  be  performed  for  each  specific  propellant 
whose  composition  will  be  monitored  and  controlled  during  production. 
High  rate  propellants  [19,20 J in  contrast  to  low  rate  propellants, 
contain  ultrafine  ammonium  perchlorate  having  a weight  median  diameter 
of  1 ym  or  smaller,  and  a more  effective  ballistic  modifier  to  enhance 
burning  rate.  Most  composite  propellants  under  development  also  use 
a state-of-the-art  HTPB  binder  system.  This  binder  provides  improved 
structural  integrity  over  earlier  binders  such  as  the  PBAA  type. 

From  an  X-ray  fluorescence  analysis  standpoint,  the  HTPB  hinder 
cannot  be  determined  because  it  has  no  detectable  element  and  the 
determination  of  ammonium  perchlorate  and  aluminum  are  complicated  hy 
the  fact  that  the  ultrafine  ammonium  perchlorate  tends  to  form 
agglomerates  during  propellant  mixing  |20],  These  agglomerates,  which 
can  vary  in  size  among  propellant  hatches,  affect  the  precision  of  the 
aluminum  determination  because  the  AIK  emission  line  intensity  is 
sensitive  to  the  agglomerate  size. 
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TABLE  37.  EQUATIONS  FOR  ESTIMATING  INGREDIENT  CONCENTRATIONS 
IN  CURED  PBAA  PROPELLANTS  (CONSTANT  PARTICLE  SIZE) 
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TABLE  39.  BEST  SETS  OF  MULTIPLE  REGRESSION  DATA  FOR  CURED 
PROPELLANT  ANALYSES  (CONSTANT  PARTICLE  SIZE) 


■1 
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A central  composite  experimental  design  [11]  was  used  to  evaluate 
applicability  of  the  X-ray  fluorescence  method  to  bigb-rate  propellant 
analysis.  For  this  report,  however,  the  design  has  been  broken  down 
Into  three  separate  experiments  consisting  of  two  factorial  designs  and 
the  full  central  composite  design.  The  compositions  of  the  calibration 
batches  for  the  first  experiment  are  shown  in  Table  40.  Each  propellant 

3 

batch  also  contains  HTPB  binder.  This  is  a full  2 factorial  design 
with  independent  ingredient  concentrations.  The  liquid  ballistic 
modifier  in  this  case  contained  sulfur;  therefore,  it  could  be  deter- 
mined in  the  propellant.  The  ammonium  perchlorate  oxidizer  was  composed 
of  a bimodal  blend  (55/45  by  weight)  of  nominal  0.5-  and  90-pm  sizes. 

3 

Experiment  -2,  as  shown  in  Table  41,  was  also  a 2 factorial  design. 

The  ingredient  percentages  were  the  same  as  for  Experiment  -1  except 
that  the  ammonium  perchlorate  consisted  of  a 75/25  by  weight  bimodal 
blend  of  0. 5-  to  90-pm  ammonium  perchlorate. 


TABLE  40.  COMPOSITIONS  OF  CALIBRATION  BATCHES  (WEIGHT  1)  FOR 
DETERMINING  INGREDIENT  CONCENTRATIONS  IN  HIGH  RATE  HTPB 


PROPELLANT 

(EXPERIMENT  -1  , 

PARTICLE  SIZE  CONSTANT) 

m 

35-pm 

Ballistic 
Mod  i f i er, 

Ammonium 

Aluminum,  X, 

X0 

Perchlorate,  X * 

1 

2 

J 

i 

12.80 

6,20 

70.65 

2 

12.80 

9.80 

68.85 

3 

13.20 

6.20 

68.85 

4 

12.80 

6.20 

72.43 

5 

15.20 

9.80 

68.85 

6 

12.80 

9.80 

72.45 

7 

15.20 

6.20 

72.45 

8 

15.20 

9.80 

70.65 

*Bimodal  blend  of  45 X - 0.5-pm  and  552  - 90-pm 
ammonium  perchlorate. 

NOTE:  Remainder  of  propellant  composition  is  HTPB  binder. 
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TABU:  41.  COMPOSITIONS  OF  CALIBRATION  BATCHES  (WEIGHT  %)  FOR 
DETERMINING  INGREDIENT  CONCENTRATIONS  IN  HIGH  RATE  HTPB 
PROPELLANT  (EXPERIMENT  -2,  PARTICLE  SIZE  CONSTANT) 


Batch 

35-pm 

Aluminum,  X^ 

Bal 1 1st  I c 
Mod i f ler 

X2 

Ammonium 
Ferchlorate*,  X^ 

1 

12.80 

6.20 

70.65 

2 

12.80 

9.80 

68.85 

3 

15.20 

6.20 

68.85 

4 

12.80 

6.20 

<r 

CM 

r-- 

5 

15.20 

9.80 

68,85 

6 

12.80 

9.80 

72.45 

7 

15.20 

9.80 

70.65 

*Blmodal  blend  of  75%  - 0.5-um  and  25%  - 90-uin 
ammonium  perchlorate. 

NOTE:  Remainder  of  propel iant  composition  is  HPTB  binder. 


The  analytical  emission  lines  and  the  instrumental  operating 
conditions  used  are  given  in  Table  42.  A fixed  time  measurement 
technique  was  used.  Pulse  height  analysis  was  required  only  for  the 
determination  of  aluminum.  Only  uncured  propellant  samples  were 
analyzed.  They  were  prepared  and  analyzed  in  the  same  manner  as 
di  •scribed  for  the  low-rate  PBAA  propellant  in  Section  11. 

2.  Results  and  Discussion 

The  average  X-ray  Intensity  ratios  for  the  three 
analyzed  ingredients  in  each  calibration  batch  are  shown  in  Tables  43 
and  44  for  Experiments  -1  and  -2,  respectively.  Regression  analyses 
were  made  using  tiiese  X-ray  intensity  responses  and  the  composition 
data  in  Tables  40  and  41. 

The  regression  data  for  the  multiple  linear  regression  model  with 
all  components  present  are  shown  in  Tables  45  and  46  for  Experiment  -1 
and  in  Table  47  for  Experiment  -2.  The  precision  for  estimating  the 
Rj's  is  less  than  it  is  for  low-rate  propellant  analysis.  The  correla- 
tion Indices  for  aluminum  and  ammonium  perchlorate  determinations 
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TABLE  42.  INSTRUMENTAL  PARAMETERS  FOR  DETERMINING  HIGH  RATE 
IITPB  PROPELLANT  INGREDIENT  CONCENTRATIONS  (PARTICLE  SIZE 
CONSTANT) 


Ingredien. 

Emission 

Line 

Anal yzing 
Crystal 

Peak 
Angle 
(deg  20) 

Fixed 
T ime 
(sec) 

Pulse 

Height 

Analysis 

Aluminum 

Aluminum  K 

PET 

145.10 

50 

Yes 

Ball  ist ic 
Modi f ier 

Sulfur  K 

a 

PET 

75.92 

50 

No 

Ammonium 

Perchlorate 

Chlorine  K 

|2 

PET 

60.55 

50 

No 

NOTE:  Chromium  target  X-ray  tube  operated  at  40  kV  and 

'30  mA  (constant  potential). 


TABLE  43.  MEAN  X-RAY  INTENSITY  RATIOS*  FOR  UNCURKD 
HIGH  RATE  HTPB  PROPELLANT  ANALYSES  (EXPERIMENT  -I, 
PARTICLE  SIZE  CONSTANT) 


Batch 

Aluminum,  R^ 

Ballistic 

Modifier, 

*2 

Ammonium 
Perchlorate,  R^ 

1 

0.288 

1.243 

0.747 

2 

0.274 

1.936 

0.726 

3 

0.353 

1.255 

0.715 

4 

0.264 

1.255 

0.76' 

5 

0.366 

1.976 

0.691 

6 

0.282 

1.974 

0.727 

7 

0.191 

1.093 

0.816 

8 

0.323 

1.946 

0.716 

*Eath  ratio  is  the  mean  of  three  sample  determinations 
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TABLE  44.  MKAN  X-RAY  INTENSITY  RATIOS*  FOR  UNCURED  HIGH 
RATE  HTPB  PROPELLANT  ANALYSES  (EXPERIMENT  -2,  PARTICLE 
SIZE  CONSTANT) 


Batch 

Aluminum,  R^ 

Ballistic 
Modi f ier, 

R2 

Ammonium 
Perchlorate,  R^ 

1 

0.112 

1.014 

0.843 

2 

0.096 

1.582 

0.817 

3 

0.143 

0.989 

0.837 

4 

0.  1 33 

1.061 

0.797 

5 

0.132 

1.541 

0.802 

6 

0.163 

1.696 

0.736 

7 

0.151 

1.575 

0.780 

*Each  ratio  is  the  mean  of  three  sample  determinations. 

Indicate  an  unacceptable  fit  of  the  model  to  the  data.  The  model  itself 
is  believed  to  be  adequate.  The  relatively  poor  fit  for  these  ingre- 
dients Is  attributed  to  the  adverse  Influence  of  the  ultrafine  ammonium 
perchlorate  agglomeration  in  the  propellant.  This  agglomeration  problem 
was  compounded  In  these  experiments  by  the  propellant  processing 
I procedure  used.  Shear  stress  during  mixing  was  kept  high  by  the  con- 

trolled addition  of  the  solids  In  an  attempt  to  hreak  up  agglomerates. 

This  was  later  found  to  make  the  agglomerate  formation  problem  worse. 

At  the  present  time,  with  improved  procedures,  high-rate  propellants 
can  be  processed  to  give  substantially  better  results  than  were  obtained 
in  these  experiments.  However,  the  overall  precision  for  high -rate 
propellant  analysis  hy  X-ray  spectrometry  is  still  less  than  that  for 
low-rate  propellant  analysis,  particularly  for  the  aluminum  determination. 

A model  that  includes  first -osder  interaction  terms  was  also  con- 
sidered in  Experiment  -1.  The  regression  analysis  data  are  presented 
in  Table  4ft.  This  model  appears  to  be  superior  to  the  one  that  con- 
tains only  main  effect  terms  (as  in  Table  45).  But,  the  goodness-of-f it 
might  be  misleading  hecause  there  is  only  one  degree-of-freedom 
I remaining  for  estimating  the  residual  error.  If  all  of  the  degrees-of- 

freedom  (seven  In  this  case)  had  been  used  for  regression,  then  the 
correlation  index  would  necessarily  have  a value  of  one.  That  Is,  the 
regression  planes  would  have  been  forced  through  all  of  the  points. 

Because  of  the  lack  of  fit  In  most  cases,  the  regression  equations  were 
not  Inverted  to  obtain  working  equations  for  estimating  ingredient 
percentages. 
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TABLE  45.  MULTIPLE  REGRESSION  DATA  FOR  UNCURED  HIGH  RATE  HTPB  PROPELLANT 
ANALYSES  (EXPERIMENT  -1,  PARTICLE  SIZE  CONSTANT) 
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TABLE  46.  MULTIPLE  REGRESSION  DATA  WITH  CROSS  PRODUCTS  FOR  UNCURED 
HIGH-RATE  HTPB  PROPELLANT  ANALYSES  (EXPERIMENT  -1,  PARTICLE  SIZE 
CONSTANT) 


K,  “ h . ,,  + t> . . X , + ...  + |>  X . + h,,X.X-  + ...  4 h..X.,X, 
I 1 0 ill  ill  1412  i 6 2 j 


Ooef f ic lent 

1 !i>:rr<l  lent 

l.cve  1 

Cue  1 t i c it'll  l 

Sb 

l 

i.' 

IT 

A 1 UR)  in  (in 

blo 

-6.  1.4491 

hi 

0.07984 j 

0.0891  1 

0.89180 

bl2 

-0.4  160  1 

0.06442 

-7.07901 

h13 

0.  1 Oh  )(» 

0.02116 

4.97619 

0.00461 

0. 9990 i 

'*  1 4 

0. 00701 

0.0009  1 

7. 18061 

b15 

-0.01202 

0.00120 

- 10.01667 

b14 

0.0011 1 

0.00080 

6.41210 

B.i  1 1 i st  Ic 
Mud  i 1 i«T 

b20 

-11 .07177 

b2 1 

1 .0811)7 

0.  17710 

2.87900 

* * 2 

-0. 2097 j 

0.27214 

-0.76914 

b2i 

0. 1 7911 

0.090  >7 

1 .9H6  19 

0.01968 

0, 99966 

b24 

0.001)94 

0.00  19  4 

I . 76142 

l>21 

-0.011.4  I 

0.00109 

-J.  22790 

*26 

0.00447 

0. 00  1 19 

1.  11818 

\nrnu»n  ( urn 
IVrch  1 or.it  t- 

h 10 

1.80441 

b 11 

-0.424  10 

0.01742 

-7.  18941 

12 

0.  ’HI  10 

0.04110 

6.87470 

’ll 

-0.0481 7 

0.01  176 

- 1. 100  7 1 

0.00  100 

0.99911 

b J4 

-0.00141 

0.00060 

-1.68111 

0.00642 

0.00077 

8.  11760 

b16 

-0.00  149 

0.00011 

-6.84  114 
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V.  DETERMINATION  OF  PARTICLE  SIZES  WITH  INGREDIENT 
CONCENTRATIONS  CONSTANT 

A.  General  Considerations 

The  X-ray  emission  line  intensities  from  elements  in  a 
propellant  sample  are  sensitive  not  only  to  the  percentages  of  the 
corresponding  ingredients  but  also  to  the  average  particle  sizes  of  the 
ingredients.  The  fact  that  in-situ  particle  size  measurements  of 
propellant  solid  ingredients  could  be  made  by  X-ray  spectrometry  was 
recognized  and  applied  on  a semiquantitati ve  basis  by  this  command 
several  years  ago  {1).  Normally,  particle  size  or  grain  effects  are 
considered  to  be  a disturbing  effect  in  X-ray  spec trometric  analysis; 
samples  are  treated  either  to  minimize  or  eliminate  the  effects.  In 
contrast,  for  propellant  analysis,  the  ability  to  measure  solids  particle 
sizes  In  finished  propellants  is  an  important  feature  of  the  X-ray 
fluorescence  method.  The  reason  for  this  is  that  the  solids  particle 
sizes  can  have  a pronounced  influence  on  propellant  ballistic,  mechanical, 
and  rheological  properties.  In  fact,  during  the  tailoring  of  propellant 
properties  the  particle  sizes  of  the  solids,  particularly  ammonium 
perchlorate  and  aluminum,  are  carefully  chosen  along  with  other  components 
to  provide  the  required  propellant  properties.  Moreover,  the  particle 
sizes,  specifically  the  particle  size  distribution,  of  these  solids  must 
be  closely  controlled  during  propellant  manufacture. 

The  importance  and  implications  of  particle  size  effects  in  X-ray 
spectrometry,  both  qualitatively  and  quantitatively,  have  been  well- 
reported  in  the  literature  (21 ,22].  These  and  other  published  methods, 
however,  are  not  di  ectly  applicable  to  composite  propellant  analysis, 
Jenkins  [23]  gives  a discussion  of  particle  size  effects  in  X-ray 
spectrometry  and  their  origin.  The  particle  size  effects  in  propellants 
arise  hecause,  as  the  particle  size  of  a component  of  a bimodal  or  multi- 
modal blend  changes,  its  concentration  in  the  analyzed  propellant  surface 
also  changes  relative  to  the  »thcr  size  components.  As  the  average 
particle  size  of  the  component  decreases,  its  concentration,  and  hence 
its  analytical  emission  line  intensity,  increases.  Conversely,  as  the 
average  particle  size  of  the  component  increases,  the  analytical  emission 
line  intensity  decreases.  It  was  reported  in  earlier  work  [24]  that  the 
C1K  emission  iine  Intensity  f'om  ammonium  perchlorate  in  propellant  is 

linearly  related  to  the  ammonium  perchlorate  weight  mean  diameter.  Also, 

the  C1K  intensity  is  a linear  function  of  the  percentage  of  a fine 
u 

ammonium  perchlorate  fraction  in  a bimodal  blend  of  fine  and  coarse  size 
f motions. 

The  effective  depth  of  penetration  of  the  emission  line  In  the 

analyzed  propellant  surface  also  plays  an  important  role  with  respect 

to  particle  size  effects,  because  it  determines  the  volume  of  propellant 

analyzed.  The  effective  depth  of  penetration  for  long  wavelength 

radiation  such  as  ClK  , SK  , and  AIK  Is  of  the  order  of  50  pm  or  less, 

u u a 
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becoming  smaller  as  the  wavelength  increases.  If  the  effective  depth 
of  penetration  is  less  than  the  diameters  of  the  particles  being 
measured,  then  particle  size  effects  are  especially  pronounced. 

Another  factor  to  he  considered  is  particle  shielding.  For  example, 
AlK^  radiation  may  have  to  pass  through  ammonium  perchlorate  particles 

before  it  is  measured.  Consequently,  as  the  ammonium  perchlorate 

particle  size  decreases,  the  measured  AIK  intensity  from  aluminum 

a 

particles  also  decreases,  even  though  the  aluminum  particle  size  is 
held  constant.  Furthermore,  the  magnitude  of  this  effect  depends  on 
the  mass  absorption  coefficient  of  the  matrix  (and  matrix  particles) 
for  the  measured  radiation.  Because  of  its  high  absorption  by  ammonium 
perchlorate  and  its  small  effective  penetration  depth,  AIK^  radiation 

is  very  sensitive  to  ammonium  perchlorate  particle  size  changes.  It 

can  be  used  to  measure  ammonium  perchlorate  weight  median  sizes  down 

to  at  least  1 pm.  The  chlorine  K radiation,  on  the  other  hand,  is  not 

« 

sensitive  to  ammonium  perchlorate  particle  size  change  below  a weight 
median  diameter  of  approximately  5 urn. 

The  purpose  of  this  experiment  was  to  develop  quantitative  proce- 
dures for  determining  ammonium  perchlorate  and  aluminum  particle  size 
variations  in  uncured  and  cured  propellants.  These  variations  can  occur 
through  a weighing  error  or  uncontrolled  alteration  of  the  particles 
during  propellant  processing.  The  percentage  of  a fine  fraction  in  a 
himodal  blend  of  sizes  was  used  here  as  a measure  of  particle  size 
variations.  The  actual  weight  (volume)  mean  diameters  of  the  particles, 
if  desired,  can  readily  be  obtained  for  a specific  propellant  formula- 
tion by  establishing  the  linear  intensity-particle  diameter  relationship. 
This  is  not  required  for  the  analysis  and  control  of  production 
propellant  compositions. 

B.  Low-Kate  PBAA  Propellant 
1.  Fxperimental 

2 

The  experimental  design  used  was  a 2 factorial 
with  two  additional  design  points  as  shown  in  Tables  48  and  49.  The 
two  additional  points  have  the  same  composition  as  the  midpoint  of 
the  design  shown  in  Table  13.  The  low  level  for  each  ingredient  was 
composed  of  the  larger  average  particle  size  and  consequently  the  high 
level  was  composed  of  the  smaller  average  particle  size.  Only  the 
average  particle  sizes  of  ammonium  perchlorate  and  aluminum  were  varied, 
with  all  ingredient  percentages  held  constant.  The  average  particle 
size  of  annonlum  perchlorate  was  varied  by  changing  the  ratio  of  nominal 
20-  and  200-pm  particle  size  fractions.  The  average  aluminum  particle 
size  was  varied  by  using  two  different  aluminum  powders  having  nominal 
sizes  of  9-  and  32-um  weight  median  diameters.  The  analytical  condi- 
tions and  instrumental  parameters  for  this  experiment  are  listed  in 
Table  50. 
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TABLE  48.  FACTORS  AND  FACTOR  LEVELS  FOR  I’ BAA 
PROP  ELI  ANT  PARTICLE  SIZE  CALIBRATION 
BATCHES  (CONSTANT  CONCENTRATIONS) 


Factor 

Symbol 

Low  Level* 

High  Level* 

Ammonium 

E 

20%  - 20  um 

60%  - 20  pm 

Perchlorate 

80%  - 200  um 

20%  - 200  pm 

Aluminum 

F 

0%  - 9 pm 

100%  - 9 pm 

100%  - 32  pm 

0%  - 32  pm 

*Weight  percent  on  a total  ingredient  hasis. 


2.  Uneured  Propellant  Results 

The  individual  X-ray  intensity  ratios  for  analyses 
of  the  calibration  standards  are  recorded  in  Table  5i;  the  average 
values  used  in  the  regression  analysis  are  recorded  in  Table  52. 

Despite  the  fact  that  only  ammonium  perchlorate  and  aluminum  particle 
sizes  were  varied,  all  of  the  ingredient  emission  ilne  intensities  were 
measured  to  evaluate  the  effects  of  the  particle  size  variations  on  ail 
inured  lent  determinations. 

The  regression  equations  for  estimating  the  R^’s  with  the  ammonium 

perchlorate,  W^,  and  aluminum,  W^,  particle  size  fine  fractions  as  the 

Independent  variahies  are  given  in  Table  53.  A good  fit  of  the  modei 
to  the  data  was  found  for  aii  ingredients.  The  magnitudes  and  signs 
of  the  partial  regression  coefficients  indicate  the  effects  of  the 
particle  size  variations  on  each  estimated  analytical  emission  line,  R^. 

A negative  coefficient  means  that  the  particular  emission  line  intensity 
decreases  as  the  average  particle  size  decreases.  Conversely,  a posi- 
tive coefficient  means  that  the  intensity  increases  as  the  average 
particle  size  decreases  (high  fine-fraction  percentage).  For  example, 
as  the  particle  size  of  ammonium  perchlorate  decreases,  the  CJK  line 

Intensity  increases;  whereas"  the  SK  , A:K  , and  FeK  ilne  intensities 
1 a o n 

decrease.  Likewise,  when  the  aluminum  particle  size  decreases,  the 

AIK  emission  line  intensity  increases,  and  the  emission  line  intensities 
a 

from  the  other  elements  decrease.  This  is  a general  phenomenon.  The 

2 

high  correlation  index,  R^,  shows  that  the  intensity-particle  size 

(weight  1 fine-fraction)  relationship  is  well  represented  by  a linear 
model. 
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TABLE  49.  COMPOSITIONS  OF  CALIBRATION  BATCHES  FOR  PARTICLE 
SIZE  ANALYSIS  (WEIGHT  7.)  (CONSTANT  CONCENTRATIONS) 
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TABLE  51.  X-RAY  INTENSITY  RATIOS  FOR  UNCURED  PBAA 
PROPELLANT  PARTICLE  SIZE  ANALYSES  (CONSTANT 
CONCENTRATIONS) 


Ferric 

Ammonium 

PBAA 

Batch 

Oxide 

Perchlorate 

Polymer 

Aluminum 

1.0368 

0.8180 

0.8456 

1.1829 

1 

1.0357 

0.8199 

0.8556 

1.1818 

1.0379 

0.8319 

0.8403 

1.1533 

1.0449 

0.8147 

0.8675 

1.1660 

1.0384 

0.3182 

0.3502 

1.1608 

9 

1.0362 

0.8203 

0.8429 

1.1766 

1.0392 

0.8251 

0.8470 

1.1694 

1.0419 

0.8219 

0.8466 

1.1517 

1.1801 

0.5904 

0.9676 

1.3864 

1.1576 

0.5778 

0.9308 

1.4256 

J 

1.1636 

0.5829 

0.9370 

1.3726 

1.1514 

0.5895 

0.9203 

1.3833 

0.9260 

1.0265 

0.7657 

0.9754 

A 

0.9318 

1.0023 

0.7723 

0.9712 

0.9386 

1.0339 

0.7667 

0.9775 

0.9218 

1.0500 

0.7475 

0.9614 

1.1848 

0.9516 

1.3721 

0.7298 

1.1520 

0.9494 

1.3181 

0.7402 

1.1692 

0.9673 

1.3526 

0.7622 

1.1749 

0.9572 

1.3424 

0.7563 

0.9503 

1.3925 

1.0069 

0.4709 

A 

0.9615 

1.3702 

0.9981 

0.4562 

0.9382 

1.4177 

0.9928 

0.4863 

0.9400 

1.4047 

0.9853 

0.4846 
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TABLE  52.  MEAN  X-RAY  INTENSITY  RATIOS  FOR  UNCURED  PBAA 

PROPELLANT  PARTICLE  SIZE  -ANALYSES  (CONSTANT  CONCENTRATIONS) 


Batch 

Ferric 

Oxide 

R1 

Ammonium 

Perchlorate 

R2 

PBAA 

Polymer 

R3 

Aluminum 

R4 

1 

1.0388 

0.8211 

0.8522 

1.1715 

2 

1.0389 

0.8214 

0.8467 

1.1646 

3 

1.1632 

0.5851 

0.9389 

1.3920 

4 

0.9296 

1.0282 

0.7630 

0.9714 

5 

1.1702 

0.9564 

1.3463 

0.7471 

6 

0.9475 

1.3963 

0.9958 

0.4745 

The  best  sets  of  multiple  regression  data  with  the  first  order 
interaction  included  are  given  in  Table  54.  These  equations,  except 

for  ammonium  perchlorate,  reduce  the  error  of  estimating  R^'s  when  compared 

with  the  equations  in  Table  53,  but  they  cannot  be  inverted  in  a straight- 
forward manner  for  estimating  WjandW^.  Tables  55  end  56  list  equations 

for  estimating  and  from  the  measured  intensltias,  R^.  Although  the 

root-mean-squere  errors  ara  small,  the  equations,  as  pointed  out  earlier, 
must  be  used  with  caution  because  the  intensities  which  ara  not  controlled 
were  used  as  the  independent  variables  for  the  leest  squares  analysis.  The 
particle  siza  estimates  in  weight  fractions  for  the  individual  calibra- 
tion batches,  using  the  best  sets  of  equations,  are  given  in  Table  57. 

A truer  test  of  the  validity  of  the  estimation  equations  would  be  to 
analyze  propellants  not  included  in  the  calibration.  It  should  be 
mentioned  that  the  equations  in  Table  53  can  be  inverted  to  give 
adequate  working  equations  for  astimatir.g  particle  sizes. 

The  various  intensity-particle  size  relationships  are  shown 
graphically  as  two-way  plots  in  Figures  11  through  16.  In  practice, 
the  in-situ  measurement  of  ammonium  perchlorate  particle  size  is  of 
primary  interest.  As  illustrated  in  Figures  11  and  12,  both  CIK^  and 

AIK  intensities  are  very  sensitive  to  ammonium  perchlorate  particle 

size  changes.  The  AlK^  radiation  is  actually  more  sensitive  to  the 

ammonium  perchlorate  size  change  than  the  CIK^  radiation  in  most 

applications.  Therefore,  the  AlK^  line  is  the  one  of  choice  for 
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TABLE  55.  EQUATIONS  FOR  ESTIMATING  AMMONIUM 
PERCHLORATE  SIZE  FRACTIONS  (W  ) IN  UNCURED 

PBAA  PROPELLANTS  (CONSTANT  CONCENTRATIONS) 


W2  X 10 

= d2Q  + d2lRl 

t-  d22R2  + d2;}R3 

+ d24R4 

Intensity 

Coefficient 

Ratios 

Level 

Coefficient 

RMSE 

d20 

-2.46759 

0.03464 

R2 

d22 

1.82788 

R4 

d24 

1.17431 

d20 

1.73962 

0.00516 

R1 

d21 

-1.31374 

R2 

d22 

0.15815 

R3 

d23 

-0.11565 

d20 

1.46859 

0.00538 

R1 

d21 

-1.41677 

R2 

d22 

0.28957 

R4 

d24 

0.14740 

a. 

NJ 

O 

5.19576 

0.00663 

R2 

d22 

-1  51767 

R3 

d23 

-1.59034 

R4 

d24 

-1.87968 

d20 

2.06579 

0.00497 

R! 

d21 

-1.18976 

R3 

d23 

-0.23482 

R4 

d24 

-0.17739 
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TABLE  56.  EQUATIONS  FOR  ESTIMATING  ALUMINUM  SIZE  FRACTIONS 
(W^)  IN  UNCURED  PBAA  PROPELLANTS  (CONSTANT  CONCENTRATIONS) 


W,  x 10  = d, „ + d,  „R.  + d , „R„  + d,„R„  + d,.R, 


4U  411 

l 4J  J 

44.  4 

Intensity 

Coefficient 

Ratios 

Level 

Coefficient 

RMSE 

d40 

-5.48833 

0.10630 

• R2 

d42 

2.80821 

R4 

d44 

3.57703 

d40 

7. 32703 

0.00550 

R1 

d41 

-4.00174 

R2 

d42 

-2.27786 

R3 

d43 

-0.35226 

d40 

6.50145 

0.00561 

R1 

d41 

-4.31555 

R2 

d,2 

-1.87756 

R4 

d44 

0.44900 

d40 

17.85463 

0.02107 

R2 

d42 

-7.38253 

R3 

d43 

-4.84428 

«4 

d44 

-5.72560 

d40 

2.62927 

0.1020 

R1 

d41 

-5.78744 

K3 

d43 

1.65222 

K4 

d44 

2.55495 
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Figure  11.  Chlorine  K intensity  es  a function  of  ammonium  perchlorate 

a 

particle  size  fractions  for  constant  aluminum  particle  size  (uncured 
PBAA  propellant). 


WEIGHT  % Of  K-pm  AMMONIUM  PERCHLORATE , W2 

Figure  12.  Aluminum  Kq  Intensity  es  e function  of  ammonium  perchlorate 

particle  size  fractions  for  constant  aluminum  particle  size  (uncured 
PBAA  propellant). 
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Figure  13.  Aluminum  intensity  «s  a function  of  aluminum  particle 

size  fractions  for  constant  ammonium  perchlorate  particle  size 
(uncured  PE^A  propellant).. 
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WEIGHT  X Of  9-pm  ALUMINUM.  W4 

Figure  14.  Chlorine  intensity  ss  s function  of  alualnua  particle 

size  frsctlons  for  constant  aaaonlum  perchlorate  particle  size 
(uncured  PBAA  propellant). 
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IRON  INTENSITY 


^igure  IS. 


Iron  Kq  intensity  ns  « function  of  smnonlua  perchlorate 


particle  size  fractions  for  constsnt  aluminum  particle  size 
(uncured  PBAA  propellant). 
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•0 


20 

WEIGHT  % Of  20-fim  AMMONIUM  PERCHLORATE.  W^ 

Figure  16.  Sulfur  Kq  Intensity  es  e function  of  aanonlun  perchlorete 

particle  size  fractions  for  constant  al  uni  nun  psrtlde  size  (uncured 
PBAA  propellant). 
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ammonium  perchlorate  particle  size  measurements  when  the  aluminum 
particle  size  is  held  constant.  Figure  12  shows  that,  because  of  a 
significant  interaction  between  aluminum  and  ammonium  perchlorate 
particle  sizes,  the  AIK^  intensity  measurement  becomes  more  sensitive 

to  ammonium  perchlorate  size  changes  as  the  particle  size  of  the  aluminum 
powder  itself  decreases.  There  also  appears  to  be  significant  inter- 
actions in  Figures  13  and  16.  Notice  that  the  FeK  and  SK  line  inten- 

a a 

sities  in  Figures  15  and  16  are  affected  less  by  the  particle  size  change 
than  the  AlK^  and  CIK^  lines.  The  aluminum  particle  size  change  has 

little  effect,  in  particular,  on  the  FeK^  emission  line  intensity,  as 

shown  in  Figure  15  because  of  the  large  penetration  depth  of  FeK^ 

radiation  relative  to  the  aluminum  particle  size,  and  the  small  absorp- 
tion coefficient  of  aluminum  for  FeK  radiation. 

a 

3.  Cured  Propellant  Results 

The  cured  propellant  results  are  similar  to  those 
for  uncured  propellant,  except  that  the  magnitudes  of  the  relative 
particle  size  effects  are  different  because  the  compositions  of  the 
analyzed  surfaces  of  cured  and  uncured  propellants  are  entirely  different 
as  has  already  been  explained.  The  composition  of  the  cured  propellant 
surface  accurately  represents  the  composition  of  the  bulk  of  the 
propellant. 

Individual  and  average  X-ray  intensity  ratio  measurements  for  the 
cured  propellant  analyses  are  recorded  in  Tables  58  and  59,  respectively. 
The  multiple  regression  equations  for  estimating  R^  with  and  as 

independent  variables  are  shown  in  Table  60.  Except  for  ammonium  per- 
chlorate, the  standard  errors,  Se,  for  estimating  R^  are  smaller  than 

those  for  uncured  propellant  analysis  (Table  53).  The  signs  of  the 
partial  regression  coefficients  are  the  same  in  each  case.  The  magnitude 
of  the  particle  size  effects  on  emission  line  intensities,  however,  is 
significantly  smaller  for  cured  propellant  analysis.  Thus,  the  measure- 
ments are  less  sensitive  to  ammonium  perchlorate  and  aluminum  particle 
size  changes  when  analyzing  cured  propellants. 

The  best  sets  of  multipie  regression  data  including  the  inter- 

action term  are  given  in  Table  61.  The  equations  for  estimating  Wj  and 
are  given  in  Tables  62  and  63,  and  the  resulting  particle  size  esti- 
mates using  these  equations  are  given  in  Table  64.  As  for  uncured 
propellant,  the  least  squares  analysis  was  made  using  the  and 

as  the  dependent  variables  and  the  measured  intensity  ratios  as  the 
independent  variables. 


Batch 


Ferric 

Ammon ium 

PBAA 

Oxide 

Perchlorate 

Polymer 

Aluminum 

1.0055 

0.9586 

0.7308 

1.1955 

1.0062 

0.9624 

0.7310 

1.1654 

1.0138 

0.9586 

0.7238 

1.1793 

1.0055 

0.9623 

0.7162 

1.1520 

1 . 0007 

0.9694 

0.7119 

1.1700 

1 . 0068 

0.9612 

0.7191 

1.1851 

1.0117 

0.9667 

0.7114 

1.1704 

0.9986 

0.9733 

0.7104 

1.1584 

1.0330 

0.9504 

0.7284 

1.2553 

1.03, .3 

0.9556 

0.7158 

1.2277 

1.03/9 

0.9577 

0.7115 

1.2274 

1.0358 

0.9567 

0.7236 

1.2227 

0.9728 

0.9793 

0.6994 

1.1015 

0.9760 

0.9760 

0.7065 

1.1110 

0.9824 

0.9770 

0.7013 

1.0899 

0.9760 

0.9765 

0.6992 

1.0998 

1.0366 

1.0460 

0.9621 

1.0817 

1.0438 

1.0325 

0.9811 

1.0948 

1.0554 

1.0325 

0.9938 

1.0856 

1.0488 

1.0417 

0.9806 

1.0827 

0.9850 

1.0869 

0.8863 

1.0265 

0.9748 

1.0944 

0.8769 

0.9976 

0.9831 

1.0903 

0.8752 

1.0250 

0.9761 

1.0916 

0.8794 

1.0038 

TABLE  59.  MEAN  X-RAY  INTENSITY  RATIOS  FOR  CURED  PBAA 
PROPELLANT  PARTICLE  SIZE  ANALYSES  (CONSTANT 
CONCENTRATIONS) 


Batch 

Ferric 

Oxide 

R1 

Ammonium 

Perchlorate 

R2 

PBAA 

Polymer 

R3 

Aluminum 

R4 

1 

1.0078 

0.9605 

0.7254 

1.1730 

2 

1.0044 

0.9676 

0.7132 

1.1710 

3 

1.0348 

0.9551 

0.7198 

1.2333 

4 

0.9768 

0.9772 

0.7016 

1.1006 

5 

1.0462 

1.0382 

0.9794 

1.0862 

6 

0.9798 

1 . 0908 

0.8794 

1.0132 

The  Intenelty-particle  size  (weight  percent  fine-fraction)  results 
for  cured  propellant  analysis  are  shown  graphically  in  Figures  17  through 
22.  The  data  can  be  interpreted  in  the  same  manner  as  was  done  for  the 
uncured  propellants.  The  most  evident  difference  between  uncured  and 
cured  propellant  results,  when  the  graphs  are  compared  to  the  same 
scale.  Is  the  lover  sensitivity  attainable  with  cured  propellant  analysis. 
Small  particle  size  changes  can  be  detected  with  a higher  degree  of 
precision  and  accuracy  when  analyzing  uncured  propellants. 

C.  High-Rate  HTPB  Propellant  Results 

The  low-rate  PBAA  propellant  experiments  were  not 
repeated  with  high-rate  propellants,  because  the  principles  involved  and 
the  qualitative  intensity-particle  size  relationships  for  the  two  types 
of  propellants  are  very  similar.  Consequently  the  experimental  ,.nd 
regression  procedures  described  for  the  low-rate  propellant  can  be 
applied  equally  wall  to  tha  high-rata  propellant.  The  major  differ- 
ence between  the  two  types  of  propellants  from  a particle  size  analysis 
standpoint  is  the  ultraflne  ammonium  perchlorate  used  In  the  hlgh-rate 
propellant.  This  ultraflne  ammonium  perchlorate  normally  has  a weight 
median  diameter  of  0.5  to  1 urn. 

It  has  been  demonstrated  that  when  ultraflne  ammonium  perchlorate 
is  processed  in  composite  propellants  of  the  type  described  here  It  forms 
agglomerates  in  the  prope' lant  (20).  Normally,  the  weight  mean 
diameter  of  the  agglomerates  is  in  the  range  of  5 to  10  pm.  The 
agglomerates  can  affect  propellant  properties  and  the  precision  and 
accuracy  of  propellant  spectrometrlc  analyses.  The  X-ray  spectrometr ic 
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TABLE  60.  MULTIPLE  REGRESSION  DATA  FOR  CURED  PBAA 
PROPELLANT  PARTICLE  SIZE  ANALYSES  (CONSTANT 
CONCENTRATIONS) 
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TABLE  62.  EQUATIONS  FOR  ESTIMATING  AMMONIUM  PERCHLORATE  SIZE 
FRACTIONS  (W2)  IN  CURED  PBAA  PROPELLANTS  (CONSTANT 

CONCENTRATIONS) 


•>  -2 

W2  x 10 

= d20  + d21Rl  H 

h d22R2  + d23R3 

+ d24R4 

Intensity 

Ratios 

Coefficient 

Level 

Coefficient 

RMSE 

d20 

16.52225 

0.13675 

R2 

d22 

-8.32892 

R4 

d24 

-6.91237 

d20 

7.95476 

0.00430 

R1 

d21 

-7.29558 

R2 

d22 

-0.61157 

R3 

d23 

0.52355 

d20 

7.10588 

0.00544 

R1 

d21 

-5.61735 

R2 

d22 

-0.34525 

R4 

d24 

-0.61735 

d20 

4.26451 

0.01507 

R2 

d22 

0.54618 

R3 

d23 

-1.75242 

R4 

d24 

-2.68393 

d20 

6.00542 

0.00877 

R1 

d21 

-3.44175 

R3 

d23 

-0.67871 

R4 

d24 

-1.41777 
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TABLE  63.  EQUATIONS  FOR  ESTIMATING  ALUMINUM  SIZE  FRACTIONS 

(W. ) IN  CURED  PBAA  PROPELLANTS  (CONSTANT  CONCENTRATIONS) 

4 


W.  X lo 
4 

' " d40  + d41Rl  ' 

f d42R2  + d43R3 

♦ d44R4 

Intensity 

Ratios 

Coefficient 

Level 

Coef  f icient 

RMSE 

d40 

25.85048 

0.19339 

R2 

d42 

-17.88178 

R4 

d44 

-6.49235 

d40 

9.47954 

0.00392 

R1 

d41 

-1.55671 

R2 

d42 

-5.70588 

R3 

d43 

-1.96760 

d40 

12.66980 

0.01208 

R1 

d41 

-7.86383 

R2 

d42 

-6.70677 

R4 

d44 

2.32029 

d40 

8.69213 

0.00558 

R2 

d42 

-5.45885 

R3 

d43 

-2.45324 

R4 

d44 

-0.57269 

d40 

-8.70761 

0.07225 

R1 

d41 

34.39916 

R3 

d43 

-13.18458 

R4 

d44 

-13.22761 
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TABLE  64.  ESTIMATED  SIZE  FRACTIONS  OF  AMMONIUM  PERCHLORATE  AND  ALUMINUM  IN  CALIBRATION 
BATCHES  OF  CURED  PBAA  PROPELLANT  (CONSTANT  CONCEN1RA LIONS) 
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WEIGHT  % Of  30-mm  AMMONIUM  PERCHLORATE.  W, 

Figure  17.  Chlorine  KQ  Intensity  ee  e function  of  aaaonlua  perchlor- 
ate particle  else  frectlons  for  constant  alualmai  particle  site 
(cured  PBAA  propellant). 
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WCX2HT  % 2*hm>  AMMONIUM  KNCMLOfUTC.W2 


Figure  18.  Alwlnun  Intenalty  u a function  of  amonlun  porchlor 

•to  particle  ilu  fractions  for  constant  alimlnun  particle  olss 
(cured  PBAA  propellant). 
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WEIGHT  % t-nm  ALUMINUM,  W^ 
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Figure  19.  Alualma  Kq  Intensity  r s e function  of  alualnua  particle 

else  frectlons  for  constant  anonlua  perchlorate  particle  site 
(cured  PIAA  propellant). 
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WEIGHT  tt  O-Min  ALUMINUM.  W4 

Figure  20.  Chlorine  Kq  Intensity  ss  a function  of  elunlnum  particle 

else  fractions  for  constant  SMeonlun  perchlorate  particle  size 
(cured  PBAA  propellent). 
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WEIGHT  % 2G-m»  AMMONIUM  PERCHLORATE,  Wj 

Figure  21.  Iron  intaneity  *«  • function  of  mmoqIub  perchlorate 

particle  else  frectione  for  conetent  alueinun  particle  else 
(cured  PMA  propellent). 


20  WEIGHT  % 20-^m  AMMONIUM  FEACHlORATf.  W,  *° 

Figure  22.  Sulfur  K intensity  «a  * function  of  aaBonlun  perchlorate 

particle  else  fraction*  for  constant  eluninua  particle  size  (cured 
FBAA  propellant). 
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technique  described  here  is  the  only  known  one  for  detecting  and  measuring 
ultrafine  ammonium  perchlorate  agglomeration  in  high-rate  propellants. 
Accordingly,  it  is  a very  valuable  tool  for  monitoring  and  controlling 
the  compositions  of  these  propellants  in  a manufacturing  situation. 

The  influence  of  ultrafine  ammonium  perchlorate  agglomeration  on 
X-ray  spectrometric  particle  size  measurements  is  illustrated  in 
Figure  23.  For  this  example  the  recommended  AIK^  radiation  as  well  as 

ClK^  radiation  were  used  to  detect  the  ammonium  perchlorate  particle  size 

change  quantitatively.  The  propellant  contained  a bimodal  blend  of 

nominal  1-  and  70-pm  ammonium  perchlorate.  The  ratio  of  this  blend  was 

purposely  varied,  but  all  ingredient  percentages  were  held  constant. 

The  plot  for  C1K  radiation  is  linear  whereas  that  for  AIK  radiation 
a a 

is  curvilinear.  The  curvilinear  relationship  is  caused  by  the  influence 

of  increasing  ultrafine  ammonium  perchlorate  agglomeration,  as  the 

ultrafine  ammonium  perchlorate  percentage  is  increased,  on  the  AIK^ 

intensity.  The  AIK^  radiation  with  its  short  effective  penetration  depth 

and  high  sensitivity  for  small  ultrafine  ammonium  perchlorate  particles 

detects  the  agglomeration  effect  as  a departure  from  linearity.  The 

C1K  radiation,  on  the  other  hand,  is  insensitive  to  the  ultrafine 
u 

ammonium  perchlorate  agglomerates,  but  is  sensitive  to  the  large  change 
in  weight  mean  diameter  of  the  ammonium  perchlorate  resulting  from 
varying  the  bimodal  blend  ratio.  Consequently,  the  relationship  using 

C1K  radiation  is  linear.  The  relationship  with  AIK  radiation  would 

a a 

also  be  linear  if  the  percentage  of  ultrafine  ammonium  perchlorate 
agglomeration  remained  constant  as  the  total  ultrafine  ammonium  per- 
chlorate percentage  increased.  In  any  event  the  AIK^  intensity  will 

increase  as  the  ultrafine  ammonium  perchlorste  agglomerate  size  increases. 
Thus,  the  X-rsy  fluorescence  method  will  readily  detect  an  abnormal 
ultrafine  ammonium  perchlorate  agglomeration  condition  in  high  rate 
propellants,  and  thereby  prevent  substandard  propellant  from  being  used. 

Ultrafine  ammonium  perchlorate  agglomeration  affects  the  precision 
of  X-rsy  spectrometric  analyses  of  high  rste  propellants  alluded  to  in 
Section  IV. B.  The  AIK^  intensity  measurement  precision  is  strongly 

affected  because  of  its  sensitivity  to  the  agglomeration  and  the  fact 
that  the  agglomeration  creates  a locally  inhomogeneous  condition  in  the 
analyzed  propellant  surface.  The  CIK^  measurement  precision  is  normally 

not  affected,  as  indicated  in  Figure  23,  unless  the  agglomeration  condi- 
tion is  very  severe.  This  was  the  case  for  the  data  in  Tables  46,  47, 
and  48.  High-rate  propellant  processing  procedures  hsve  now  been 
improved  to  the  point  where  C1Kq  intensity  measurements  can  be  made  with 

a high  degree  of  precision,  but  the  ulcrsfine  ammonium  perchlorate 
agglomeration  and  its  influence  on  AIK  intensity  measurements  still 
persists. 
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WEIGHT  % 1-mmULTRAFINC  AMMONIUM  PERCHLORATE 

Figure  23.  X-ray  fluorascenca  analysis  of  aomonlum  perchlorate 
particle  else  in  high  rate  propellant. 


VI.  SIMULTANEOUS  DETERMINATION  OF  INGREDIENT 
CONCENTRATIONS  AND  PARTICLE  SIZES 

A.  Low  Rate  PBAA  Propellant 

1.  Experimental 

It  would  be  desirable  In  the  X-ray  spectrometric 
analysis  of  propellants  to  be  able  to  determine  not  only  Ingredient 
percentages  and  particle  size  variations  but  also  a simultaneous 
variation  of  these  parameters.  Such  a simultaneous  variation  might  be 
expected  to  occur  during  an  actual  propellant  manufacturing  situation. 
This  part  of  the  project  was  done,  therefore,  to  demonstrate  the 
methodology  required.  No  attempt  was  made  to  analyze  unknown  propellants 
In  which  Ingredient  percentages  and  particle  sizes  were  simultaneously 
varied,  but  extension  of  the  principles  to  this  application  will  be 
apparent.  The  primary  limitation  in  this  case  is  the  total  number  of 
variables  that  can  be  determined.  A.s  explained  earlier,  the  total 
number  is  equal  to  the  number  of  different  ingredients  (emission  lines) 
that  can  be  determined.  Consequently,  it  is  often  necessary  to  have 
separate  and  independent  information  about  the  propellant  composition 
to  supplement  the  X-ray  fluorescence  analysis  and  provide  a complete 
determination. 

A 1/8  fraction  of  a 2^  factorial  design  was  chosen  for  the  prepar- 
ation of  calibration  batches  required  for  the  multiple  regression 
analysis.  The  factors  and  factor  levels  are  given  in  Table  65,  and 
the  design  data  and  defining  contracts  are  given  in  Table  66.  The 
particle  size  of  ammonium  perchlorate  was  varied  by  changing  the  ratio 
of  a bimodal  blend  of  nominal  20-  and  200-pm  size  fractions;  whereas 
the  particle  size  of  aluminum  was  varied  by  changing  the  blend  ratio 
of  nominal  9-  and  32-ym  sizes.  The  actual  compositions  of  the 
resulting  nine  calibration  batches,  including  the  midpoint  of  the 
design,  are  shown  in  Table  67.  The  instrumental  parameters  for  t*~is 
determination  are  recorded  in  Table  68. 

2.  Uncured  Propellant  Results 

Four  replicate  samples  of  each  calibration  batch 
were  analyzed.  The  individual  X-ray  intensity  ratios  for  each  of  the 
four  ingredients  are  listed  in  Table  69.  The  mean  intensity  ratios 
used  in  the  regression  analysis  are  listed  in  Table  70. 

Multiple  regression  data  for  uncured  propellant  analysis  with 
terms  for  all  ingredient  concentrations  and  particle  sizes  present 
are  shown  in  Table  71.  The  relatively  small  standard  errors  and  large 
correlation  indices  indicate  thst  the  model  is  a good  estimate  of  the 
true  relationship.  The  best  sets  of  multiple  regression  equations  are 
shown  in  Table  72.  Except  for  the  ferric  oxide  determination,  the 
particle  size  variables,  and  W^,  are  in  all  sets.  The  aluminum 

particle  size  as  previously  shown,  has  very  little  effect  on  iron 
radiation;  hence,  this  variable  is  not  in  the  best  set  for  ferric  oxide. 
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TABLE  65.  FACTORS’  AND  FACTOR  LEVELS  FOR  CALIBRATION  BATCHES 
(VARIABLE  CONCENTRATIONS  AND  PARTICLE  SIZES) 


Fac  tor 

Symbo 1 

Low  Level  (g) 

High  Level  (g) 

Ferric  Oxide 

A 

4.5 

5.5 

Ammonium  Perchlorate 

B 

660.0 

700.0 

PBAA  Polymer 

C 

125.0 

145.0 

Aluminum 

D 

150.0 

170.0 

Ammonium  Perchlorate 
Particle  Size 

E 

*30.01?  - 20  urn 
70.0%  - 200  pm 

40.0%  - 20  pm 
60.0%  - 200  pm 

Aluminum 
Particle  Size 

F 

*80%  - 9 pm 
20%  - 32  pm 

100%  - 9 pm 
0%  - 32  pm 

*Weight  percent  on  a total  ingredient  basis. 


An  analysis  of  variance  for  these  results  is  given  in  Table  73. 

This  analysis  is  similar  to  that  given  in  Section  IV. A,  and  can  be 
interpreted  in  the  same  manner.  The  residual  errors  in  this  experiment 
are  significant  at  the  5%  significance  level  for  ferric  oxide,  ammonium 
perchlorate,  and  PBAA  polymer  determinations. 

3.  Cured  Propellant  Results 

Individual  X-ray  intensity  ratios  for  cured 
propellant  determinations  are  recorded  in  Table  74;  the  mean  values 
are  recorded  in  Table  75.  The  multiple  regression  equations  with  ail 
independent  variables  present  in  the  model  are  shown  in  Table  76.  The 
results  compare  favorably  with  those  for  uncured  propellant  analysis. 

The  best  sets  of  multiple  regression  equations  for  cured  propellant 
analysis  are  shown  in  Table  77.  Terms  for  aluminum  and  ammonium 
perchlorate  particle  size  effects  are  present  in  all  sets.  Like  uncured 
propellant  analyses,  the  best  sets  of  equations  provide  only  a small 
Improvement  in  the  estimation  of  intensity  ratios  when  compared  with 
the  complete  linear  model. 

An  analysis  of  variance  for  the  cured  propellant  results  is  given 
in  Table  78.  The  residual  errors  for  ferric  oxide  and  ammonium  per- 
chlorate determinations  are  significant  at  the  5%  significance  level. 
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Defining  contrasts:  I,  ADE,  BCE,  ACF,  BDF,  ABCD,  ABEF,  CDEF 


TABLE  68.  INSTRUMENTAL  PARAMETERS  FOR  COMBINED  INGREDIENT  PARTICLE 
SIZE  AND  CONCENTRATION  DETERMINATIONS  (PBAA  PROPELLANT) 
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TABLE  70.  MEAN  X-RAY  INTENSITY  RATIOS  FOR  UNCUREO  PBAA 
PROPELLANT  ANALYSES  (VARIABLE  CONCENTRATIONS  AND 
PARTICLE  SIZES) 


Batch 

Ferric 

Oxide 

R1 

Ammonium 

Perchlorate 

R2 

PBAA 

Polymer 

R3 

Aluminum 

R4 

1 

1.1240 

0.8980 

0.8219 

0.9906 

2 

0.9415 

0.8064 

0.9026 

1.1127 

3 

1.0044 

0.8134 

0.9368 

1.0127 

4 

1.1628 

0.9153 

0.9859 

0.8953 

5 

0.9178 

0.8844 

0.8047 

0.9818 

6 

1.1504 

0.7699 

0.9998 

1.0660 

7 

1.2012 

0.7016 

0.8427 

1.1896 

8 

1.0042 

0.7940 

0.9804 

1.0223 

9 

1.0740 

0.8058 

0.9109 

1.0473 

B.  High  Rate  Propellant  Results 

The  composlt Iona  of  the  calibration  batches  for  Exper- 
iment 3 of  the  hlgh-rate  propellant  analyses  are  given  In  Table  79. 

Thes  ■ compositions  in  combination  with  those  In  Tables  40  and  41  form 
a complete  central  composite  design  with  four  variables.  The  fourth 
variable  Is  the  percentage  of  0.5-ym  ammonium  perchlorate  in  the 
propellant. 

J 

The  mean  X-ray  Intensity  for  each  of  the  three  Ingredients,  calcu- 
lated from  duplicate  sample  determinations,  1s  given  in  Table  80.  All 
samples  were  analyzed  using  the  instrumental  conditions  shown  in 
I Table  42. 

A complete  quadratic  polynomial  model,  as  shown  in  Tahle  81,  was 
fitted  to  the  data.  The  model  gives  a good  fit  to  the  data  for  each 
i of  the  three  ingredients,  but  It  cannot  be  inverted  In  a stralghtfor- 

| ward  manner  for  estimating  ingredient  percentages  and  particle  sizes. 

A polynomial  model  was  used  for  this  experiment  because  of  the  rela- 
tively wide  ingredient  concentration  and  particle  size  ranges  used. 

A linear  model  would  not  be  an  accurate  representation  of  the  true 
curvilinear  relationships  Involved. 
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TABLE  71.  MULTIPLE  REGRESSION  DATA  FOR  UNCURED  PBAA 
PROPELLANT  ANALYSES  (VARIABLE  CONCENTRATIONS 
AND  PARTICLE  SIZES) 


Ann* 'ii  i ti 
l vt.  hi. * 


i i'vt  i i ■ i 

> m 

1 l'Vi'1 

h * fi  X f li  X + I'  X.  ♦ h + h . V*  , * 1>  U , 
Hi  ill  i . . i i i i • ■ i i th  • 


>•>•1  ti>  i*  ni 


HI 

I.4VJD0  t.’.DhHJl 

u.omo;  ii.  I uh  v> 

U.nh.M’  II.  It. ’11%  1) , S 

lUhoin  ii.ii.Mh  o.4;shh 

-O  HUSH.’  I).  ’Jill  y*  - i. 

o,  uihj.m  o.  iinn  ■#»  n.ii  m’h 


n. 

*i  * 

0. 

44SHM 

- 1. 

I n. 

-0. 

'.l.'.x 

-n, 

1 h / 

-i). 

111.'. 

*ni  is 

1 1 . 

; hh;  l 

- \ , 

141  ( 

-MHO 

n.  1 1 m 


o.  r>n 

U.n/KKH 
*n.  t H >♦» n J 
-IJ,  l»J*4  » 


H 


l>.  441)  t4 


o ^ I i.  :u.‘, 


ii.iiiji;  I <i.n;.«;  | 


U.ll4h 

(I  HU*»  ) U>.  *. 


1*  >4l) 

u.twu  ih 

(1 . i M i . 

U,4W«.  j 

n,n*  ’in* 

ti.  * 

1 >. *.*t’l J * 

ii.  im^»> 

O.IJH-,  J 

<>.  i .*i 

ti.  i.’H ’4 

1 , l MAIN 

1 1 , 1 K IM  1 S 

i).  mil  * i 

- . 

(1,  IKJ41  7 

U.ixxhf  7 

-4  till 

117 


TABLE  74.  X-RAY  INTENSITY  RATIOS  FOR  CURED  PBAA  PROPELLANT 
ANALYSES  (VARIABLE  CONCENTRATIONS  AND  PARTICLE  SIZES) 


Ferric 

Oxide 

Ammonium 

Perchlorate 

P8AA 

Polymer 

Aluminum 

1.1277 

0.9819 

0.7333 

1.0989 

1.1119 

0.9874 

0.7366 

1.0846 

1.1076 

0.9839 

0.7400 

1.0760 

1.1043 

0.9853 

0.7359 

1.0615 

0.9615 

0.9249 

0.8224 

1.2352 

0.9562 

0.9314 

0.8137 

1.2210 

0.9351 

0.9446 

0.8012 

1 . 1913 

0.9370 

0.9366 

C. 8025 

1.2072 

0.9750 

0.9577 

0.7858 

1.1722 

0.9814 

0.9648 

0.7834 

1.1499 

0.9916 

0.9568 

0. 797d 

1.1529 

0.9659 

0.9620 

0.7861 

1.1426 

1.1647 

0.9511 

0.8845 

1.0853 

1.1689 

0.9480 

0.8813 

1.0910 

1.1506 

0.9493 

0.8918 

1.0838 

1.1524 

0.9468 

0.8914 

1.0852 

0.9000 

0.9572 

0.7117 

1.2052 

C.8949 

0.9572 

0.7202 

1.1959 

0.9097 

0.9633 

0.7024 

1.1818 

0.8995 

0.9686 

0.6987 

1.1607 

1.0913 

0.9374 

0.8327 

1.1864 

1.0907 

0.9313 

0.8405 

1.2016  j 

1.0894 

0.9363 

0.8188 

1.2038  j 

1,0936 

0.9431 

0.8126 

1.1773 

1.1469 

0.9246 

0.6984 

1.3055 

1.1623 

0.9163 

0.7134 

1.3025 

1.1534 

0.9304 

0.7086 

1.3030 

1.1534 

0.9237 

0.7032 

1.2864 

0.9726 

0.9360 

0.8457 

1.0974 

0.9776 

0.9348 

0.8484 

1.0953 

0.9760 

0.9343 

0.8409 

1. 1.107 

0.9676 

0.9398 

0.8318 

1.0977 

1.0511 

0.9457 

0.7764 

1.1744 

1.0541 

0.9407 

0.7918 

1,1842 

1.0414 

0.9493 

0.7720 

1.1618 

1.0434 

0.9551 

0.7772 

1.1460 

TABLE  75.  MEAN  X-RAY  INTENSITY  RATIOS  FOR  CURED  PBAA  PROPELLANT 
ANALYSES  (VARIABLE  CONCENTRATIONS  AND  PARTICLE  SIZES) 


Batch 

Ferric 

Oxide 

R1 

Ammonium 

Perchlorate 

R2 

PBAA 

Polymer 

R3 

Aluminum 

R4 

1 

1.1129 

0.9846 

0.7364 

1.0802 

2 

0.9474 

0.9344 

0.8100 

1.2137 

3 

0.9785 

0.9603 

0.7883 

1.1544 

4 

1.1592 

0.9488 

0.8872 

1.0863 

5 

0.9010 

0.9616 

0.7082 

1.1859 

6 

1. 0°12 

0.9370 

0.8262 

1.1923 

7 

1. 1540 

0.9238 

0.7059 

1.2994 

8 

0.9734 

0.9362 

0.8417 

1.1003 

9 

1.0475 

0.9477 

0.7794 

1.1666 

VII.  PRECISION  AND  ACCURACY 

Precision  data  for  low-rate  PBAA  propellant  analyses  were 
given  in  the  analysis  of  variance  Tables  28,  36,  73,  and  78.  To 
optimize  the  precision  of  the  procedure  for  a particular  application, 
however,  the  analyst  needs  to  assess  the  magnitudes  of  the  various 
sources  of  error.  The  methodology  for  resolving  the  error  variances 
with  specific  applications  to  the  uncured  and  cured  PBAA  propellant 
analyses  reported  here  Is  detailed  in  Appendix  A.  Possible  sources  of 
error  based  on  the  experimental  procedure  described  in  Section  II  are 
as  follows: 

a)  Counting  error. 

b)  Sampling  error  for  replicate  samples. 

c)  Error  between  pairs  (or  groups)  of  samples  analyzed 
at  different  times. 

d)  Error  associated  with  the  reference  standard  measurement 
and  determination  of  X-ray  intensity  ratios, 

e)  Instrumental,  mechanical,  and  electronic  variations. 
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TABLE  76.  MULTIPLE  REGRESSION  DATA  FOR  CURED  PBAA  PROPELLANT 

ANALYSES  (VARIABLE  CONCENTRATIONS  AND  PARTICLE  SIZES') 

R,  1 b,,.  + b X,  b h.  ,X,  + b X + h , X b h.rW  + b,,W, 

i iO  ill  1 21  2 13  3 144  i)2  164 


ItiKtvil  i * *n t 


AmmiMi  i ura 
iYrrlii-iraii- 


t'  HAA 
I’ulvni-r 


A 1 ura 1 mi  ra 


CiH' I f i r I rnt 

1.1‘VC  1 


Ciicf  t i r iun  t 


4.72240 


1.76240 


-0.04725 


-0.03524 


-0.04470 


-0.00188 


0.00076 


-0.85551 


0.02143 


0.02172 


0.00467 


0.00441 


0.00171 


-0.00040 


-0.7  464  5 


0.0  J2  14 


0.012  10 


0.07041 


-0.00/44 


-0.00078 


-0.00144 


6.  11  160 


-0.01612 


-0.06261 


-0.0681' 


0.006  IS 


-0. 004 1 4 


0.00048 


0.00042 


0.0887  4 


0.0600? 


0.06212 


0.06224 


0.04407 


0.0  1468 


0.00046 


0.00028 


0.024  16 


0.02526 


1 .88045 


0.24146 


0.  16157 


0. 14466 


0.14406 


2.0  3524 


0.00042  -0.65238 


0.  548  I 1 


0.  10764 


1.71428 


0.04146  -0.18  0-6 


1 . I1  285 


-5. 14284 


0.0J544  -0.44/40 


-2.5/014 


0.02514  -2.70504 


O.  251  18 


0.00014  I -12. .0588 


0. 01X11  7 


0.0104 1 0.44665 


0.08724  20.20288 


0.05405  -0.80016 


0.06106  -0. 57713 


0.06121  -0.84335 


0.00083  -2.26506 


0.01112  0.90873 


0.00/40  0.44641 


0.004  1 I 0.'t484  7 
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TABLE  77.  BEST  SETS  OF  MULTIPLE  REGRESSION  DATA  FOR  CURED  PBAA 
PROPELLANT  ANALYSES  (VARIABLE  CONCENTRATIONS  AND  PARTICLE 
SIZES) 


K.  = b . X . + b.-X.  + h.  .X.,  + b . X , + b W + b..W. 

i lit  122  131  144  152  16  4 


1 iiKrfd  i k-nt 


Oooff  ik’  ik-ill 
l.rVk‘1 


(Toi*  f 1 1 <■  i k-n  t 

Sb 

t 

t. 36970 

1.78820 

0.067  34 

26. 55479 

-0.01326 

0.00375 

-3. 5 )600 

-0.01950 

0.00489 

-3.9877) 

-0.00206 

0.00069 

-2.98550 

0.00070 

0.00035 

2.00000 

0.08382 

0.01227 

0.00213 

5.76056 

0.00178 

0.0005! 

3.49019 

-0.000)7 

0.00025 

-I .48000 

0.47804 

0.05818 

0.00222 

26.20720 

-0.020)3 

0.0022  ) 

-9.11659 

-0.00072 

0.00040 

-1 .80000 

-0. 00140 

0.00020 

-7.00000 

0.00715  0.90569 


0.00562  0.99586 


TABLE  79.  RESIDUAL  COMPOSITIONS*  FOR  DETERMINING  COMBINED 
INGREDIENT  PARTICLE  SIZES  AND  CONCENTRATIONS  IN  HIGH- 
RATE  HTPB  PROPELLANT  (WEIGHT  7.,  EXPERIMENT-3) 


Batch 

35-ym 

Aluminum 

X1 

Ballistic 

Modifier 

X2 

Ammonium 

Perchlorate** 

X3 

1 

14.00 

8.00 

70.65 

2 

14.00 

8.00 

70.65 

3 

14.00 

8.00 

70.65 

4 

14.00 

4.40 

70.65 

5 

14.00 

11.60 

70.65 

6 

11.60 

8.00 

70.65 

7 

16.40 

8.00 

70.65 

8 

14.00 

8.00 

67.05 

9 

14.00 

8.00 

72.25 

10 

14.00 

8.00 

70.65 

11 

14.00 

8.00 

70.65 

*These  compositions  and  those  In  Tables  40  and  41 
form  a central  composite  experimental  design. 

The  fourth  Independent  variable  (X^)  is  weight  2 
0.5-yia  ammonium  perchlorate. 

**Bimodal  blend  of  602  - 0.5-ym  and  402  - 90-pra 
ammonium  perchlorate. 

NOTE:  Remainder  of  propellant  Is  HTPB  binder. 


TABLE  80.  MEAN  X-RAY  INTENSITY  RATIOS  FOR  UNCURED  HIGH  RATE 
1ITPB  PROPELLANT  ANALYSES  (EXPERIMENT-3) 


Batch 

Aluminum 

R1 

Ballistic 

Modifier 

*2 

Ammonium 

Perchlorate 

R3 

1 

0.187 

1.425 

0.784 

2 

0.179 

1.417 

0.779 

3 

0.193 

1.412 

0.795 

4 

0.211 

0.804 

0.830 

5 

0.191 

2.272 

0.767 

6 

0.162 

1.448 

0.799 

7 

0.202 

1.491 

0.763 

8 

0.278 

1.502 

0.733 

9 

0.170 

1.446 

0.804 

10 

0,544 

1.858 

0.658 

11 

0.144 

1.242 

0.815 

The  theoretical  counting  error  for  an  individual  emission  line 
peak  intensity  measurement  is  fixed  by  the  total  number  of  counts 
collected.  A tabulation  of  theoretical  counting  errors  for  various 
fixed  counts  is  given  in  Table  82.  Thus,  the  analyst  can  control  the 
counting  error.  The  resolution  of  the  counting  error  from  the  other 
sources  of  error  in  these  experiments  is  not  straightforward,  as  shown 
in  Appendix  A,  because  of  the  manner  in  which  the  reference  standard 
and  propellant  samples  were  analyzed.  The  measurement  scheme  that  will 
be  used  and  the  counting  error  that  can  be  tolerated  will  depend  on 
the  relative  magnitudes  of  the  counting  error  and  other  sources  of 
error  as  well  as  the  amount  of  time  that  can  be  allotted  to  the 
measurement. 

Although  a detailed  discussion  of  the  various  sources  of  error  for 
FRAA  propellant  analysis  is  given  in  Appendix  A,  it  is  worthwhile  to 
compare  lust  the  sampling  errors  for  low-  and  high-rate  propellant 
analyses  as  shown  in  Table  83.  The  sampling  errors  in  Table  83  are  the 
estimated  relative  standard  deviations  for  individual  sample  analyses. 
The  sampling  error  includes  the  counting  error,  propellant  ’.nhoraogeneity . 
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TABLE  81.  (Concluded) 


-0.00006  0.00001  -4.2808 


TABLE  82.  X-RAY  COUNTING  STATISTICS-RELATIVE 
STANDARD  DEVIATION,  lOO/Vn” 


No.  of 
Counts  (n) 

a 

95% 

Conf idence 
Limit 

99% 

Confidence 

Limit 

3 a 

100 

10.0 

19.6 

25.8 

30.0 

200 

7.07 

13.9 

18.2 

21.2 

500 

4.47 

8.76 

11.5 

13.4 

1000 

3.16 

6.19 

8.15 

9.48 

2000 

2.24 

4.39 

5.78 

6.72 

5000 

1.41 

2.76 

3.64 

4.23 

10,000 

1.00 

1.96 

2.58 

3.00 

20,000 

0.707 

1.39 

1.82 

2.12 

50,000 

0.447 

0.876 

1.15 

1.34 

100,000 

0.316 

0.619 

0.815 

0.948 

200,000 

0.224 

0.439 

0.578 

0.672 

500,000 

0.141 

0.276 

0.364 

0.423 

1,000,000 

0.100 

0.196 

0.258 

0.300 

2,000,000 

0.0707 

0.139 

0.182 

0.212 

5,000,000 

0.0447 

0.0876 

0.115 

0.134 

10,000,000 

0.0316 

0.0619 

0.0815 

U.0948 

•ffacta,  and  error  caused  by  ehorc-term  instrumental  variations.  The 
error  associated  with  measurement  of  the  reference  standard  is  not 
included.  The  counting  error  in  the  table  is  the  theoretical  value 
from  Table  82. 

Except  for  the  determination  of  PBM  polymer  in  cured  propellant, 
the  sampling  errors  for  all  ingredients  in  the  low-rate  propellant  are 
less  than  1Z.  This  shows  that  the  propellant  was  very  homogeneously 
mixed  and  that  the  sample  preparation  procedure  used  provided  for  very 
repeatable  sample  analyses.  The  repeatability  of  ammonium  perchlorate 
and  aluminum  determinations  was  not  as  good  for  two  ultrahigh  burning 


TABLE  83.  ESTIMATED  SAMPLING  ERRORS  FOR  PBAA  AND  HTPB  PROPELLANT  ANALYSES 
~ BY  X-RAY  SPECTROMETRY  (RESTIVE  STANDARD  DEVIATIONS) 


rate  propellants  tested.  This  is  attributed  to  the  ultrafine  ammonium 
perchlorate  agglomeration  problem  that  was  discussed  in  Section  II. B. 
However,  the  ammonium  perchlorate  sampling  error  for  high-rate  HTPB 
propellant  is  very  small.  In  fact,  it  is  significantly  better  than  that 
for  the  low-rate  propellant.  This  excellent  repeatability  for  ammonium 
perchlorate  determinations  is  typical  for  high-  and  ultrahigh-rate 
propellants  made  with  current  improved  processes.  The  sampling  errors 
for  aluminum  determinations  in  high-  and  ultrahigh-rate  propellants,  on 
the  other  hand,  are  typically  In  the  range  of  4%  to  5%.  Because  this  is 
for  a single  aluminum  determination,  the  precision  for  analysis  of 
a particular  propellant  batch  can  be  substantially  improved  by  making 
several  sample  replications  and  averaging  the  results.  The  lower  pre- 
cision for  high-rate  propellant  aluminum  determin.it ions  is  due  to  the 
high  sensitivity  of  AlK  radiation  to  the  inhomogeneous  ammonium  per- 
chlorate agglomerates. 

The  true  accuracy  of  the  X-ray  fluorescence  method  was  not  directly 
evaluated  by  analyzing  primary  propellant  standards.  Nevertheless,  based 
on  extensive  experience  with  the  calibration  procedure  used,  which  pro- 
vides an  accurate  representation  of  the  true  intensity-concentration 
relationships,  and  the  residual  errors  obtained  (Tables  25,  38,  57,  and 
64),  it  is  expected  that  low-rate  propellant  ingndients  can  be  deter- 
mined with  a relative  accuracy  of  IX  or  better.  The  accuracy  of  high- 
rate  propellant  determinations  will  be  somewhat  less,  particularly  for 
aluminum,  because  of  the  ultrafine  ammonium  perchlorate  agglomeration 
problem. 

Application  of  the  recommended  X-ray  fluorescence  analysis  method 
to  a specific  propellant,  as  for  example  in  production,  will  result  in 
the  estimation  of  ingredient  percentages  using  an  expression  such  as 
Equation  (9).  The  estimated  percentages  will  differ  from  the  actual 
(nominal)  percentages  because  of  the  experimental  error  of  the  method. 
Consequently,  the  analyst  must  determine  whether  there  is  a high 
probauility  that  the  estimated  ingredient  percentages  could  have  been 
generated  from  a propellant  formulation  having  the  nominal  or  expected 
composition.  This  determination  can  be  made  by  placing  Joint  confidence 
interval  estimates  on  the  actual  ingredient  percentages  (concentrations). 
The  methodology  for  doing  this  is  developed  and  Illustrated  for  uncured 
PBAA  propellant  analyses  in  Appendix  B. 


VIII.  CONCLUSIONS 

The  X-ray  fluorescence  spectrooetrie  method  described  in  this 
report  has  bean  demonstrated  to  be  an  excellent  procedure  for  monitoring 
and  determining  the  compositions  of  solid  composite  propellants.  The 
recommended  method  Is  nondestructive,  rapid,  and  capable  of  a high  degree 
of  precision  and  accuracy.  Typically,  the  composition  of  a composite 
propellant  batch  can  be  determined  within  15  ic  30  min,  thus  enabling 
substandard  batches  either  to  be  discarded  or  corrected  prior  to  casting 
into  motors.  The  method  applies  equally  v*  11  to  cured  and  uncured 
propellant  samples. 
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The  X-ray  fluorescence  method  is  unique  in  its  ability  to  determine 
ingredient  percentages  and  in-situ  particle  sizes  as  well  as  combina- 
tions of  these  parameters.  The  estimated  relative  standard  deviation 
for  the  determination  of  ingredients  in  low-rate  propellants  and 
ammonium  perchlorate  in  high-rate  propellants  is  IX  or  less.  The 
estimated  relative  standard  deviation  for  determining  aluminum  in 
high-rate  propellants  is  larger  (4Z  to  5X)  because  of  problems  caused  by 
ultrnfine  ammonium  perchlorate  agglomeration.  The  ability  of  the  method 
to  detect  and  quantify  the  ultrafine  ammonium  perchlorate  agglomeration 
is  a very  valuable  feature  for  high  burning  rate  propellant  applications. 
Because  of  the  influence  of  ultrafine  ammonium  perchlorate  agglomeration  on 
propellant  properties,  this  capability  for  agglomeration  analysis  alone 
justifies  the  application  of  the  method  tohigh-rate  propellant  manufacture. 
With  four  to  six  sampie  replications,  the  accuracy  for  determining  all 
ingredient  percentages  is  of  the  order  of  1%  to  27«  relative  or  better. 


The  X-ray  fluorescence  method  does  have  some  practical  limitations 
with  respect  to  propellant  ’.nalysis.  It  is  an  elemental  emission  method 
capable  of  detecting,  from  a practical  standpoint,  elements  of  atomic 
No.  12  (magnesium)  and  higher.  Therefore,  organic  propellant  ingredients 
such  as  hydroxyl-terminated  polymers,  plasticizers,  and  some  ballistic 
modifiers  cannot  be  detected.  Consequently,  a complete  analysis  requires 
that  these  ingredients  be  determined  by  an  alternative  procedure.  It 
should  be  emphasized,  however,  that  a substandard  propellant  composition 
can  be  detected  by  the  X-ray  method  if  the  percentage  of  a detectable 
element  changes,  even  if  the  change  is  caused  by  an  error  in  the  per- 
centage of  a nondetectable  ingredient. 


IX.  RECOMMENOEO  IMPLEMENTATION 

The  X-ray  fluorescence  method  described  here  has  potential 
application  to  all  types  of  solid  composite  propellants  used  in  Army 
missile  systems.  The  method  is  ideai  for  application  to  propeii  nt 
production  because  of  its  speed,  precision,  and  nondestructive  nature, 
VIPER  propellant,  which  is  in  engineering  development,  is  an  exceilent 
candidate  at  this  time,  because  of  the  anticipated  large  production 
rate  and  the  high  propellant  cost.  The  method  has  been  recommended  to 
the  VITER  Project  Office.  Their  personnel  have  been  apprised  of  Its 
salient  features.  The  method  should  also  be  considered  for  application 
to  the  propellant  manufacturing  processes  for  the  PATRIOT  and  PERSHING 
missile  systems. 

This  MTT  project  should  logically  progress  to  a Manufacturing 
Methods  and  Technology  (MM&T)  project  to  demonstrate  specific  applica- 
tion of  the  developed  test  method  to  a propellant  manufacturing  process. 
The  MM&T  project  would  perhaps  best  be  conducted  by  the  appropriate 
missile  system  propulsion  subcontractor.  Such  a project  should  include 
the  purchase,  installation,  and  demonstration  of  automated,  computer 
controlled  instrumentation  that  would  be  more  suitahle  for  a production 
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application  than  the  manual  instrumentation  used  in  this  MTT  project. 
Th*'  necessary  automated  instrumentation  is  commercially  available  at 
a cost  of  approximately  $100,000. 

Finally,  the  X-ray  fluorescence  test  procedure  will  be  prepared  in 
the  proper  format  and  submitted  for  possihle  inclusion  in  MIL-STD-286B 
as  a tentative  method.  MIL-STD-286B  currently  contains  no  test  methods 
for  finished  composite  piopell.mts. 
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Appendix  A.  RESOLUTION  OF  ERROR  VARIANCES 


Section  IV  clearly  outlined  that  the  experimental  technique  involved 
taking  two  pairs  of  observations  at  each  concentration  (batch)  of  the 
mixture  for  hoth  uncured  and  cured  PBAA  propellant  samples,  one  pair 
from  the  first  sample  and  another  pair  from  the  second  sample.  An 
important  point  is  that  within  a sample,  the  measurement  in  seconds  on 
the  reference  standard  is  taken  only  once.  Hence  the  numerator  of 
both  ratios  will  be  the  same.  However,  the  reference  standard  is 
measured  again  for  the  observations  of  the  second  sample.  The  intensity 
data  in  seconds  were  given  in  Table  17.  The  four  ratios  for  the  uncured 
samples  from  which  the  average  ratios  in  Tahle  19  were  calculated  are 
shown  in  Tahle  18.  The  first  number  in  each  block  in  Table  17  is  the 
reading  in  seconds  for  the  reference  standard  and  the  next  two  the 
readings  for  the  unknown  in  question,  all  for  the  first  sample  from  that 
batch.  Directly  below  the  first  sample  measurement  are  three  similar  measure- 
ments for  the  second  sample  from  the  same,  batch.  Tables  29  and  30  give  the  sec- 
onds and  intensity  ratios,  respectively,  for  the  cured  samples . The  number 
of  counts  varied  from  ingredient  to  ingredient . These  are  shown  in  Table  16 . 

It  is  of  interest  to  the  experimenter  to  know  the  relative  orders 
of  magnitude  of  the  different  sources  of  variance  in  the  experiment. 

It  is  expected  that  the  different  batches  would  cause  the  greatest 
variation  because  the  concentrations  are  changing  across  batches.  The 
other  sources  are  the  sample  to  sample  variation  and  within-sample 
variation.  Theoretically,  the  latter  consists  of  two  components: 

(1)  variance  due  to  the  counting  error,  and  (2)  propellant  Inhomogeneity 
and  other  errors  — electronic,  etc.  These  error  variances  shall  be 
2 2 

called  as  and  a , respectively. 

1 . Count ing  Error 

The  number  of  radiation  counts  used  In  an  experiment  will  be 
called  n.  I t Is  generally  assumed  that  the  probability  model  describing 
the  numher  of  counts  produced  by  the  counting  device  in  time  t is  the 
Poisson  distribution  [231  with  probability  function 


(x  « 0,1,2,...) 


(A-l) 


where  X is  the  counting  rate  in  counts  per  unit  time.  Parrish  (26)  in 
reviewing  the  general  problem  of  counting  error  says  "Two  measurements 
of  a constant  Intensity  in  which  counting  is  performed  during  equal 
times  t,  will  not  in  general  yield  the  same  number  n of  counts,  owing 
to  a random  distribution."  He  mentions  approximating  this  with  a 
Gaussian  distribution  with  mean  N and  standard  deviation  o ■ kR,  where 
N is  an  average  of  the  number  of  counts  la  time  t obtained  from  a large 
number  of  experiments.  This  corresponds  to  the  normal  approximation 
to  the  Poisson  distribution  for  the  number  of  counts  In  a given  time  t 
for  the  fixed  time  procedi  ).  The  procedure  used  here  is  the  fixed 
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count  procedure.  The  measurement  used  here  is  the  ratio  of  the  time 

t that  it  takes  any  n counts  to  be  taken  on  the  standard  to  t , the 
s s u 

time  that  it  takes  for  n^  counts  to  be  made  on  the  component  of  unknown 

concentration.  Thus,  if  the  counting  variance  is  to  be  separated  from 
the  other  variations  in  the  experiment,  the  variance  of  this  ratio  is 
required. 

Either  the  numerator  or  denominator  of  that  ratio,  i.e.,  the  time 
t for  n counts  is  considered.  If  the  number  of  counts  in  time  t follows 
the  Poisson  law  as  given  by  Equation  (A-l),  and  if  x denotes  the  random 
variable  representing  the  time  to  the  nth  count,  then 


P (X  < 0)  = 0 

r 


(A- 2) 


and 


P (0  < X < t) 


I 


k — X t 
(At)  6 
k! 


(A- 3) 


The  first  follows  because  this  time  cannot  be  negative.  The  second 
equation  merely  says  that  the  probability  of  getting  n or  more  counts 
in  time  t is  the  probability  of  requiring  t or  less  time  for  n counts. 
Thus,  Equations  (A-2)  and  (A-3)  describe  the  distribution  function  for 
this  random  variable  X.  Equations  (A-2)  and  (A-3)  can  then  be  differ- 
entiated with  respect  to  t in  order  to  obtain  the  probability  density 
function. 

vO  CD  09 

1 V (At)1*  -At  -At  V (At)k~lj  A -At  V (At)k-A 

Jt  Z k!  6 " e L (k-1) ! " e L,  k! 

k-n  k*n  k-n 


P(t) 


, nn-1  -At 

A t e 

r(n) 


(A— 4 ) 


This  describes  the  well-known  Gamma' density  function  with  parameters 
n-1  and  1/A,  Thus  the  waiting  time  for  n counts  is  a Gamma  distribution. 
The  distribution  of  t /t_  ia  required  where  t.  and  t_  are  both  Gamma 
variables. 

At  the  outset,  it  will  be  assumed  that  the  measurements  involve 
different  counting  rates  A^  and  and  different  numbers  of  counts 
n^  and  n^.  Thus, 
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p (t L>  “ 


n -1  -X 

(tL)  e 

•'  (tij  j 


1C1 


P(t2)  = 


n2  n2  1 _X2t2 
*2  (t2^  e 

r(n2) 


(t1  > 0;  n1,X1  > 0) 


(t2  > 0;  n2,X2  > 0) 


If  the  independence  of  t^  and  t„  is  assumed, 


P(t1,t2) 


x\;2  v1  v1  -<Yi+w 

1 2 t1  t e 

r(ni)r(n2) 


Letting  u = t /t  , the  required  ratio  and  z = t„,  the  following  is 
obtained : 


"D(t  ,t  ) " 

P(u,z)  - p(tltt2)  j • (A-5) 


After  evaluating  the  Jacobian  and  simplifying,  Equation  (A-5)  becomes 
n -1  n +n  -1  -(Az  u + A-z) 

p(u,z)  = C u z e (A-6) 


where 


C 


r(ni)r(n2) 


The  density  function  for  u is  required,  so  if  z is  Integrated  out,  t he 
following  is  obtained 


p(u) 


nl~l  f nj+n2-l  -(A^zu+A2z) 

C u I z e dz 


0 


n*  n«  n. -1 

XL  L\2  u r(nL  + n2) 


r(n1)r(n2)(x1u  + x2) 


nl+n2 


(A-7) 
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The  moments  of  this  distribution  can  then  be  found  in  the  usual  way. 
The  mean  is  as  follows: 


E (u)  = E 


X1  n2  ' 1 


If  n^  and  n^  are  both  large,  this  quantity  reduces  to  Aj/A^,  the  ratio 

of  the  counting  rates.  To  obtain  the  variance  of  this  ratio,  first  the 
crude  second  moment  must  be  found: 


E(u  ) = K 


at 

I 


n +1 


0 (A^u  + A^) 


nl+n2 


du 


where 


K >=  C-T(n1  + n2). 


E(u2) 


ac 

I 


n +1 


I n,+n„ 

*(■•¥) 


du 


<»2> 


nl+n2 


If  A^u/A2  ■ X,  then  the  preceding  integral  becomes 


E(u2) 


(a2) 


/a  \1+*  ? ni+l 

K _/^2  \ t X 1 
n +n9  \A  / / n 

1 2 X 17  *0  (1  + X)  1 


dx 


(A-8) 


After  integrating  und  simplifying,  Equation  (A-8)  becomes 


2.  X2  (ni  + 1)nl 
t(u  ) - -r  


x2  (n2  - l)(n2  - 2) 


t 

$■ 

n 
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For  the  variance  of  u the  following  is  obtained: 


2 

a u 


F.(u2) 


[K(u)]2 


(n1  + Dnj 

(n2  - 1) (n2  - 2) 


(n2  - l)(n1  + l)(n1)  - n2(n2  - 2) 

(n2  - l)2(n2  - 2) 


n1(n1  - 1)  + n1n2 
(n2  - l)2(n2  - 2) 


(A— 9 ) 


If  n^  = n2  = n,  i.e.,  if  the  same  number  of  counts  is  used  for  the 

standard  as  for  the  unknown,  and  if  it  is  further  assumed  that  n is 
large.  Equation  (A-9)  will  reduce  to 


2 „ 
o u - 


(A-10) 


Equations  (A-9)  and  (A-10)  give  the  variance  of  the  ratio  t^/t2  under 

the  conditions  specified.  It  is  interesting  to  note  from  Equation  (A-9) 
that  if  the  experimenter  is  willing  to  allow  n^  j*  n2>  he  can  reduce  the 

variance  by  making  n2  *>n^.  Thus  it  appears  that  in  terms  of  counting 

precision  It  is  best  to  use  a larger  number  of  counts  on  the  component 

of  unknown  concentration,  rather  that  splitting  up  a large  number  of 

counts  equally  among  the  standard  and  unknown.  For  example,  if 

n^  * n2  ■ 10,000  is  used,  u2u  A2/A2  whereas  if  n^  * 1000  and 

n2  - 20,000  is  used,  o2u  - A2/*2  0.0000025. 

2.  1. Inear  Model 


It  was  necessary  to  arrive  at  a method  of  using  the  data  in 
Tables  17,  18,  29,  and  30,  along  with  the  theoretical  information  obtained 
In  this  appendix  on  the  counting  variance  to  make  an  overall  evaluation 
of  the  error  variances.  If  it  is  assumed  that  the  ratio  data  (for  the 
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cured  and  uncured  propellant  samples)  for  the  Kth  Ingredient  follows 
the  model 


+ S 


00 


+ C 


00 


j(i>  ljl 


+ c 


(k) 

ljl 


where 


fkl  2 (kl 

B;  = batch  effect,  with  variance  a 
1 P 

(k)  2(k) 

S , , , . = effect  of  the  sample  within  batches  with  variance  a 
j(i)  K s 

(k)  2(k) 

C...  = counting  error  with  variance  o 
ljl  c 


(k) 

e i = within  sample  error  not  including  counting  error,  e.g., 

J inhomogeneity  and  sample  instrumental  error,  with  variance 
2(k) 


It  is  particularly  important  to  determine  what  portion  of  the  within- 
sample  variation  is  due  to  counting  variance  and  what  portion  is  due 
to  the  remaining  errors  because  the  experimenter  can  control  the 
counting  variance.  The  ratio  data  in  Tables  18  and  30  were  analyzed 
as  a hierarchial  [27]  (subsampling)  classification  not  in  order  to 
make  any  particular  tests  but  to  estimate  the  variance  components. 
Tables  A-l  and  A-2  show  the  mean  squares  for  the  uncured  and  the  cured 
propellant  data,  respectively. 

Certain  linear  combinations  of  the  mean  squares  are  unbiased 
estimates  of  the  variance  components.  It  is  reasonably  easy  to  show 
that  the  following  expressions  represent  the  expected  mean  squares. 

The  k superscript  will  be  dropped  at  this  stage. 


E(hse)  ■ o\  + 


(MS~  Is  error  mean  square) 

fc. 


E<HSB(A)>  ■ 

E (MS , ) - o2 
A e 


2 2 2 

o + o + 2o 
c c s 

2 2 2 

+ o + 2o  + 4oD 

C 8 8 


(MS_ ^ ^ is  mean  square  samples 
B(A; 

within  batches) 

(MS  is  mean  squares  batches) 
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TABLE  A-l.  VARIANCE  ANALYSIS*  FOR  UNCURED  PBAA  PROPELLANT 


i 


I 


* t 


Source 

DF 

SS 

MS 

Ferric  Oxide 

Batch 

11 

0.40205711 

0.03655064 

Samples 

12 

0.00196559 

0.00016371 

Error 

24 

0.00141607 

0.00005900 

Ammonium  Perchlorate 

Batch 

11 

0. 05074630 

0.00461330 

Samples 

12 

0.00112890 

0.00009407 

Error 

24 

0.00108068 

0.00004502 

PBAA  Polymer 

Batch 

11 

0.17377193 

0.01579744 

Samples 

12 

0.00099994 

0.00008332 

Error 

24 

0.00095275 

0.00003969 

Aluminum 

Batch 

11 

0.12441220 

0.01131020 

Samples 

12 

0.00159362 

0.00013280 

Error 

24 

0.00196859 

0.00008202 

*This  is  a nested  analysis  in  which 
"samples"  are  actually  "samples  within 
batches"  and  error  represents  the 
variation  "between  observations  within 
samples,"  and  batches. 

DF  * Degrees-of-freedom 
SS  ■ Sum  of  squares 
MS  ■ Mean  square 


Therefore,  the  estimates  are  given  by  the  following: 


2 2 

MS„  estimates  a 4-  a 

E c c 


MSB(A)  " 2 

— estimates  o 

2 s 


msa  - MSB(A) 


estimates  o 


B 


(A-ll) 
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TABLE  A-2.  VARTANCH  /\NALYS1S*  FOR  CURED  PBAA  PROPELLENT  SAMPLES 


Source 

DF 

SS 

MS 

Ferric  Oxide 

Batch 

9 

0.32633644 

0.03629200 

Samples 

10 

0.00187111 

0.00018011 

Error 

20 

0.00147753 

0.00007388 

Ammonium  Perchlorate 

Batch 

9 

0.01225774 

0.00136197 

Samples 

10 

0.00050185 

0.00005018 

Error 

20 

0.00068033 

0.00003401 

PBAA  Polymer 

Batch 

9 

0.13362281 

0.01484697 

Samples 

10 

0.00194809 

0.00019480 

Error 

20 

0.00196067 

0.00009803 

Aluminum 

Batch 

9 

0.12455100 

0.01383900 

Samples 

10 

0.00189946 

0.00018994 

Error 

20 

0.00138341 

0.00006917 

*Thls  Is  a nested  analysis  In  which  "samples" 
are  actually  "samples  within  batches"  and 
error  represents  the  variation  "between 
observations  within  samples",  and  batches. 

DF  “ Decrees  of  freedom 
SS  ■ Sura  of  squares 
MS  - Mean  square 


The  results  In  Tables  A-l  and  A-2  can  then  he  used  to  arrive  at  estimates 
of  these  variance  components  for  each  of  the  Ingredients  for  both  the 
uncured  and  cured  propellant.  However,  an  elaboration  must  be  made 

2 

concerning  o , the  actual  theoretical  counting  variance.  Using 

c -) 

Equation  (A- 10),  a*-  can  be  computed  because  in  all  cases  n^  - n7  * n 

with  n being  greater  than  20,000.  MS_  represents  the  vlthln-sample 

t 

variation  In  the  data,  l.c.,  the  variation  between  two  replicates  within 
a sample.  However,  it  was  emphasized  earlier  that  within  a sample  the 
same  measurement  on  the  reference  standard  was  used  in  the  numerator 
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of  hoth  ratios.  Therefore,  while  the  theoretical  counting  variance 
derived  in  Paragraph  1 is  the  true  one,  it  is  not  the  one  which  is 
being  estimated  by  MS^  because  the  latter  does  not  take  into  account 

counting  error  associated  with  the  standard.  Thus,  to  be  able  to 

2(k) 

isolate  an  estimate  of  o , an  alteration  in  the  counting  variance 
must  he  made  to  account  for  the  experimental  procedure  used  here. 


The  variable  v = l/t^  is  considered  where  t^  follows  a Gamma  distri- 
bution. The  distribution,  mean,  and  variance  of  this  random  variable 
will  be  found.  For  expedience  the  n2  and  notation  will  not  be  used 

as  before.  Instead,  the  usual  notation  for  a Gamma  density,  namely, 
parameters  u and  (1  will  be  used.  Therefore,  a = n^  - 1 and  1/A  = I . 

It  follows  that 


P(t2) 


r(u+l)Ba+1 


R(v) 


1 -l/vt> 
- e 
v 

u+1 


It 


r(j+i)f 

The  value  of  the  l.icohian  Is  1/v' 

i+2 


,,  W 

K(v)  - 


- 1 / r v 


T(a+1)B 


u+1 


(A-12) 


Thus  g(v)  can  be  simplified  to 


The  first  two  moments  of  this  distribution  will  now  he  found: 


K(v> 


Lt? 

/ 


-u-1  -1/  v 
v e dv 


r<a+l)6 


a+i 


'(u+1) 6 


at* 


. -a-1  - 1 / fiv 

(6v)  e 

-i-l 


dv 


If  2 ■ l/6v,  then  the  preceding  expression  hecomes 


K(v) 


1 

I- ( u+1)  6 


a> 


a- 

z 


1 -z 
e 


dz 


r(a)  L 

F(a+1)6  “ 6a 
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2.  1 I -ct  -l/8v  . 

h(v  ) = "VVT  I v e dv 

r(a+i>e  1 Q 

-a  -1/Bv 


/ 


I’(n+])B 


1 f (Bv)"°te~ 


dv 


(A-l 3) 


As  before,  z = 1/Bv.  Equation  (A-13)  then  becomes 


E(v  > ' r(„+i) 


W ZVZ  I'  b?) 


dz 


= r(g-i>  „ 


r(<+l)B2  32a(«-l) 


j2v  = E(v2)  - [e(v)1  2 = — 


(A-14) 


Substituting  n2  - a + 1 and  Aj  - l/B  into  Equation  (A-14)  the  following 
is  obtained 

,2 


a v ■ — 


(n2  - l)~(n2  - 2) 


(A-l 5) 


(A-16) 


for  the  case  where  n2  = n and  n is  large.  Equation  (A-lb)  gives  an 
expression  for  the  variance  of  l/tj.  Actually  Var(t^/t2)  is  needed 

2 2 3 

where  t^  is  considered  to  be  a constant;  Var(tj/t;)  ■ tj(X2/n  ).  On 
the  average  for  n^  • n,  ■ n,  it  is  expected  that  t^  ■ n/A^.  Thus,  an 
approximation  to  the  required  variance  is: 

2 

v 

V 


'nr  r2 


(A-l 7) 


It  must  be  emphasized  here  that  this  is  not  the  counting  variance.  The 
latter  was  given  earlier  in  Paragraph  1.  This  is,  however,  the  variance 
which  is  estimated  by  MS^  in  these  experiments. 
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3. 


Estimation  of  Variance  Components 


The  expressions  in  Equation  (A— 11)  were  used  to  compute  the 
hatch  and  sample  variance  components.  Equation  (A-10)  was  used  to 

2 

compute  the  true  counting  variance  a"  and  Equation  (A— 17)  was  used  to 


2 2 2 
•ompute  o'  ; was  computed  hy  MS^.  - o^  . 


The  occurrence  of  and  in  formulae  for  the  variance 

presents  a hit  of  a problem.  Stri  tly  speaking,  these  A*s  are  not  a 
qualitative  function  of  the  ingredient  hut  a quantitative  one.  That 
is,  the  A's  depend  on  the  concentration  of  each  ingredient.  The  X's 
were  estimated  for  a particular  ingredient  by  taking  the  counting  rate 
for  that  ingredient  averaged  over  all  the  data.  The  X^'s  were  found 

in  the  same  way.  Table  A- 3 gives  the  estimated  variance  components  for 
each  ingredient  for  the  cured  and  uncured  propellants.  Also  included 
2 

is  a",  the  true  counting  variance. 

c 

2 

As  expected,  a *7  is  always  large.  The  estimate  of  within-sample 

^ *2 

variance  (excluding  counting)  cj  is  approximately  the  same  for  eacli 

component  for  both  uncured  and  cured  propellant  with  the  exception  of 
PBAA  polymer  in  the  uncured  propellant  where  this  variance  seems  to  be 

'2 

exceptionally  small.  Otherwise  the  order  of  magnitude  of  n"  seems  to 

-5  L 

be  in  the  vicinity  of  3,7  * 10  . In  comparing  this  variance  with  the 

theoretical  counting  variance  (which  is  the  proper  comparison  to  make, 

2 2 

i.e,,  o should  be  compared  with  a ),  it  is  noted  that  for  n in  the 
*.  c 

2 2 

range  of  20,000  to  30,000  a and  > _ seem  to  he  ahoiit  the  same  order  of 

magnitude.  That  is,  for  n In  this  range  the  within-sample  variance  is 
approximately  50 Z due  to  counting  and  50Z  due  to  other  errors.  An 

? 

increase  from  50,000  to  500,000  decreases  a"  hy  a factor  of  10.  As 

c 

previously  mentioned,  it  is  expected  that  the  counting  variance  is  he  t 
reduced,  not  by  increasing  the  total  number  of  counts  but  by  taking 
more  counts  on  the  unknown  and  fewer  on  the  standard.  It  must  he 
emphasized  though  that  the  overall  variance  within  samples  is  small 
and  it  is  doubtful  that  reassigning  n.  and  n would  reduce  the  overall 

I *• 

2 2 

variance  within  samples,  l.e.,  by  more  than  a factor  ot  two. 

Of  course  it  is  assumed  that  all  of  this  reduction  is  accounted  for  by 
the  decrease  of  0 . 
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TABLE  A-3.  VARIANCE  COMPONENTS  AND  ESTIMATES  OF  VARIANCE  COMPONENTS 
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Appendix  B.  CONFIDENCE  INTERVAL  ESTIMATE  OF  CONCENTRATION: 
EVALUATION  OF  PROPELLANT  PRODUCTION 

The  purpose  of  this  appendix  is  to  develop  and  illustrate  a method 
for  determining  how  well  the  point  estimates  described  in  Section  IV 
estimate  the  concentrations.  The  obvious  way  of  attempting  thi”  type 
of  determination  is  via  the  route  of  joint  confidence  interval  estimates 
on  the  actual  concentrations.  There  are  certain  nominal  concentrations 
that  the  production  process  will  attempt  to  attain.  It  renains  then  to 
determine  if  there  is  strong  reason  (i.e.,  high  probability)  to  believe 
that  the  intensity  ratios  from  the  sample  nr  the  estimates  of  the  con- 
centrations from  Equations  (15)  could  have  been  generated  from  mixtures 
with  the  nominal  concentrations.  This  procedure  represents  a "go"  or 
"no  go"  type  of  situation,  i.e.,  based  on  the  confidence  intervals,  it 
is  concluded  that  the  true  concentrations  either  are  or  are  not  what 
they  are  supposed  to  be.  The  underlying  theory  on  which  these  confi- 
dence intervals  are  based  is  found  in  the  following  paragraphs.  The 
theory  and  development  is  general.  The  application  was  made  here  only 
to  uncared  PBAA  propellant  samples  where  ingredient  particle  sizes  were 
held  constant. 


Box  and  Hunter  [28]  discuss  the  problem  of  finding  Joint  confidence 
interval  estimates  on  the  solution  of  a set  of  simultaneous  equations 
when  the  coefficients  are  subject  to  error.  Their  work,  was  actually 
part  of  a more  specific  problem  of  finding  a confidence  region  for  a 
stationary  point  in  response  surface  analysis.  However,  the  theory 
can  be  applied  to  the  problem  discussed  here  because  the  estimates  of 
the  concentrations  are  found  by  solving  a set  of  simultaneous 
regression  equations. 


It  is  supposed  that  there  are  m simultaneous  equations  of  the 


type : 


m 

I 

j-0 


bijxj 


(i  - 1,2,. ...m) 


( B— 1 ) 


where  the  b^  are  subject  to  error.  The  quantities 


m 

I 

J-0 


Vj 


i 


(i  - 1,2,. . . ,m) 


(B-2) 


are  considered  where  the  £ are  the  values  of  the  X's  that  would  satisfy 
Equation  (B-l)  if  the  actual  regression  coefficients  were  used  in  place 
of  the  b^.  For  this  work,  the  t's  represent  the  actual  concentrations; 

thus  Equation  (B-2)  is 


♦ 
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where 


to 


2 

j=i 


Vj 


Then  It  is  desirable  to  attach  joint  confidence  intervals  on 

?2 » * and  • 

If  a vector  of  the  6's  is  considered,  say  _6,  as  having  a multi- 
variate normal  distribution  with  mean  vector  <3  and  variance-covariance 

matrix  E(6_  6^)  = V,  then  6.'  V ^ _6  follows  a X*'  distribution  with  m 
degrees-of-f reedom  (m  = 4 in  this  case).  The  remarks  made  here  rely 
on  the  assumption  that  the  c in  Equation  (7)  follow  a normal  dis- 

2 

tribution  with  zero  mean  and  some  variance  a . For  estimates  of  the 
elements  of  V,  the  following  is  obtained: 


Var  (R1  - R^)  = 


1 + 


Cov  (Rj  ~ Rj»  \ ~ \)  = s 


2 2 Wil-  sii-" 

n + 2 2 Wl  j 

h 1 J 


h 1 
1 + 


Sik'H 


(B-3) 


(B-4) 


where  s is  the  sample  estimate  of  the  error  variance  and  s^ 


is 


the  sample  estimate  of  the  covariance  between  the  e.,  and  the  c,  , 
^ ij  kj 


(j  * l,2,...,n)  in  Equation  (7), 


hi 


is  the  (hi)  element  of  the 


inverse  of  the  matrix  of  corrected  sums  of  squares  and  cross  products 
of  the  y's  for  the  calibration  sample.  Replacing  V by  its  estimate, 
and  dividing  by  the  appropriate  degrees-of-f reedom,  the  following  ratio 
is  obtained: 


n - 8 


22 

1 k 


5lSkv 


ik 


ik 


which  is  distributed  as  F with  4 and  n-8  degrees-of-f reedom.  Here  w' 
is  the  (ik)  element  of  the  inverse  of  the  matrix  W,  the  matrix  of 
residual  sum  of  squares  and  cross  products  of  the  R’s  from  the  original 
data.  The  following  is  obtained: 


2 Vj  - 2 vj  ■ ri  • \ 

J J 


(B-5) 
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If  6^  is  replaced  in  the  preceding  F statistic  by  the  expression  given 
in  Equation  (B-5) , the  following  is  obtained: 


F4,n-8 


n - 8 


4 


I II I <v vai - v bijbki“Ik 

i j k 1 _ 


H 

II  <V  S'  (ii  - h>  "ji 

J 1 

H 


(B-6) 


where 


Ml 


is  the  (jl)  element  of  the  matrix: 


Q = 


B'  W 


B 


B is  the  matrix  of  regression  coefficients  for  the  set  of  4 regression 
lines.  Equation  (B-6)  represents  joint  confidence  intervals  on  the 
actual  concentrations.  That  is,  given  values  of  the  estimates  X^,  X0, 

X^,  and  X^ , particular  values  of  Cj.  Cj*  ^3'  ar*d  ^4  cnn  substituted 

into  Equation  (B-6)  and  if  the  resulting  expression  is  less  than 
F . a (upper  tall  point)  then  the  C's  fall  inside  the  100  (1-a)  2 

01  * 4 ^ Il“o  ^ 

confidence  region.  The  W and  Q matrices  for  the  uncured  PBAA 
propellant  data  are: 


w-1  - 

’ 7214.8162 

2554.8459 

4014.0983 

-3201.5790 

-1714.2325 

2679.7650 

3439.8046" 

2122.4663 

-1867.4456 

- 

2825.4942. 

Q - 

'24274.424 

-15.343 

2.4899 

-274.4115 

2.8058 

10.8662 

219.582 

0.77389 

-36781 

5.3264  . 

The  obvious  f/  s to  consider  are  the  nominal  concentrations.  If 
the  X's  are  close  to  the  f.'s,  then  this  F value  in  Equation  (B-6)  will 
be  dost  to  zero,  the  nominal  values  would  be  covered  by,  say  a 952 
confidence  region  and  hence  it  would  be  concluded  that  the  mixture  con- 
centrations do  not  deviate  from  the  target  values.  In  case  the 
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est hum <’s  deviate*  significantly  from  the  target,  i.e.,  if  the  F value 
is  significant,  the  target  T.'s  fall  outside  the  95%  confidence  region 
and  it  is  concluded  that  the  mixture  composition  is  not  what  It  should 
he.  It  can  be  said  that  there  is  0.05  probability  that  the  set  of 
estimates  calculated  could  have  heen  generated,  for  which  the  95%  joint 
confidence  interval  would  not  cover  the  true  concentrations.  Thus, 
there  is  good  protection  against  wrong  conclusions  that  the  process  is 
not  producing  the  desired  concentrations. 

Very  often  the  difficulty  with  this  kind  of  procedure  is  that  the 
width  or  area  covered  by  the  confidence  region  is  very  large.  That  is, 
it  is  often  impractical  because  it  could  always  be  concluded  that  the 
production  process  is  attaining  ingredient  concentrations  that  do  not 
differ  significantly  from  those  desired.  This  means  that  the  sensi- 
tivity of  the  confidence  interval  method  is  low,  or  in  terms  of  hypoth- 
esis testing  language,  the  power  (in  classifying  the  concentrations 
as  not  differing  from  the  nominal  ones  when,  in  fact,  there  are  differ- 
ences) is  very  low.  This  problem  did  not  appear  in  this  work.  The 
confidence  regions  were  narrow  and  the  results  have  much  practical  use. 

The  nominal  concentrations  of  uncured  PHAA  propellant  ingredients 

for  this  example  correspond  to  the  midpoint  of  the  design  (Table  15) 

and  are.*  K,  = 0.5,  r ~ b8.0,  = 13.5,  and  C.  = 16.0.  All  of  these 

1 l i s 

are  in  weight  percent.  This  leaves  2.0%  for  the  remaining  binder 

components.  If  a sample  is  taken  from  a batch  of  propellant  of  unknown 

ingredient  concentrations  and  analyzed  for  each  ingredient  and  i f Equations 

(15)  result  in  the  estimated  concentrations  = 0.5,  = 68.0, 

X^  = 14.5,  and  = 17.0,  then  substituting  into  Equation  (B-6)  yields 

an  F value  of  7.45.  The  numerator  and  denominator  degrees-of-f reedom 

are  both  4 because  n = 12  for  the  original  experiment.  The  upper  95% 
point  is  6.39.  This  means  that  the  nominal  concentrations  are  not 
covered  by  the  95%  confidence  region;  thus  it  is  concluded  that  this 
analyzed  mixture  has  true  concentrations  that  deviate  significantly 
from  the  nominal  ones.  The  estimates  did  not  deviate  a great  deal 
from  the  C's;  yet  the  procedure  was  able  to  detect  the  difference. 

As  another  example,  if  X,  ■ 0.49,  X„  ■ 68.0,  X = 13.5,  and  X,  = 15.5, 

an  F value  of  5.02i  will  result.  This  value  is  less  than  tne  95% 

point;  thus  it  is  concluded  that  the  concentrations  do  not  deviate  a 
significant  amount  from  the  nominal  ones. 

It  should  be  noted  here  that  Equation  (B-6)  should  contain  not 
the  actual  f.'s  and  X'.s,  but  the  corrected  C's  and  X's,  i.e.,  each 
corrected  for  the  average  concentration  of  that  component  in  the 
original  calibration  mixtures  shown  in  Tabic*  15.  These  averages  are: 


H - 0.494 
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C,  = 68.26 
7 = n.37 
= 16.02 


Therefore,  for  the  previous  example,  the  following  should  be  inserted 
into  Equation  (B-6): 

Xx  = 0.49  - 0.494,  X2  = 68.0  - 68.26,  etc. 

= 0.5  - 0.494,  ^2  = 68  - 68.26,  etc. 

As  mentioned  earlier,  the  C's  in  Equation  (B-6)  are  the  elements 
of  the  Inverse  matrix  of  sum  of  squares  and  products  using  the  original 
calibration  mixes.  This  matrix  is 


C 


34.72988 

i. 00642 

1. 13464 

0.83355 

4.29894 

4.89776 

4.38275 

4.20892 

4.29237 

4.59324 

Several  charts  were  prepared  which  Illustrate  the  use  of  this 
method  for  evaluating  a propellant  mix  of  unknown  ingredient  concen- 
trations. Figures  B-l  through  B-7  are  given  in  which,  for  the  C's 
held  at  the  nominal  values,  contours  of  constant  probability  I*  were 
plotted  for  variable  values  of  the  X's,  the  estimates.  The  contours 
represent  constant  (i  - u)  probability  corresponding  to  a confidence 
level  whose  confidence  region  contains  the  nominal  C's  exactly  at 
the  houndarv.  For  example,  for  a set  of  estimates  Xj,  X^,  X^,  and  X^ , 

a value  of  P of,  say,  0.1  means  that  the  nominal  C's  are  at  the 
boundary  of  a 90%  confidence  region.  A good  rule  of  thumb  might  be 
to  consider  a probability  of  0.95  as  being  significant,  i.e.,  if 
p > 0.05,  then  the  deviation  hetween  the  estimates  and  the  nominal 
concentrations  is  considered  to  he  significant. 
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Flgura  B-2 . Probability  contours  for  ^ and  X 
when  X,  - 13. 5X,  k.  - 15. OX. 


Flgura  B-5.  Probability  contours  for  and  X 
whan  X,  - 68.0X,  X,  - 13. 5X. 


Figure  a-7. 

■ 


Probability  contours  for  and  X 
13.07,,  X,  - 16. OZ. 
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